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PREFACE 


This  instruction  report  documents  the  computer  program  CSETT  that  can  be 
used  to  compute  induced  stresses  under  irregular  loads  and  the  resulting  con¬ 
solidation  settlements.  Funds  to  modify  the  computer  program  and  to  write  the 
user's  guide  were  provided  to  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES)  by  the  Civil  Works  Directorate  of  the  Office,  Chief  of  Engi¬ 
neers,  US  Army  (OCE) ,  under  the  Geotechnical  Aspects  of  the  Computer-Aided 
Structural  Engineering  (GCASE)  Project. 

Specifications  for  the  program  were  provided  by  members  of  the  GCASE  Task 
Group : 

Mr.  Phillip  Napolitano,  New  Orleans  District  (Chairman) 

Dr.  Roger  Brown,  South  Atlantic  Division 

Mr.  Larry  Cooley,  Vicksburg  District 

Mr.  Frank  Coppinger,  North  Atlantic  Division 

Mr.  Lavane  D.  Dempsay,  St.  Paul  District 

Mr.  Ed  Demsky,  St.  Louis  District 

Mr.  Charles  Trahan,  LMVD 

Mr.  Thomas  Wolff,  St.  Louis  District 

Mr.  Reed  Mosher,  WES 

Dr.  N.  Radhakrishnan,  WES 

Mr.  Bill  Strohm,  WES 

Mr.  Dennis  Williams,  WES 

Mr.  Richard  Davidson,  OCE 

Mr.  Don  Dressier,  OCE 

Mr.  Rixby  Hardy,  OCE 

CSETT,  by  Mr.  Alexis  E.  Templeton,  US  Army  Engineer  District,  Vicksbutg 
(VXD),  is  a  modified  version  of  a  computer  program  developed  by  Mr.  James  G. 
Flock  and  Dr.  Micheal  E.  Pittman,  New  Orleans  District  (NOD),  during  the 
period  November  1978  to  March  1979.  The  original  program  was  developed  under 
the  direction  of  Mr.  Leonard  Manson,  Chief,  Engineering  Systems  and 
Programming  Section,  NOD. 

This  report  was  written  by  Mr.  Templeton.  The  report  and  program 
modifications  were  made  under  the  guidance  and  supervision  of 
Messrs.  Robert  L.  Fleming,  Chief,  Analytical  Section,  Foundation  and  Materials 
Branch,  and  Larry  A.  Cooley,  Chief,  Foundation  and  Materials  Branch,  VXD. 
Technical  assistance  was  provided  by  Messrs.  Mosher  and  Williams,  Automation 
Technology  (AT)  Center,  WES.  Appendix  B  was  provided  by  Mr.  Flock  and 
Dr.  Pittman.  The  work  was  managed  and  coordinated  by  Dr.  Radhakrishnan, 

Chief,  AT  Center,  WES,  and  CASE  Project  Manager.  Mr.  Hardy  and,  later. 
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Hr.  Davidson,  Geotechnical  Branch,  Civil  Works  Directorate,  were  the  OCE 
points  of  contact. 

Commanders  and  Directors  of  WES  during  the  development  of  the  program  and 
the  publication  of  this  instruction  report  were  Colonel  Nelson  P.  Conover,  CE, 
and  Colonel  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  Frederick  R. 
Brown. 
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Program  Information 


Program  Description 

CSETT,  denoted  10019  in  the  Conversational  Oriented  Real-Time  Programming 
System  (CORPS)  library,  computes  consolidation  settlement  for  two-  and  three- 
dimensional  irregular  loads  using  Terzaghi's  one-dimensional  consolidation 
theory.  Solutions  are  provided  for  ultimate  and  time-rate  settlements.  In¬ 
duced  stresses  are  computed  by  integration  of  either  the  Boussinesq  or  Wester- 
gaard  point  load  equation  for  irregular  loads. 


Operation  Overview 

Data  can  be  entered  at  the  user  terminal  interactively  during  execution  of  the 
program  or  by  a  predefined  data  file.  Interactive  input  follows  a  question- 
answer  sequence.  Data  files  are  created  using  command  words  and  saved  in  per¬ 
manent  files  for  later  use  in  the  program.  The  data  file  format  is  described 
in  WES  report  "Computer  Program  for  Determining  Induced  Stresses  and  Consoli¬ 
dated  Settlements  (CSETT)." 


CoETT  has  interactive  edit  capabilities,  and  numerous  error  checks  are  in¬ 
cluded  to  reduce  input  errors.  Echoprint  of  input  data  is  available.  Output 
consists  of  settlements  at  each  location  and  time  specified.  Output  and  echo- 
prints  are  written  to  the  user's  terminal  or  permanent  file. 


Coding  and  Computer  Systems 

CSETT  is  written  in  FORTRAN  and  is  operational  on  the  following  systems: 
si.  US  Army  Engineer  Waterways  Experiment  Station  (WES)  Honeywell 


DPS/1. 


V 


b.  Control  Data  Corporation  CYBERNET  System,  CYBER  170/760. 

Accessing  CSETT 

To  log-on  and  run  CSETT,  follow  the  initiation  procedures  below.  User's 
responses  are  underlined. 

WES  DPS/1 

064200  Dial  the  computer  and  hit  a  carriage  return  at  the  tone.  The  computer 
will  respond  with  the  following: 

C0EWES  HIS  TIMESHARING  ON  10/17/83  at  10.328  CHANNEL  2244  TS1 

USER  ID  —  ROKAOMP  (Your  user  identification) 

PASSWORD  —  XXXXXX  (Your  password) 

USERS  =  043  TSS=145K  MEM-USED=28  SYS=-128K  PR0=2  OOO-WAIT-OOOK 

*F0RT  NEW  (Places  system  in  proper  subsystem) 

*RUN  WESLIB/C0RPS/I0019,R  (Executes  CSETT) 

CYBERNET  System 


Dial  the  computer  and  hit  a  carriage  return  at  the  tone.  The  computer  will 
respond  with  the  following: 


83/10/17.  11.07/05  AA313IA 


EASTERN  CYBERNET  CENTER  SN487  NOS  1 . 4/531 . 281 /18AD 
FAMILY:  KOE  (Family  name  you  will  be  working  under) 

USER  NAME:  CER000  (Your  user  identification) 

PASSWORD  XXXXXX  (Your  password) 

TERMINAL:  515,  NAMIAF 

RECOVER/CHARGE:  CHARGE , EROXX , CEROOO  (Your  charge  account  and  user 
identification) 

/  OLD,CORPS/UN=CECELB  (Calls  up  CORPS) 

/  BEGIN, .CORPS,  10019  (Executes  CSETT) 


CORPS  Computer  Program  Librarj 


The  CORPS  contains  many  other  useful  programs.  These  may  be  catalogued  from 


CORPS  by  the  LIST  command.  The  execute  command  for  CORPS  on  the  WES  system 


RUN  WESLIB/CORPS/CORPS , R 

ENTER  COMMAND  (HELP, LIST, BRIEF, MESSAGE, EXECUTE, OR  STOP) 
*?LIST  or (BRIEF  or  MESSAGE,  etc.) 

On  the  CYBERNET  System,  the  commands  are: 

OLD, CORPS /UN=CECELB 

ENTER  COMMAND  (HELP, LIST, BRIEF, MESSAGE, EXECUTE  OR  STOP) 
*?LIST  or  (BRIEF  or  MESSAGE,  etc.) 


3 


♦ 


CONTENTS 


PREFACE  . 

CONVERSION  FACTORS,  INCH-POUND  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT  . 

PART  I:  INTRODUCTION  . 

Purpose  of  Program  CSETT  . 

Organization  of  Report  . 

PART  II:  THEORETICAL  CONSIDERATIONS  . 

Stress  Induction  . 

Ultimate  Settlement  . 

Time-Rate  of  Consolidation  . 

PART  III:  MECHANICS  OF  SOLUTION  PROCEDURE  .  .  .  . 

Introduction  . 

Stratification  . 

In  Situ  Properties  . 

Load  Definition  . 

PART  IV:  INPUT  GUIDE  . 

Source  of  Input  . 

Data  Format  . 

Data  Entry  from  Terminal  ....  . 

Data  Entry  from  a  File . . . 

Sections  of  Input  . 

Units  . 

Input  Description  . 

Abbreviated  Input  Guide  . 

PART  V:  EXAMPLE  PROBLEMS  . 

Example  Problem  1  . 

Example  Problem  2  . 

Example  Problem  3  . 

Example  Problem  4  . 

Example  Problem  5  . 

REFERENCES  .  . 

APPENDIX  A:  EXAMPLE  RUNS  WITH  HAND  VERIFICATIONS  . 

APPENDIX  B:  GENERAL  EQUATIONS  FOR  STRESS  INDUCTION 

C 


3 


CONVERSION  FACTORS,  INCH-POUND  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Inch-pound  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

feet 

kips  (1000  lb  force) 
kips  (force)  per  square  foot 
pounds  (force)  per  square  foot 
pounds  (mass)  per  cubic  foot 
square  feet  per  year 


By 

0.3048 

4.448222 

47.880263 

47.880263 

16.018463 

0.0939 


_ To  Obtain _ 

meters 

kilonewtons 

kilopascals 

pascals 

kilograms  per  cubic  meter 
square  meters  per  year 
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USER'S  GUIDE:  COMPUTER  PROGRAM  FOR  DETERMINING 
INDUCED  STRESSES  AND  CONSOLIDATION  SETTLEMENT 
(CSETT)* 

PART  I:  INTRODUCTION 


Purpose  of  Program  CSETT 


1.  The  computer  program  CSETT  is  designed  to  compute  induced  stresses 
under  general-shaped  loads  and  to  evaluate  resulting  consolidation  settlements 
in  underlying  clay  strata.  The  program  is  capable  of  modeling  the  following 
loading  and  foundation  conditions. 

«i.  Loads  defined  in  two  or  three  dimensions. 

Jj.  Ramp  loads  and  instantaneous  loads. 

c^.  Single  or  double  drainage  conditions. 

d^  Stratified  compressible  and  incompressible  materials. 

e^.  Normally  consolidated  or  overconsolidated  deposits. 

f_.  Soil  rebound  due  to  excavations  occurring  prior  to  loading. 


Organization  of  Report 


2.  The  remainder  of  this  report  is  divided  into  the  following  parts: 

_a.  Part  II  presents  a  brief  description  of  the  theory  and  assump¬ 
tions  utilized  in  the  program. 

_b.  Part  III  describes  the  mechanics  of  the  solution  procedure. 

c_.  Part  IV  presents  instructions  for  the  input  and  execution  of  the 
program. 

jd.  Part  V  presents  example  problems  to  illustrate  data  input 
and  program  capabilities. 

e.  Appendix  A  presents  example  solutions  with  hand  verif ications. 

*  Three  pages  entitled  "Program  Information"  have  been  hand-inserted  inside 
the  front  cover  of  this  report.  They  present  general  information  on  the 
program  and  describe  how  it  can  be  accessed.  If  procedures  used  to  access 
this  and  other  CORPS  library  programs  should  change,  recipients  of  this 
report  will  be  furnished  a  revised  version  of  the  "Program  Information." 


Appendix  B  describes  the  derivation  of  the  general  equations  for 
vertical  stress  induction  for  the  two-  and  three-dimensional 
load  cases. 


)  ’ 
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PART  II:  THEORETICAL  CONSIDERATIONS 


3.  The  theory  and  assumptions  employed  in  the  program  CSETT,  including 
stress  induction,  ultimate  settlements,  and  time-rate  of  settlement,  are  dis¬ 
cussed  below. 


Stress  Induction 


4.  Stress  induction  equations  used  in  CSETT  for  computing  vertical 
stresses  due  to  imposed  loads  were  derived  from  the  Boussinesq  and  Westergaard 
point  load  formulae  integrated  over  general-shaped  regions.  These  equations 
are  applicable  to  any  two-dimensional  load  defined  by  straight  line  segments 
and  for  three-dimensional  rectangular  or  irregular  loads  of  uniform  pressure 
(Flock  and  Pittman  1979). 

5.  The  Boussinesq  and  Westergaard  stress  induction  equations  are  based  on 
the  theory  of  elasticity.  The  assumptions  used  in  the  development  of  these 
equations  are  as  follows: 

a.  The  soil  is  a  semi-infinite  half  space. 

b.  The  soil  is  elastic  and  homogeneous  and  obeys  Hooke's  law. 

c.  The  soil  is  weightless. 

d.  The  soil  volume  remains  constant. 

e.  There  is  stress  continuity. 

f.  Stress  distribution  is  symmetrical  about  the  vertical  axis. 

£.  The  loads  are  perfectly  flexible. 

h.  The  loads  are  applied  to  a  horizontal  surface. 

i.  The  loads  maintain  a  constant  shape. 

In  addition,  the  Boussinesq  theory  assumes  that  the  soil  is  isotropic,  and  the 
Westergaard  theory  assumes  there  are  no  lateral  deflections. 

6.  In  most  cases  the  assumptions  are  not  so  restrictive  as  to  prohibit 
the  use  of  these  theories.  However,  the  assumption  of  a  perfectly  flexible 
load  is  valid  only  for  earth  structures  and  flexible  foundations  and  is  not 
valid  for  structures  with  rigid  foundations.  Therefore,  when  analyzing  rigid 
structures,  the  engineer  should  give  particular  attention  to  determining  the 
applicability  of  these  theories. 

7.  Flock  and  Pittman's  discussion  of  the  development  of  the  stress 


induction  equations  derived  from  Boussinesq  and  Westergaard  formulae  is  pre¬ 
sented  in  Appendix  B. 


Ultimate  Settlement 


8.  The  ultimate  settlement  for  a  soil  layer  of  a  given  thickness  and 
loading  condition  is  dependent  on  the  stress-volume  relationship  of  the  soil. 
CSETT  allows  the  user  to  define  this  relationship  either  in  terms  of  pressure 
and  void  ratio  or  in  terms  of  pressure  and  strain.*  In  addition,  the  stress- 
volume  data  may  be  input  by  means  of  a  piecewise  linear  curve  or  a  bilinear 
curve  as  described  below. 

£.  The  piecewise  linear  curve  is  defined  by  a  series  of  stress- 
strain  points  over  the  range  of  stresses  to  be  considered  in 
the  analysis.  The  stress-strain  behavior  is  assumed  to  vary 
linearly  between  defined  points. 

jj.  The  bilinear  curve  is  defined  by  a  compression  slope  (C^)  for 
stresses  greater  than  the  in  situ  stress  and  a  recompression 
slope  (C  )  for  stresses  less  than  the  in  situ  stress. 

For  each  method  a  rebound-recompression  slope  (CR)  defines  the  behavior  of  the 
material  during  unloading  and  recompression  (i.e.,  at  stresses  less  than  the 
in  situ  stress).  See  Figure  I. 

9.  Assuming  the  stress-strain  curve  represents  the  field  response  of 
the  soil  to  loading  and  unloading,  the  user  can  determine  the  change  in  strain 
resulting  from  an  increased  effective  stress  due  to  loading  directly  from  the 
curve  by  taking  the  difference  between  the  strain  associated  with  in  situ 
stress  and  the  strain  resulting  from  the  increase  in  effective  stress  due  tc 
loading  as  shown  below: 

change  in  strain  (At)  =  in  situ  strain  -  final  strain  (1) 

10.  The  ultimate  settlement  AH  is  then  defined  as 

AH  =  — ^ —  H  (2a) 

1  +  e 

o 


*  Throughout  this  report  the  pressure-void  ratio  or  pressure-strain  relation¬ 
ship  will  be  referred  to  as  the  stress-strain  relationship. 


or 


AH  =  AeH 


(2b) 


where 

Ae  =  change  in  void  ratio 

e  *  in  situ  void  ratio 
o 

Ae  =  change  in  strain 
H  *  thickness  of  compressible  stratum 

11.  The  stress-strain  relationships  are  generally  based  on  laboratory 
compression  tests.  Conventional  interpretation  of  these  tests  considers  only 
primary  consolidation  (that  compression  occurring  during  which  time  the  rate 
of  compression  is  controlled  by  the  resistance  to  flow  of  water  out  of  the 
soil  voids)  and  neglects  secondary  compression  and  creep  (compression  occur¬ 
ring  after  excess  hydrostatic  pressures  have  dissipated).  CSETT  considers 
only  primary  consolidation  and  is  not  recommended  for  analyzing  settlements  in 
soils  which  exhibit  significant  secondary  compression. 

Time-Rate  of  Consolidation 


12.  The  basic  equations  used  in  time-rate  of  consolidation  computations 
are  derived  from  Terzaghi's  classical  theory  of  one-dimensional  consolidation 
(Taylor  1948).  The  following  assumptions  are  used  in  the  derivation: 

a.  Homogeneous  soil. 

b.  Saturated  material. 

c.  Negligible  compression  of  soil  grains  and  water. 

d.  Continuous  behavior  of  the  soil. 

e.  One-dimensional  compression. 

f.  One-dimensional  flow. 

The  validity  of  Darcy's  law. 

h.  Constant  soil  properties. 

i.  Constant  rate  of  change  of  void  ratio  with  pressure. 

The  time-rate  of  consolidation  theory  provides  a  means  of  evaluating  at  any 
time  the  degree  of  consolidation  (U)  that  has  occurred  within  the 
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consolidating  strata.  CSETT  computes  time-rate  solutions  for  instantaneous 
loading  and  time-dependent  loading  assuming  constant  initial  excess 
hydrostatic  pressure. 

13.  For  the  case  of  instantaneous  loading  the  degree  of  consolidation 
at  any  time  is: 


m=°° 

U  =  1  -  M2  EXP(-M2T) 
m=0 


(3) 


where 

w  _  if  (2m  +  1) 

M - ^ - 

m  =  0,  1,  2,  3,  ... 

c  t 

T  =  time  factor  *  — x 
H 

c  =  coefficient  of  consolidation 
v 

t  =  time  at  which  U  is  desired 
H  =  length  of  longest  drainage  path 

14.  For  time-dependent  loading  (Figure  2) ,  the  degree  of  consolidation 
at  any  time  is: 


For  T  <  T 

c 


00 

u-fc) 

m*0 

1  -  EXP(-M2T)J 

For  T  i  T 

c 


(4) 


/  X  ”  .  V  I 

] 

£eXP(M2Tc)  -  1 

J  exp(-m2t)I 

€ 


u  -  1 


(5) 


L 


where 


c  t 
v  c 


t  =  time  at  end  of  loading  (Figure  2) 

All  other  terms  are  as  defined  previously. 

Further  discussions  on  the  degree  of  consolidation  for  instantaneous  and  time- 
dependent  loading  may  be  found  in  Taylor  (1948)  and  Olson.* 

15.  Any  sequential  combination  of  instantaneous  and/or  time-dependent 
loads  may  be  defined  using  the  two-dimensional  soil  load,  two-dimensional  pres¬ 
sure  load  or  three-dimensional  rectangular  load  options.  When  applying  the 
three-dimensional  irregular  load  option,  the  user  assumes  all  loads  to  be 
applied  as  specified  for  the  first  load  defined  (see  "Load  Definitions," 

pp  19-25) . 

16.  The  length  of  the  longest  drainage  path,  H  ,  is  determined  for 
each  stratum  according  to  the  drainage  condition  specified  for  that  particular 
stratum  (Figure  3) . 

a.  If  the  stratum  is  designated  double  drained,  H  equals  1/2  the 
stratum  thickness. 

b.  If  the  stratum  is  designated  single  drained,  H  is  taken  as 
the  total  stratum  thickness. 

c.  If  a  stratum  is  designated  contiguous,  H  is  found  by  combin¬ 
ing  the  thicknesses  of  the  contiguous  strata  and  treating  them 
as  single  or  double  drained  as  specified  for  the  upper  stratum 
in  the  group  of  contiguous  strata. 

17.  It  should  be  emphasized  that  the  analyses  of  this  program  are  based 
on  the  assumptions  of  one-dimensional  compression  and  one-dimensional  drain¬ 
age.  For  loading  cases  that  will  result  in  significant  horizontal  displace¬ 
ments  and/or  horizontal  drainage  the  settlements  and  time-rate  of  settlement 
will  be  underestimated  by  CSETT. 


*  R.  E.  Olson.  1977.  Unpublished  Notes,  "Consolidation  under  Time- 
Dependent  Loading,"  University  of  Texas,  Austin,  Texas. 
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PART  III:  MECHANICS  OF  SOLUTION  PROCEDURE 


Introduction 


18.  The  mechanics  of  the  solution  procedure  employed  in  CSETT  are  de¬ 
scribed  below.  The  user  should  be  familiar  with  available  options  to  ensure 
that  loading  and  consolidation  conditions  are  modeled  correctly. 

Stratification 

19.  Each  stratum  must  be  designated  as  compressible  or  incompressible. 
Compressible  strata  are  designated  as  having  single  drainage  or  double  drain¬ 
age  or  being  contiguous  with  adjacent  strata.  If  a  stratum  is  designated  as 
contiguous,  its  consolidation  characteristics  (i.e.,  length  of  drainage  path) 
are  determined  as  if  it  were  part  of  the  stratum  with  which  it  is  contiguous. 

20.  Each  compressible  stratum  is  internally  subdivided  into  layers 

3  ft*  thick  or  less  to  increase  the  accuracy  of  settlement  calculations.  The 
original  stratification  as  supplied  by  the  user  is  retained  for  output 
purposes . 

In  Situ  Properties 

21.  The  stress-strain  relationship  supplied  by  the  user  is  utilized  in 
computing  the  in  situ  properties  of  each  stratum.**  The  procedure  varies  de¬ 
pending  on  the  type  of  stress-strain  relationship  supplied  (see  paragraph  8) 
and  whether  an  excavation  precedes  loading. 

22.  If  the  stress-strain  relationship  is  defined  by  a  piecewise  linear 
curve  and  no  excavation  precedes  loading,  the  in  situ  properties  are  deter¬ 
mined  as  follows: 

a.  The  in  situ  overburden  stress  (PqN)  is  found  by  summing  the 
stress  contribution  of  each  material  from  the  ground  surface  to 
the  midheight  of  the  stratum  of  interest  (Figure  4). 

b.  The  in  situ  strain  is  determined  by  interpolating  between  ap¬ 
propriate  points  on  the  stress-strain  curve  (Figure  3) . 


*  A  table  of  factors  for  converting  inch-pound  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 

**  Unless  noted  otherwise  "stratum"  refers  to  internally  generated  stratifi¬ 
cation  described  in  paragraph  20. 


INTERNALLY 

GENERATED 

INPUT  STRATIFICATION 

STRATIFICATION 


WHERE.  P-..,  -  IN  SITU  OVERBURDEN  STRESS 
FOR  STRATUM  N. 


Figure  4.  In  situ  stress  for  the  case  of  no  excavation 


LOG  jQ  STRESS,  P 

Figure  5.  In  situ  strain  for  strata  2,  3,  and 
4  from  piecewise  linear  stress-strain  relation¬ 
ship  and  the  case  of  no  excavation  (soil  profile 
is  shown  in  Figure  4) 


23.  If  the  stress-strain  relationship  is  defined  as  a  piecewise  linear 
curve  and  an  excavation  precedes  loading,  the  in  situ  properties  are  deter¬ 
mined  as  follows: 

a.  The  pre-excavation  overburden  stress  (P^)  is  found  by  summing 
the  contribution  of  all  material  from  the  pre-excavation  ground 
surface  to  the  midheight  of  the  stratum  of  interest  (see 
Figure  6) . 


INPUT 

STRATIFICATION 

EXCAVATION 


SOIL  ONE 


SOIL  TWO 


INTERNALLY 

GENERATED 

STRATIFICATION 


P  =  p  Ap 

POSTN  ON  “  N 


EQ.  8 


WHERE,  PnfJ  =  PRE-EXCAVATION  IN  SITU  OVERBURDEN 
STRESS  FOR  STRATUM  N 

Pp-.-M  =  POST  EXCAVATION  OVERBURDEN 
STRESS  FOR  STRATUM  N 

APn=  STRESS  RELIEF  FOR  STRATUM  N  DUE 
TO  EXCAVATION 


Figure  6.  Calculation  of  pre-excavation  and  post-excavation 

overburden  stresses 


b.  The  pre-excavation  strain  (e  )  is  found  by  interpolating 

between  appropriate  points  on  the  stress-strain  curve  (Figure  7,. 


LOG  10  STRESS,  P 


Figure  7.  Pre-excavation  and  post-excavation  strains 
from  piecewise  linear  stress-strain  relationship  for 
soil  profile  shown  in  Figure  6 

c.  The  stress  relief  (AP.,)  due  to  the  excavation  determined  bv 
the  Boussinesq  or  the  Westergaard  formula.  The  post-excavation 
stress  (PposTN^  t*ien  ^ounc^  by  subtracting  this  stress  relief 
from  the  pre-excavation  overburden  (see  Figure  6). 

d.  The  post-excavation  strain  ^pQgpfl)  is  found  by  computing  the 
change  in  strain  along  the  recompression  slope  (C  )  due  to  the 
stress  relief  and  adding  it  to  the  pre-excavation  strain  (e^) 

(Figure  7). 

24.  If  the  stress-strain  relationship  is  defined  as  a  bilinear  curve 
and  no  excavation  precedes  loading,  in  situ  properties  are  determined  as 

follows:  ! 

f 


i 
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b.  The  pre-excavation  strain  (c^)  is  found  by  interpolating 

between  appropriate  points  on  the  stress-strain  curve  (Figure  7) 


POST-EXCAVATION  CONDITIONS 
^  UNLOADING 


postn  r- —  —  T 
e  [  O. 
w  on  r  —  —  -*♦ 


■  PRE-EXCA  VA  TION  CONDITIONS 


LOADING 


Figure  7.  Pre-excavation  and  post-excavation  strains 
from  piecewise  linear  stress-strain  relationship  for 
soil  profile  shown  in  Figure  6 

c.  The  stress  relief  (AP^)  due  to  the  excavation  ir.  determined  by 
the  Boussinesq  or  the  Westergaard  formula.  The  post-excavation 
stress  (PposTfp  t^len  found  by  subtracting  this  stress  relief 
from  the  pre-excavation  overburden  (see  Figure  6). 

d.  The  post-excavation  strain  (GpQgpjg)  is  found  by  computing  the 

change  in  strain  along  the  recompression  slope  (C  )  due  to  the 

stress  relief  and  adding  it  to  the  pre-excavation  strain  (c„.,) 

ON 

(Figure  7). 

24.  If  the  stress-strain  relationship  is  defined  as  a  bilinear  curve 
and  no  excavation  precedes  loading,  in  situ  properties  are  determined  as 
follows: 
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The  in  situ  overburden  stress  (P^)  is  found  by  summing  the 
contributions  of  all  materials  from  the  ground  surface  to  the 
midheight  of  the  stratum  of  interest  (Figure  A). 

_b.  The  in  situ  strain  (e^)  is  found  by  calculating  the  strain 
corresponding  to  the  in  situ  stress  along  the  bilinear  curve. 

It  should  be  noted  that  if  the  in  situ  overburden  stress  is 
less  than  the  stress  defining  the  intersection  of  the  recompres¬ 
sion  slope  and  the  compression  slope,  the  in  situ  strain  will 
be  located  along  the  recompression  slope  (Figure  8). 


LOG jq  STRESS  P 

Figure  8.  Determination  of  in  situ  strain  for 
strata  2,  3,  and  A  from  bilinear  stress-strain 
relationship  and  the  case  of  no  excavation  (soil 
profile  is  shown  in  Figure  A) 

25.  If  the  stress-strain  relationship  is  defined  as  a  bilinear  curve  and 
an  excavation  precedes  loading,  the  following  procedure  is  used  to  determine 
the  in  situ  properties. 


a.  The  pre-excavation  in  situ  stress  (P..T)  is  found  by  summing  the 
contributions  of  all  materials  from  the  pre-excavation  ground 
surface  to  the  midheight  of  the  stratum  of  interest  (Figure  6). 

b.  The  pre-excavation  strain  (£qN)  is  found  by  calculating  the 
strain  corresponding  to  the  pre-excavation  stress  on  the  bilin¬ 
ear  curve  (Figure  9). 


LOG  10  STRESS.  P 


Figure  9.  Determination  of  pre-excavation  and 
post-excavation  strains  for  strata  3  and  4  from 
bilinear  stress-strain  relationship 

c.  The  stress  relief  due  to  the  excavation  is  determined  by  the 
Boussinesq  or  the  Westergaard  formula.  The  post-excavation 
stress  (ppQSTN)  is  found  by  subtracting  this  stress  relief  from 
the  pre-excavation  stress  (Figure  6). 

d.  The  post-excavation  strain  (epQgTN)  is  found  by  computing  the 
change  in  strain  along  the  recompression  slope  due  to  the 
stress  relief  and  adding  it  to  the  pre-excavation  strain  (see 
Figure  9) . 

Load  Definition 

26.  Four  options  are  available  for  defining  loads: 


a_.  Two-dimensional  soil  loads  which  are  load  density  profiles  de¬ 
fined  in  the  X-Z  plane  (Figure  10).  The  load  is  assumed  to 
extend  infinitely  in  the  Y  direction,  plane-strain  conditions. 

1).  Two-dimensional  pressure  loads  which  are  pressure  diagrams  de¬ 
fined  in  the  X-Z  plane  (Figure  11).  Again,  the  load  is  assumed 
to  extend  infinitely  in  the  Y-direction,  plane-strain 
conditions . 

c_.  Three-dimensional  rectangular  loads  which  are  rectangular  re¬ 
gions  of  uniform  pressure  defined  in  the  X-Y  plane  (Figure  12). 
ji.  Three-dimensional  irregular  loads  which  are  piecewise  linear 

regions  of  uniform  pressure  defined  in  the  X-Y  plane  (Figure  13). 


Figure  11.  Example  of  two-dimensional  pressure  loads 


E3  LOAO  AREA  ONE.  PEFF  «  100  PSF 

I - 1  LOAD  AREA  TWO.  PEFF  -  0  0  PSF 

I — I  LOAO  AREA  THREE.  PErF  ISO  PSF 

Figure  13.  Example  of  three-dimensional 
irregular  loads 
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27.  Two-dimensional  pressure  loads  and  two-dimensional  soil  loads  may 
be  combined  ("mixed")  to  define  a  loading  case.  When  mixed  loads  are  defined, 
they  should  be  numbered  consecutively  with  all  pressure  loads  preceding  the 
soil  loads.  The  three-dimensional  load  options  cannot  be  used  concurrently 
with  two-dimensional  loading. 

28.  The  superposition  techniques  employed  in  CSETT  require  that  certain 
guidelines  be  followed  when  defining  loads.  These  guidelines  will  be  ex¬ 
plained  in  the  following  paragraphs  and  are  illustrated  in  the  example  prob¬ 
lems  in  Part  V  and  in  Appendix  A.  It  should  be  noted  that  the  load  number 
associated  with  the  load  profile  is  not  indicative  of  the  time  of  load  applica¬ 
tion.  The  load  profile  number  indicates  the  order  required  by  the  program  to 
read  the  load  data.  The  time  of  load  application  is  supplied  as  a  separate 
variable  for  each  load. 

29.  Two-dimensional  soil  load  profiles  are  defined  by  X-Z  coordinates 
bounding  layers  of  constant  unit  weight.  Coordinate  pairs  defining  each  load 
profile  should  be  input  in  a  left  to  right  order.  When  ordering  sets  of  pro¬ 
files,  the  user  places  the  topmost  layer  first  and  follows  with  subsequent 
layers  in  order.  To  define  a  uniform  loading  condition  extending  to  infinity 
in  all  directions,  two  points  at  x-coordinates  of  -9999  and  9999  are  entered 
along  with  the  appropriate  load  elevation.  The  explanation  below  and  the 
examples  in  Appendix  A  should  aid  the  user  in  coding  difficult  loading  cases. 

30.  Consider  the  levee  and  berm  sections  shown  in  Figure  10. 

a.  Load  profile  1  is  defined  by  coordinates  A,  B,  C,  D,  E,  F,  G, 

H,  I,  and  J  and  is  assigned  a  unit  weight  of  Y,  where 

berm 

points  A  and  J  have  x-coordinates  of  -9999  and  +9999, 
respectively. 

_b.  Load  profile  2  is  defined  by  coordinates  A,  K,  E,  F,  L,  and  J 
and  is  assigned  a  unit  weight  of  T^eyee  • 

31.  The  induced  stresses  due  to  each  load  are  determined  as  follows: 

a_.  The  induced  stress  at  the  point  of  interest  is  first  computed 

assuming  that  the  entire  area  is  enclosed  by  profile  1  and  that 

the  ground  surface  contains  material  with  a  unit  weight  y, 

berm 

b.  The  actual  induced  stress  due  to  the  berms  is  determined  by 

computing  the  stress  induced  by  a  material  of  unit  weight 
Yberm  t*ie  area  enclosed  by  load  profile  2  and 
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subtracting  this  stress  from  that  determined  in  £  above. 

£.  The  stress  induced  by  the  levee  is  then  computed  using  a 

material  of  unit  weight  Y^gvee  within  the  area  enclosed  by 
load  profile  2. 

32.  Two-dimensional  pressure  loads  are  defined  by  coordinates  and  pres¬ 
sures  input  in  a  left  to  right  order  (Figure  11).  To  define  a  uniform  load 
extending  to  infinity  in  all  directions,  two  points  at  x-coordinates  of  -9999 
and  +9999  with  the  appropriate  pressure  should  be  supplied.  Multiple  pressure 
profiles  may  be  entered.  Each  is  considered  individually  to  determine  induced 
stresses.  The  total  induced  stress  for  a  given  location  is  the  sum  of 
stresses  computed  for  individual  loads. 

33.  Three-dimensional  rectangular  loads  are  defined  by  supplying  the 
X-Y  coordinates  of  the  lower  left  corner  and  upper  right  corner  (plan  view)  of 
each  rectangular  region  of  uniform  pressure  (Figure  12).  The  stress  induced 
by  each  region  is  calculated  separately.  The  total  induced  stress  for  a  given 
location  is  the  sum  of  the  stresses  contributed  bv  individual  loads. 

34.  The  three-dimensional  irregular  load  option  allows  the  user  to  de¬ 
fine  multiple  areas  by  piecewise  linear  perimeters  in  the  X-Y  plane  (plan 
view)  with  uniform  pressure  loads.  The  stress  at  each  location  is  determined 
by  summing  the  individual  stresses  computed  for  each  load  area.  Loads  must  be 
defined  in  a  clockwise  order,  and  the  pressure  assigned  to  each  area  must 
represent  the  resultant  pressure  after  considering  any  overlying  loads.  A 
superposition  code  has  been  included  as  an  optional  means  of  controlling  pres¬ 
sures  assigned  to  each  area  during  the  stress  computations.  When  a  super¬ 
position  code  is  assigned  to  a  load  area  the  pressure  assigned  to  that  area 
during  stress  induction  computations  is 


P 


eff 


(NL) 


P.  (NL)  -  P. 

input  input 


(NS) 


(9) 


where 

NL  =  load  number  of  the  load  being  considered 

NS  =  superposition  code  which  is  the  load  number  of 
the  load  immediately  surrounding  the  load  being 
considered 


P, 

input 

P, 

input 


(NL) 

(NS) 


=  pressure  assigned  to  load  NL  by  user 
=  pressure  assigned  to  load  NS  by  user 
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P  ,  (NL)  =  pressure  used  to  compute  stresses  induced  by 
ett  load  NL 

The  following  example  illustrates  the  uses  of  the  three-dimensional  irregular 
load  option  with  and  without  using  the  superposition  code  (Figure  13). 

35.  If  the  loads  are  to  be  input  using  superposition  codes: 

a.  Load  one  is  assigned  a  pressure  of  100  psf  and  a  superposition 
code  of  zero  (0) . 

b.  Load  two  is  assigned  a  pressure  of  0.0  psf  and  a  superposition 
code  of  one  ( 1) . 

c.  Load  three  is  assigned  a  pressure  of  150  psf  and  a  superposition 
code  of  two  (2) . 

36.  The  induced  stress  at  the  point  being  considered  is  computed  as 
described  below. 

a.  First  a  stress  is  computed  assuming  a  100-psf  pressure  over 
the  entire  area  bounded  by  the  perimeter  of  load  one.  Note 
that  this  area  includes  areas  occupied  by  loads  two  and  three 
whose  associated  pressures  are  not  100  psf. 

b.  Next,  stresses  induced  by  a  pressure  of 

Peff  (2)  =  PinPut  (2)  "  PinPut  (1)  *  °  '  100  =  ~100  ”sf 
acting  over  the  area  bounded  by  the  perimeter  of  load  two  are 

added  to  the  stresses  computed  in  Step  a.  This  effectively 

superimposes  a  zero  (0)  load  region  within  the  100-psf  load 

region  one. 

c.  Lastly,  stresses  induced  by  a  pressure  of 

P  ,,  (3)  =  P.  _  (3)  -  P.  (2)  =  150  -  0  =  150  psf 

eff  input  input 

acting  over  the  area  bounded  by  the  perimeter  of  load  three  are 
added  to  the  stresses  computed  in  Step  b.  The  resulting  stress 
is  the  total  induced  stress  resulting  from  the  three  loads. 

37.  If  the  loads  are  to  be  input  without  using  superposition  codes: 

a.  Load  one  is  assigned  a  pressure  of  100  psf  and  a  superposition 
code  of  zero  (0) . 

b.  Load  two  is  assigned  a  pressure  of  -100  psf  and  a  superposition 
code  of  zero  (0). 

c.  Load  three  is  assigned  a  pressure  of  150  psf  and  a  superposi¬ 
tion  code  of  zero  (0). 

38.  The  induced  stress  at  a  point  being  considered  is  computed  by 


summing  the  values  computed  using  the  assigned  pressure  associated  with  each 
area.  Note  that  when  superposition  codes  are  not  used,  the  user  is  manually 
accomplishing  the  same  results  generated  internally  with  the  superposition 


PART  IV:  INPUT  GUIDE 


Source  of  Input 

39.  Input  data  may  be  supplied  interactively  from  the  terminal  or  from  a 
predefined  data  file.  Data  supplied  interactively  from  the  terminal  may  be 
entered  by  following  a  prompting  sequence  or  by  using  the  same  keyword  format 
as  that  used  to  create  a  predefined  data  file. 


Data  Format 


40.  All  input  data  (interactive  or  data  file)  are  read  in  free-field 
format . 

a.  Data  items  must  be  separated  by  one  or  more  blanks  (comma  separa¬ 
tors  are  not  permitted) . 

b.  Integer  numbers  must  be  in  the  form  NNNN. 

c.  Real  numbers  may  be  in  the  form  ±  xxxx,  ±xx.xx,  or  ±  xx.xxE+rr. 

Data  Entry  from  Terminal 

41.  Data  are  entered  from  the  terminal  using  the  prompt  or  no-prompt  mode. 
In  the  prompt  mode,  the  user  is  queried  for  the  required  data.  In  the  no-prompt 
mode,  data  are  input  by  using  keywords  followed  by  the  associated  data.  A  de¬ 
scription  of  the  keyword  format  is  described  in  the  section  titled  Input 
Description.  Note  that  the  format  for  the  no-prompt  mode  is  the  same  as  that 
for  a  data  file  except  that  line  numbers  are  not  used  when  entering  data  from 
the  terminal. 


Data  Entry  from  File 

42.  Using  a  data  file  allows  the  user  to  create  a  permanent  file  of  the 
input  data  prior  to  executing  the  program.  To  create  a  data  file,  follow  the 
guidelines  provided  in  the  section  Input  Description. 

43.  In  addition  to  the  general  format  requirements  given  in  paragraph  40, 
the  following  pertain  to  data  being  entered  using  the  Input  Description  guide. 
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a.  Each  line  of  a  data  file  must  begin  with  a  nonzero,  positive  line 
number,  denoted  by  LN.  Line  numbers  are  not  used  when  entering 
data  from  the  terminal. 

b.  A  line  of  input  may  require  a  keyword.  The  acceptable  abbrevia¬ 
tion  for  the  keyword  is  indicated  by  underlined  capital  letters, 
e.g.,  the  acceptable  abbreviation  for  the  keyword  BOUSsinesq  is 
BOUS. 

c.  Alphanumeric  data  items  are  enclosed  in  quotes. 

d.  Items  designated  by  upper  case  letters  and  numbers  without  quotes 
indicate  numeric  data  values.  Numeric  data  values  are  either 
real  or  integer  according  to  standard  FORTRAN  variable  naming 
conventions. 

e.  Data  items  enclosed  in  brackets  []  may  be  omitted  if  default 
values  are  to  be  used.  Data  items  in  parentheses  ()  indicate 
that  a  special  note  follows. 

_f.  Input  data  are  divided  into  the  sections  listed  in  paragraph  44 
below.  Each  section  consists  of  a  keyword  line  and  may  be 
followed  by  one  or  more  data  lines. 

g.  Each  input  line  is  limited  to  70  characters  including  LN. 


Sections  of  Input 

44.  Input  data  are  divided  into  the  following  sections: 

a.  Title. 

b.  Applied  loads. 

£.  Excavation. 

d.  Soil  stratification. 

e.  Stress-strain  properties. 

£.  Type  and  extent  of  stress  analysis. 

£.  Time  sequence  for  consolidation  calculations. 

h.  Location  of  planes  of  interest. 

i.  Termination. 


Units 


45.  All  data  must  be  In  feet,  pounds,  and  years. 
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Input  Description 


46.  Title 


a.  Keyword  line — one  (1)  per  analysis 
[LN]  "TITLe" 

b.  Data  lines — one  (1)  per  analysis 
b(l)  Contents 

[LN]  title 
b(2)  Definition 

title  =  alphanumeric  information  up  to  70  characters  including 


47.  Applied  Loads 

a.  2-Dimensional  Pressure  Load 

a(l)  Keyword  line — one  (1)  per  load  (max  of  20  loads) 

(a)  Contents 

[LN]  "2DPRessure"  NL  N  TIM1  TIM2 

(b)  Definitions 
''2DPRessure"  =  keyword 
NL  =  load  number 

N  »  number  of  points  defining  load  NL  (max  of  50) 

T1M1  =  beginning  time  of  load  application  (yrs) 

TIM2  =  ending  time  of  load  application  (yrs) 
a(2)  Data  lines — one  (1)  to  fifty  (50)  points  per  load 

(a)  Contents 

[LN]  XI  PI  X2  P2  ... 

(b)  Definitions 

XI  =  x-coordinate  of  point  1  defining  pressure  load  (ft) 
PI  =  pressure  corresponding  to  location  1  defining 
pressure  load  (psf)  etc., 

b.  2-Dimensional  Soil  Load  Profile 

b(l)  Keyword  line — one  (1)  per  load  (max  of  20  loads) 

(a)  Contents 

[LN]  "2DS0il"  (NL)  N  T1M1  T1M2  GAM 

(b)  Definitions 
"2DS011"  -  keyword 

(NL)  -  load  profile  number 


N  =  number  of  points  defining  load  NL  (max  of  50) 

T1M1  =  beginning  time  of  load  application  (yrs) 

T1M2  =  ending  time  of  load  application  (yrs) 

GAM  =  effective  unit  weight  of  soil  load  (pcf) 

(c)  Discussion 

As  noted  previously,  the  load  profile  number,  NL,  indi¬ 
cates  a  superposition  sequence  and  should  not  be  mis¬ 
taken  for  the  chronological  order  of  load  application. 

£(2)  Data  line — One  (1)  to  fifty  (50)  lines  per  load 

(a)  Contents 

[LN]  XI  Z1  X2  Z2  ... 

(b)  Definitions 

XI  =  x-coordinate  of  point  1  defining  soil  load  NL  (ft) 
Z1  =  z-coordinate  of  point  1  defining  soil  load 
NL  (ft),  etc. 

£.  3-Dimensional  Rectangular  Load 

£(1)  Keyword  line-one  (1)  per  load  (max  of  25  loads) 

(a)  Contents 

[LN]  "3DREct angular"  NL  T1M1  T1M2  PRES 

(b)  Defintions 
"3DREctangular"  =  keyword 
NL  =  load  number 

T1M1  =  beginning  time  of  load  application  (yrs) 

T1M2  =  ending  time  of  load  application  (yrs) 

PRES  =  pressure  of  load  region  (psf) 

£(2)  Data  lines — one  (1)  line  per  load 

(a)  Contents 

[LN]  XL  YL  XR  YR 

(b)  Definitions 

XL  =  x-coordinate  (ft)  defining  lower  left  corner  of 
load  region  (PLAN  VIEW) 

YL  *  y-coordinate  (ft)  defining  lower  left  corner  of 
load  region  (PLAN  VIEW) 

XR  =  x-coordinate  (ft)  defining  upper-right  corner  of 
load  region  (PLAN  VIEW) 
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YR  =  y-coordinate  (ft)  defining  upper  right  corner  of 
load  region  (PLAN  VIEW) 
d.  3-Dimensional  Irregular  Load 

d(l)  Keyword  line — one  (1)  per  load  (max  of  20  loads) 

(a)  Contents 

[LN]  "3DIRregular"  NL  (NS)  N  T1M1  T1M2  PRES 

(b)  Definitions 
"3DIRregular"  =  keyword 
NL  =  load  number 

(NS)  =  superposition  code 

N  =  number  of  points  defining  load  region  (max  of  50) 
T1M1  =  beginning  time  of  load  application  (yrs) 

T1M2  =  ending  time  of  load  application  (yrs) 

PRES  =  pressure  of  load  region  (psf) 

(c)  Discussion 

The  superposition  code  is  the  load  number  of  the 
surrounding  load.  For  external  loads,  the  super¬ 
position  code  is  zero  (0).  For  interior  loads,  the 
superposition  code  is  the  number  corresponding  to  the 
load  region  immediately  outside  (surrounding)  the  load. 
d(2)  Data  lines — one  (1)  to  fifty  (50)  points  per  load 

(a)  Contents 

[LN]  XI  Y1  X2  Y2 . . . 

(b)  Definitions 

XI  *  x-coordinate  of  point  1  defining  load  region  (ft) 
Y1  *  y-coordinate  of  point  1  defining  load  region  (ft) 
etc. 

(c)  Discussion 

All  loads  must  be  entered  in  a  clockwise  order. 

48.  Excavation 

a.  Keyword  line — one  (1)  per  analysis 
a(l)  Contents 

[LN]  "EXCAvation"  N  GAM 


a(2)  Definitions 

"EXCAvation”  -  keyword 


N  =  number  of  points  defining  excavated  profile  (max  of  50) 
GAM  =  effective  unit  weight  of  excavated  material  (pcf) 

Data  line — one  (1)  to  fifty  (50)  lines 
b^(l)  Contents 

[LN]  XI  ZI  X2  Z2. . . 

]j(2)  Definitions 

XI  =  x-coordinate  of  point  1  defining  excavation  profile  (ft) 
Zl  =  z-coordinate  of  point  1  defining  excavation  profile 
(ft)  etc. 

Comments 

For  excavations  assumed  to  extend  infinitely  in  all 
directions,  the  excavation  profile  must  be  defined  with  two 
points  at  the  pre-excavation  ground  surface  elevation  at 
infinite  x-coordinates .  This  format  is  also  used  to  define 
excavations  for  3-D  Load  Cases  (XI  =  -9999  X2  =  9999). 

49.  Soil  Stratification 

£.  Keyword  line — one  (1)  per  stratum 
£(1)  Contents 

[LN]  "SOIL"  N  EL  ("NDC")  GM  [CR]  [CV]  [ST] 
a^2)  Definition 

"SOIL"  =  keyword 

N  =  stratum  number 

EL  =  elevation  of  top  stratum  (ft) 

(NDC)  *=  drainage  condition 

GM  =  effective  unit  weight  of  stratum  (pcf) 

[CR]  =  recompression  index  (not  required  if  NDC  =  N) 

2 

[CV]  =  coefficient  of  consolidation  (ft  /yr)  (not  required 
if  NDC  =  N) 

[ST]  =  Poisson's  ratio  (default  =  0.32) 
jj(3)  Discussion 

(a)  Each  soil  stratification  line  must  be  immediately 
followed  by  the  stress-strain  information  for  that 
stratum  (see  paragraph  50)  unless  an  incompressible 
stratum  is  input.  No  stress-strain  information  is 
required  for  incompressible  strata. 
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50. 


(b)  The  drainage  condition  indicates  the  drainage  assump¬ 
tion  for  time-rate  of  consolidation  calculations.  A 
value  of  S  indicates  single  drainage,  D  indicates 
double  drainage,  and  C  indicates  the  stratum  should  be 
combined  with  the  previous  stratum.  A  value  of  N  indi¬ 
cates  an  incompressible  stratum. 

_b.  Data  lines  -  none 
Stress-Strain  Properties 
a_.  Void  Ratio-Pressure  Curve 

a^(l)  Keyword  line — one  (1)  per  stratum 
(a)  Contents 


[LN]  "VOID"  NP 
(b)  Definitions 

"VOID"  =  keyword 

NP  =  number  of  points  defining  curve  (max  of  10  pts) 
a_(2 )  Data  lines — one  (1)  to  five  (5)  lines 
(a)  Contents 


[LN]  0R1  ABI  0R2  AB2. . . 

(b)  Definitions 

OR I  =  void  ratio  value  of  first  point  on  curve 

ABI  =  pressure  value  of  first  point  on  curve  (psf), 

etc. 

1>.  Strain-Pressure  Curve 

_b(l)  Keyword  line — one  (I)  per  stratum 

(a)  Contents 

[LN]  "STRAin"  NP 

(b)  Definitions 
"STRAin"  -  keyword 

NP  =  number  of  points  defining  curve  (max  of  10  pts) 
_b(2)  Data  line — one  (1)  to  five  (5)  lines 
(a)  Contents 


[LN]  OR I  ABI  0R2  AB2 . . . . 

(b)  Definitions 

0R1  =  strain  value  of  first  point  on  curve 

ABI  =  pressure  value  of  first  point  on  curve  (psf) 
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£.  Index  Properties 

£(1)  Keyword  line — one  (1)  line  per  stratum 

(a)  Contents 

[LN]  "INDEx"  CC  PO  EO 

(b)  Definitions 
"INDEx"  =  keyword 

CC  =  compression  index 
PO  =  in  situ  pressure  (psf) 

EO  =  in  situ  void  ratio 
£(2)  Data  lines  -  none 

51.  Type  of  Stress  Analysis  -  one  (1)  line  per  analysis 
a_.  Boussinesq  Analysis 

£(1)  Keyword  line — one  (1)  per  analysis 

(a)  Contents 

[LN]  "BOUSsinesq"  DMAX 

(b)  Definitions 
"BOUSsinesq"  =  keyword 

DMAX  =  maximum  depth  below  top  of  strata  1  to  which 
analysis  will  be  extended  (ft) 
a(2)  Data  lines  -  none 
j).  Westergaard  Analysis 

_b(l)  Keyword  line — one  (I)  per  analysis 

(a)  Contents 

[LN]  "WESTergaard"  DMAX 

(b)  Definitions 
"WESTergaard"  =  keyword 

DMAX  =  Maximum  depth  below  top  of  stratum  1  to  which  the 
analysis  will  extend 
jj(2)  Data  lines  -  none 

52.  Time  Sequence  for  Consolidation  Calculations 
a.  Geometric  Progression 

£(1)  Keyword  line — one  (1)  per  analysis 
(a)  Contents 

[LN]  "TIMI"  [VTS ] [ALF ] 


\ 
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(b)  Definitions 

"TIHI"  =  keyword 

[VTS]  =  initial  time  period  (default  =  1  week  = 

0.0192  yrs) 

[ALF]  =  multiplicative  factor  for  the  geometric  pro¬ 
gression  algorithm  (default  =  2.0) 

(2)  Data  lines  -  none 

(3)  Discussion  -  a  total  of  thirteen  (13)  computation  times  are 
generated 

b^.  Specified  Times 

_b(l)  Keyword  line — one  (1)  line  per  analysis 

(a)  Contents 

[LN]  "TIMS”  T 1  T2 . . . 

(b)  Definitions 
"TIMS"  =  keyword 

T1  =  first  time  for  which  time-rate  of  consolidation 
will  be  determined  (yrs). 

T2  =  second  time  for  which  time-rate  of  consolidation 
will  be  determined  (yrs),  etc. 
d(2)  Data  lines  -  none 
53.  Location  of  Planes  of  Interest 
a_  2 -Dimensional  Cases 

a/1)  Keyword  line — one  (1)  per  analysis 

(a)  Contents 

[LN]  "OUTPut"  XL  XR  SP 

(b)  Definition!, 

"OUTPut"  -  keyword 

XL  =  leftmost  point  defining  plane  of  interest  (ft) 

XR  =  rightmost  point  defining  plane  of  interest  (ft) 

SP  =  spacing  along  plane  of  interest  where  output  is 
desired  (ft) 
a(2)  Data  lines  -  none 
J).  3-Dimensional  Case 

b^(l)  Keyword  line — one  (1)  to  the  number  of  planes  at  which 
output  is  desired 
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(a)  Contents 

[LN]  "OUTPut"  XI  Y1  X2  Y2  SP 

(b)  Definitions 
"OUTPut"  =  keyword 

XI  =  x-coordinate  (ft)  defining  lower  (or  left)  corner 


of  plane  of  interest  (plan  view) 

Y1  =  y-coordinate  (ft)  defining  lower  (or  left) 
corner  of  plane  of  interest  (plan  view) 

X2  =  x-coordinate  (ft)  defining  upper  (or  right)  cor¬ 
ner  of  plane  of  interest  (plan  view) 

Y2  «  y-coordinate  (ft)  defining  upper  (or  right)  cor¬ 
ner  of  plane  of  interest  (plan  view) 

SP  =  spacing  along  plane  of  interest  where  output  is 
desired  (ft) 

_b(2)  Data  lines  -  none 

54.  Termination  of  Input  Sequence 

a.  Keyword  line-one  (1)  line  required 
£(1)  Contents 

[LN]  "END" 

Abbreviated  Input  Guide 

55.  Title — one  (I)  line  identifying  problem 
a..  Keyword  line  —  one  (1)  per  analysis 

[LN]  "TITLe"  title 
_b.  Data  line — one  (1)  per  analysis 
[LN]  title 

56.  Applied  Loads 

£.  2-Dimensional  Pressure  Load 

£(1)  Keyword  line — one  (1)  per  load  (max  of  20  loads) 

[LN]  "2DPRessure"  NL  N  T1M2  T1M2 
a(2)  Data  lines — one  (1)  to  fifty  (50)  lines  per  load. 

[LN]  XI  PI  X2  P2. . . . 
b^.  2-Dimensional  Soil  Load  Profile 

]>(1)  Keyword  line — one  (1)  per  load  (max  of  20  loads) 

[LN]  "2DS0il"  (NL)  N  T1M1  T1M2  GAM 


l)(2)  Data  line — one  (1)  to  fifty  (50)  lines  per  load 
[LN]  XI  Z1  X2  Z2 . . . . 

£.  3-Diraensional  Rectangular  Load 

£(1)  Keyword  line — one  (1)  line  per  load  (max  of  25  loads) 
[LN]  "_3DREct  angular"  NL  T1M1  T1M2  PRES 
£(2)  Data  lines — one  (1)  line  per  load 
[LN]  XI  YI  XL  YL 

jl.  3-Dimensional  Irregular  Load 

<i(l)  Keyword  line — one  (1)  line  per  load  (max  of  20  loads) 
[LN]  "3DIRregular"  NL  (NS)  N  T1M1  T1M2  PRES 
cl(2)  Data  lines — one  (1)  to  fifty  (50)  lines  per  load 
[LN]  XI  Yl  X2  Y2 . . . 

57.  Excavation  (2-Dimensional  Analyses  Only) 

£.  Keyword — one  (1)  line  per  analysis 

[LN]  "EXCAvation"  N  GAM 
_b.  Data  line — one  (1)  to  fifty  (50)  lines 
[LN]  XI  Z1  X2  Z2... 

58.  Soil  Stratification 

a.  Keyword — one  (1)  line  per  stratum 

[LN]  "SOIL"  N  EL  ("NDC")  GM  [CR]  [CV]  [ST] 
ib.  Data  lines  -  none 

59.  Stress-Strain  Properties 

a_.  Void  Ratio — Pressure  Curve 

£(1)  Keyword — one  (1)  line  per  stratum 
[LN]  "VOID”  NP 

a(2)  Data  lines — one  (1)  to  five  (5)  lines 
[LN]  0R1  AB1  0R2  AB2... 

Strain-Pressure  Curve 
_b ( 1 )  Keyword — one  (1)  line  per  stratum 
[LN]  "STRAin"  NP 

Jb(2)  Data  line — one  (1)  to  five  (5)  lines 
[LN]  OR I  AB1  0R2  AB2... 

£.  Index  Properties 

£(1)  Keyword — one  (1)  line  per  stratum 
[LN]  "INDEx"  CC  PO  E0 
c(2)  Data  lines  -  none 


60.  Type  of  Stress  Analysis — one  (1)  line  per  analysis 

Boussinesq  Analysis 

Ml)  Keyword — one  (1)  line  per  analysis 
[LN]  "BOUSsinesq"  DMAX 
M2)  Data  lines  -  none 
J5.  Westergaard  Analysis 

Ml)  Keyword — one  (1)  line  per  analysis 
[LN]  "WESTergaard"  DMAX 
M2)  Data  lines  -  none 

61.  Time  Sequence  For  Consolidation  Calculations 
a_.  Geometric  Progression 

a(l)  Keyword — one  (1)  line  per  analysis 
[LN]  "TIMI"  [VTS ]  [ALF] 

M2)  Data  lines  -  none 
M  Specified  Times 

Ml)  Keyword — one  (1)  line  per  analysis 
[LN]  "TIMS"  T1  T2  T3... 

M2)  Data  lines  -  none 

62.  Location  of  Plane  of  Interest 
a.  2-Dimensional  Case 

Ml)  Keyword-one  (1)  line  per  analysis 
[LN]  "OUTPut"  XL  XR  SP 
M 2)  Data  lines  -  none 
_b.  3-Dimensional  Case 

Ml)  Keyword — one  (1)  to  the  number  of  planes  at  which  output  is 
desired 

[LN]  "OUTPut"  XI  Y1  X2  Y2  SP 
M2)  Data  lines  -  none 

63.  Termination  of  Input  Sequence 
[LN]  "END" 


PART  V:  EXAMPLE  PROBLEMS 

64.  The  following  example  problems  are  solved  using  the  program  CSETT. 
These  examples  illustrate  the  input  requirements  for  the  program. 


Example  Problem  1 

65.  In  this  problem  the  ultimate  settlement  and  time-settlement  rates  re¬ 
sulting  from  a  9-ft-high  uniform  fill  are  determined.  The  soil  profile  and  soil 
properties  are  shown  in  Figure  14.  The  input  data  are  presented  in  Figure  15. 
The  clay  stratum  is  assumed  to  be  drained  only  at  its  upper  boundary.  The 
soil  is  normally  consolidated.  Results  of  this  analysis  can  be  compared  to 
those  presented  by  Hough  (1957). 

£L  9 

- - - - 1 - 1 - - - - - 

V,tt  -  ns  Kf 

...  1 1 1 1  n  1 1 1 1 1 

SAND  ^SATm,2*pC* 


7sub  ■ 50  PCF 
Cc  -  0.40 

•0-».4 

C*-  120  FT*/YR 


Figure  14.  Soil  profile  and  properties  for  Example  Problem  1 


010  TITLE 

020  example  problem  prom  HOUOH 

0  >1S  (DIMENSIONAL  SOIL  LOAD  DATA  ML  N  TIM1  TIM2  GAM) 

040  -9999  9  9999  9  (X.  Z  _  i  -  1  N) 

050  SOIL  1  0  N  125  (SOIL  DATA  N  EL  nAc  CM) 

060  SOIL  2  -14  S  50  .4  120  (SOIL  DATA  N  EL  NDC  CM  CR  CV) 

070  INDEX  .4  2000  1.4  (INDEX  PROPERTIES  CC  PO  EO) 

0S0  SOIL  3  -24  N  52.5  (SOIL  DATA  N  EL  NDC  CM) 

090  BOIJS  25  (BOUSSINESQ  SOLUTION  DMAX) 

100  TIMI  (GEOMETRIC  PROGRESSION) 

110  OUTPUT  0  3  10  (OUTPUT  XL  XR  SP) 

120  END  (TERMINATION) 


Figure  15.  Input  data  file  for  Example  Problem  1 
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PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/83  TIME:  16:  3:22 


I.  INPUT  DATA 


1.  TITLE  -  EXAMPLE  PROBLEM  FROM  HOUGH 


2.  BOUSSINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  25.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 


PROFILE  NUMBER 


: NUMBER  OF  POINTS®  2 
BEGINNING  TIME  OF  APPLICATION  =  0. 
ENDING  TIME  OF  APPLICATION  =  0. 
EFFECTIVE  UNIT  WEIGHT  OF  SOIL  LOAD® 


YRS. 

YRS. 

115.00  PCF 


POINT  NO.  X  Y 

(FT.)  (FT.) 

1  -9999.00  9.00 

2  9999.00  9.00 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 
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7.  EXCAVATION  DATA 
NONE 


3.  SOIL  DATA 


STRATA  EL.  OF  TOP  DRAINAGE  EFF  UNIT  RECOMPR.  COEF.OF  POISSON'S 
NO.  OF  STRATUM  CONDITION  WEIGHT  INDEX  CONSOL.  RATIO 

(FEET  NGVD )  (PCF)  (SQFT/YR) 

1  0.  N  125.00 

2  -14.00  S  50.00  0.40000  120.00000  0.32000 


9.  STRESS-STRAIN  DATA 


STRATUM  NO.  1 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  2 


COMPRESSION  INDEX*  0.40000 
RECOMPRESS I ON  INDEX*  0.40000 
I NS ITU  VOID  RAT 10=  1.40000 
I NS ITU  OVERBURDEN*  2000.00  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


A  GEOMETRIC  PROGRESSION  WITH  AN  INITIAL  TIME  PERIOD  OF  0.019 
YEARS  AND  A  MULTIPLICATIVE  FACTOR  OF  2.0000  WILL  BE  USED  IN 
THE  TIME  RATE  OF  CONSOLIDATION  CALCULATIONS. 


11.  OUTPUT  CONTROL  DATA 


XXL=  0.  FT. 

XUL=  5.0000  FT. 

DELX=  10.0000  FT. 


WILL  OUTPUT  GO  TO  TERMINAL,  FILE,  OR  BOTH^> 
ENTER  T,  F,  OR  B 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/83  TIME:  16:  4:24 


II.  OUTPUT  SUMMARY. 


1.  TITLE-  EXAMPLE  PROBLEM  FROM  HOUGH 


POSITION:  X=  0. 

■a****************** 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 

STRATA  MID-DEPTH  IN-SITU  DELTA  ULTIMATE 
NO.  OF  STRATA  OVERBURDEN  SIGMA  SETTLEMENT 
(FEET)  (LB/SQ  FT)  (LB/SQ  FT)  (FEET) 


1  7.00  875.00  1035.00  0. 

2  19.00  2000.00  1035.00  0.303 


3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

ULT 

0.02 

0.04 

0.  08 

0.  15 

0.31 

0.61 

NO 

(YRS. ) 

(YRS. ) 

(YRS. > 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.303 

0.052 

0.074 

0.  103 

0.  146 

0.205 

0.262 

TOTALS: 

0.303 

0.052 

0.074 

0.  103 

0.  146 

0.205 

0.262 

(SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

STRATA 

1.23 

2.45 

4.91 

9.  82 

19.64 

39.  28 

78.  55 

NO 

( YRS. ) 

(YRS. ) 

(YRS. > 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.296 

0.303 

0.  303 

0.  303 

0.  303 

0 . 303 

0.  303 

TOTALS: 

0.296 

0.303 

0.303 

0.  303 

0.  303 

0.  303 

0.303 

Solution  by  Hough  (1957) 
AH  *  3.6  in. 
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Example  Problem  2 

66.  In  this  example  the  ultimate  settlement  and  time-settlement  rates  re 
suiting  from  a  uniform  load  of  1500  psf  are  determined.  The  soil  profile  con¬ 
sists  of  two  homogeneous  clays  above  an  impervious  layer  separated  by  a  free- 
draining  layer  of  incompressible  sand.  The  soil  profile  and  soil  properties 
are  shown  in  Figure  16.  The  input  data  file  is  shown  in  Figure  17.  Results  o 
this  analysis  can  be  compared  to  those  presented  by  Means  and  Parcher  (1963). 


010  TITLE 

020  EXAMPLE  PROBLEM  FROM  MEANS  !<  PARCHER 

030  2DPR  1200  (2-DIMENSIONAL  PRESSURE  LOAD  DATA  NL  N  TIMI  TIM2) 

040  -9999  1500  9999  1500  (X±  Z.  -  I  =  1^) 

050  SOIL  1  0  D  116.2  .36  39.42  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 
060  INDEX  .36  2000  1.02  (INDEX  PROPERTIES  CC  PO  EO) 

070  SOIL  2  -4  C-  53.3  .36  39. 42  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

080  INDEX  .36  2000  1.02  (INDEX  PROPERTIES  CC  PO  EO) 

090  SOIL  3  -10  N  60  (SOIL  DATA  N  EL  NDC  GM) 

100  SOIL  4  -12  S  61.6  .24  23.91  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 
110  INDEX  .24  3200  .73  (INDEX  PROPERTIES  CC  PO  EO) 

120  SOIL  5  -25  N  60  (SOIL  DATA  N  EL  NDC  GM) 

130  BOUS  30  (BOUSSINESQ  SOLUTION  DMAX) 

140  TIMI  (GEOMETRIC  PROGRESSION) 

150  OUTPUT  0  10  20  (OUTPUT  XL  XR  SP) 

160  END  (TERMINATION) 

* 


Figure  17.  Input  data  for  Example  Problem  2 


1 PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/83  TIME:  15:44:44 


I.  INPUT  DATA 


1.  TITLE  -  EXAMPLE  PROBLEM  FROM  MEANS  PARCHER 


2.  BOUSSINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  30.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 


LOAD  NUMBER  1  :  NUMBER  OF  POINTS=  2 

BEGINNING  TIME  OF  APPLICATION  =  0.  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 


POINT  NO.  X 

(FT.  ) 


PRESSURE 
( PSF ) 


1  -9999.00 

2  9999 . 00 


1500.00 

1500.00 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


!  7.  EXCAVATION  DATA 

[  .  NONE 

i  c 
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.&  «  W 


8.  SOIL  DATA 


STRATA  EL.  OF  TOP  DRAINAGE  EFF  UNIT  RECOMPR.  COEF.OF 


• 

o 

z 

OF  STRATUM 
(FEET  NGVD) 

CONDITION 

WEIGHT 
( PCF ) 

INDEX 

CONSOL. 

(SQFT/YR) 

1 

0. 

D 

116. 20 

0 . 36000 

39. 42000 

o 

-J 

-4 . 00 
-10. 00 

C 

N 

53.80 

60.00 

0.36000 

39. 42000 

4 

-12. 00 

»Z" 

61.60 

0. 24000 

23.91000 

9.  STRESS-STRAIN  DATA 

STRATUM  NO.  1 


COMPRESSION  I NDEX  =  0.36000 
RECOMPRESSION  INDEX*  0.36000 
I NS ITU  VOID  RAT 10=  1.02000 
I NS ITU  OVERBURDEN-  2000.00  PSF 


STRATUM  NO.  2 


COMPRESSION  INDEX=  0.36000 
DECOMPRESSION  INOEX=  0.36000 
I NS ITU  VOID  RAT 10=  1.02000 
I NS ITU  OVERBURDEN*  2000.00  PSF 


STRATUM  NO. 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  4 


COMPRESSION  INDEX=  0.24000 
RECOMPRESSION  INDEX*  0.24000 
I NS ITU  VOID  RATIO*  0.73000 
INSITU  OVERBURDEN*  3200.00 


PSF 
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POISSON 

RATIO 


0.32000 
0 . 32000 

0.32000 


■kr< 


.  4-. 


T  *  »  9t 


I 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


A  GEOMETRIC  PROGRESSION  WITH  AN  INITIAL  TIME  PERIOD  OF  0.0192 
YEARS  AND  A  MULTIPLICATIVE  FACTOR  OF  2.0000  WILL  BE  USED  IN 
THE  TIME  RATE  OF  CONSOLIDATION  CALCULATIONS. 


11.  OUTPUT  CONTROL  DATA 


XXL=  0.  FT. 

XUL=  10.0000  FT. 

DELX=  20.0000  FT. 


PROGRAM  CSETT  -  VERTICAL 
DATE:  7/  5/33 


II.  OUTPUT  SUMMARY. 


STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 

TIME:  15:41:  4 


1.  TITLE-  EXAMPLE  PROBLEM  FROM  MEANS  ?<  PARCHER 


POSITION:  X=  0. 

******************* 


SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 


MID- DEPTH  IN-SITU 

OF  STRATA  OVERBURDEN 


DELTA 

SIGMA 


ULTIMATE 

SETTLEMENT 


(FEET) 

(LB/SQ  FT) 

(LB/SQ  FT) 

(FEET) 

1 

2.00 

232. 40 

1 500. 00 

0.561 

2 

7.00 

626. 20 

1 500 . 00 

0.524 

3 

1 1 . 00 

347.60 

1500. 00 

0. 

4 

13.50 

1 303 . 00 

1500. 00 

0.575 

I,  *1  ;f  «*. 


3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

ULT 

0.02 

0.04 

o.os 

0.  15 

0.  31 

0.61 

NO 

( YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

0.561 

0.111 

0.  156 

0.220 

0.  31 1 

0.424 

0.520 

O 

0.524 

0.  103 

0.  145 

0.  205 

0.289 

0.395 

0.484 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.575 

0.  034 

0.047 

0.068 

0.  096 

0.  136 

0.  192 

TOTALS: 

1 . 660 

0.248 

0.348 

0.  493 

0.696  < 

H.955 

1.  196 

( SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

STRATA 

1.23 

2.45 

4.91 

9.82 

19.64 

39.28 

78.55 

NO 

(YRS. ) 

(YRS.  ) 

(YRS.  ) 

(YRS.  ) 

(YRS.  ) 

(YRS. ) 

(YRS. ) 

1 

0.  558 

0.561 

0.561 

0.561 

0.561 

0.561  0.561 

•**; 

4. 

0.520 

0.524 

0.524 

0.524 

0.524 

0.524  0.524 

3 

0 .  0 . 

0. 

0. 

0. 

0 .  0 . 

4 

0.271 

0.377 

0.491 

0.562 

0.575 

0.575  0.575 

OTALS: 

1 . 349 

1 . 462 

1 . 576 

1.647 

1 . 660 

1 . 660  1 . 660 

Table  I.  Ultimate  Settlements* 


Elev.  Top 
Layer 

Initial 
Void  Ratio 

Consolidation 

Coefficient 

Induced 

Stress 

psf 

Final 

Stress 

psf 

Ultimate 

Settlement 

in. 

0 

1.356 

0.36 

1500 

1732 

6.43 

4 

1.201 

1500 

2126 

6.25 

10 

— 

— 

— 

— 

— 

12 

0.826 

0.24 

1500 

2594 

3.52 

18  0.787 

Total  Settlement  = 

19.60  in. 

1.63  ft 

1500 

2993 

3.40 

*  From  Means  and  Parcher 


Example  Problem  3 

67.  In  this  example  the  ultimate  settlement  and  time-settlement  rates 
resulting  from  a  A-  by  A-ft  footing  carrying  a  2000-psf  load  are  determined. 
The  soil  profile  and  soil  properties  are  presented  in  Figure  18.  Settlements 
are  computed  under  the  center  of  the  footing.  The  input  data  are  presented 
in  Figure  19. 


EL  0 


EL  -5 

v  SAND  ySAT=t25PCF 

EL  -25 

1 

ySUB  =  52.5  PCF 

CLAY  Cc  =  0.50 

PQ  =  1I50PSF 
eQ  -  1.30 

Cv*  25  F  T2/YR 

n 

MM 

SAND 

Figure  18.  Soil  profile  and  soil  properties  for  Example  Problem  3 


010  TITLE 

020  SETTLEMENT  UNDER  A  4X4  FOOTING 

0.30  3DRE  1  0  0  2000  (3-DIMENSIONAL  RECTANGULAR  LOAD  DATA  NL  1IMI  TTM2  PRES) 
<10  0  0  4  4  (XI  YI  XL  YL) 
jO  SOIL  1  0  N  125  (SOIL  DATA  N  EL  NDC  GM) 

060  SOIL  2  -5  D  52.5  .5  25  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

070  INDEX  .5  1150  1.3  (INDEX  PROPERTIES  CC  PO  EO) 

080  SOIL  3  -25  N  57.5  (SOIL  DATA  N  EL  NDC  GM) 

090  BOUS  30  (BOUSSINESQ  SOLUTION  DMAX) 

095  TIM  I  (GEOMETRIC  PROGRESSION) 

100  OUTP  2  2  4  2  5  (OUTPUT  XI  YI  X2  Y2  SP) 

110  END  (TERMINATION) 


* 


Figure  19.  Input  data  for  Example  Problem  3 
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PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  6/33  TIME:  9:11:  1 


I.  INPUT  DATA 


1.  TITLE  -  SETTLEMENT  UNDER  A  4X4  FOOTING 


2.  B0U3SINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  30.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 


LOAD  NUMBER  1  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YR 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 

PRESSURE  =  2000.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  *  0. 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  »  Y  =  4.00,  4.00 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


7.  EXCAVATION  DATA 
NONE 


8.  SOIL  DATA 


STRATA 

EL.  OF  TOP 

DRAINAGE 

EFF  UNIT 

RECGMPR . 

COEF . OF 

POISSON'S 

NO. 

OF  STRATUM 

CONDITION 

WEIGHT 

INDEX 

CONSOL. 

RATIO 

(FEET  NGVD ) 

( PCF ) 

(SQFT/YR) 

1 

■-> 

0. 

-5.  00 

N 

D 

125.00 
52.50  <: 

> . 50000  2 

5 . 00000 

0. 32000 

9.  STRESS-STRAIN  DATA 

STRATUM  NO.  1 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  2 


COMPRESSION  INDEX*  0.50000 
RECOMPRESSION  INBEX=  0.50000 
I NS ITU  VOID  RATIO*  1.30000 
I NS ITU  OVERBURDEN*  1150.00  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


A  GEOMETRIC  PROGRESSION  WITH  AN  INITIAL  TIME  PERIOD  OF  0.0192 
YEARS  AND  A  MULTIPLICATIVE  FACTOR  OF  2.0000  WILL  BE  USED  IN 
THE  TIME  RATE  OF  CONSOLIDATION  CALCULATIONS. 


11.  OUTPUT  CONTROL  DATA 


PLANE  1  XA= 

YA= 
XB= 
YB= 
DELX* 


2.0000  FT. 
2.0000  FT. 
4.0000  FT. 
2.0000  FT. 
5.0000  FT. 


WILL  OUTPUT  GO  TO  TERMINAL,  FILE,  OR  BOTH? 
ENTER  T,  F,  OR  B 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  6/33  TIME:  9:11:43 


1 1 .  OUTPUT  SUMMARY . 


1.  TITLE-  SETTLEMENT  UNDER  A  4X4  FOOTING 


POSITION:  X=  2.0  Y=  2.0 

***#**#**#•*#■#***#*#■•«■*#■**•**•»■■«•* 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 


MID-DEPTH  IN— SITU 

OF  STRATA  OVERBURDEN 
(FEET)  (LB/SQ  FT) 


DELTA  ULTIMATE 

SIGMA  SETTLEMENT 
(LB/SQ  FT)  (FEET) 


1 

r> 


2.50 

15.00 


312.50  1243.25  0. 

1150.00  103.14  0.205 


3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

ULT 

0.02 
(YRS. ) 

0.04 
(YRS.  ) 

0 . 08 
(YRS. ) 

0.  15 

(YRS.  ) 

0. 31 
(YRS. ) 

0 . 6 1 

(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.205 

0.017 

0.024 

0.032 

0.044 

0.063 

0.090 

TOTALS: 

0.205 

0.017 

0.024 

0.032 

0.044 

0.063 

0.  090 

(SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

STRATA 

1.23 

2.45 

4.91 

9.82 

19.64 

39.28 

73.55 

NO 

( YRS .  ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.  127 

0.  167 

0.  196 

0.  204 

0.205 

0.205 

0.205 

TOTALS: 

0.  127 

0.  167 

0.  196 

0.204 

0.205 

0.205 

0.205 

DO  YOU  WANT  A  PLOT  OF  STRESS  BULBS  ? 

ENTER  Y  OR  N.  j 


m 
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Example  Problem  4 


68.  In  this  example  the  ultimate  settlement  and  time-settlement  rates 
resulting  from  a  building  load  are  determined.  The  structure  is  40  by  100  ft 
in  plan  and  applies  a  uniform  load  of  2000  psf.  Settlements  are  computed 
under  the  center  of  the  building.  The  soil  profile  and  soil  properties  are 
presented  in  Figure  20.  The  input  data  for  this  example  is  shown  in  Fig¬ 
ure  21.  Results  of  this  analysis  can  be  compared  to  those  presented  by 
Winterkorn  and  Fang  (1975). 


EL  0 

P-  2000  PSF 

Til  HUM  f 

EL  -10 

EL  -20 

▼ 

ySAT‘"8PCE 

EL  -60 

Wj? 

I 

W  >”f’cr 

Cc  •  U4S 

c¥  -  S2.se  ftJ/yr 

•0-o.ra 

- 3W - mx - mm 

Figure  20.  Soil  profile  and  soil  properties 
for  Example  Problem  4 


OIO  TITLE 

020  EXAMPLE  FROM  WINTERKORN  S<  FANG 

030  30 IR  1  0  5  0  0  2000  (3-D  IRREGULAR  LOAD  DATA  NL  NS  N  TIMI  TIM2  PRES) 

040  0  O  0  40  100  40  100  0  0  0  (x(  y.  _  i  =  ljN) 

050  SOIL  1  ON  112  (SOIL  DATA  N  EL  NDC  CM) 

060  SOIL  2  -10  N  55. 5  (SOIL  DATA  N  EL  NDC  GM) 

070  SOIL  3  -20  S  59.5  .45  52.56  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

030  INDEX  .45  1970  .73  (INDEX  PROPERTIES  CC  PO  EO) 

090  SOIL  4  -30  C  59.5  .45  52.56  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

100  INDEX  .45  2550  .73  (INDEX  PROPERTIES  CC  PO  EO) 

110  SOIL  5  -40  C  59.5  .45  52.56  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

120  INDEX  .45  3430  .78  (INDEX  CC  PO  EO) 

130  SOIL  6  -60  N  60  (SOIL  DATA  N  EL  NDC  GM) 

140  BOUS  65  (BOUSSINESQ  SOLUTION  DMAX) 

150  TIMI  (GEOMETRIC  PROGRESSION) 

160  OUTP  50  20  60  20  20  (OUTPUT  XI  Y1  X2  Y2  SF) 

170  END  (TERMINATION) 


* 


Figure  21.  Input  data  for  Example  Problem  4 
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PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/83  TIME:  16:  9:22 


I.  INPUT  DATA 


1.  TITLE  -  EXAMPLE  FROM  WINTERKORN  &  FANG 


2.  BOUSSINESG  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  65.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 


LOAD  NUMBER  1 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 
ENDING  TIME  OF  APPLICATION  »  O.  YRS. 

PRESSURE  =  2000.000  PSF 


YRS. 


POINT  NO. 


1 

y 


3 

4 

5 


X 

(FT.  ) 


0. 

0. 

100. 00 
100. 00 
0. 


52 


Y 

(FT.  ) 


0. 

40.00 


40.00 

0. 


7.  EXCAVATION  DATA 
NONE 


S.  SOIL  DATA 


STRATA  EL.  OF  TOP  DRAINAGE  EFF  UNIT  RECOMPR.  COEF.OF  POISSON' 
NO.  OF  STRATUM  CONDITION  WEIGHT  INDEX  CONSOL.  RATIO 

(FEET  NGVD)  (PCF)  (SQFT/YR) 

10.  N  112.00 

2  -10.00  N  55.50 

3  -20.00  S  59.50 

4  -30.00  C  59.50 

5  -40.00  C  59.50 


9.  STRESS-STRAIN  DATA 


STRATUM  NO.  1 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  2 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  3 


COMPRESSION  INDEX*  0.45000 
RECOMPRESSION  INDEX*  0.45000 
INSITU  VOID  RATIO*  0. 7S000 
I NS ITU  OVERBURDEN*  1970.00  PSF 


STRATUM  NO.  4 


COMPRESSION  INDEX*  0.45000 
RECOMPRESSION  INDEX*  0.45000 
INSITU  VOID  RATIO*  0.78000 
INSITU  OVERBURDEN*  2550.00  PSF 


0 . 45000  52 . 56000  0 . 32000 
0.45000  52.56000  0.32000 
0 . 45000  52 . 56000  0 . 32000 


STRATUM  NO.  5 


COMPRESSION  INDEX®  0.45000 
RECOMPRESSION  INDEX®  0.45000 
I NS ITU  VOID  RATIO®  0.78000 
I NS ITU  OVERBURDEN®  3480.00  P3F 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


A  GEOMETRIC  PROGRESSION  WITH  AN  INITIAL  TIME  PERIOD  OF  0.019 
YEARS  AND  A  MULTIPLICATIVE  FACTOR  OF  2.0000  WILL  BE  USED  IN 
THE  TIME  RATE  OF  CONSOLIDATION  CALCULATIONS. 


11.  OUTPUT  CONTROL  DATA 


PLANE 


XA=  50.0000  FT. 

YA=  20.0000  FT. 

XB®  60.0000  FT. 

YB®  20.0000  FT. 

DELX®  20.0000  FT. 


WILL  OUTPUT  GO  TO  TERMINAL,  FILE,  OR  BOTH? 
ENTER  T,  F,  OR  B 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/S3  TIME:  16:11:12 

1.  TITLE-  EXAMPLE  FROM  WINTERKORN  k  FANG 


POSITION:  X=  50.0  Y=  20.0 

«***»««***«#**#*«*«*«#******* 

2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN-SITU 
OVERBURDEN 
(LB/SQ  FT) 

DELTA 
SIGMA 
(LB/SQ  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

5.00 

560.00 

1977.80 

0. 

2 

15.00 

1397.50 

1778.60 

0. 

3 

25.00 

1972.50 

1451.13 

0.607 

4 

35.00 

2567. 50 

1 148.80 

0.407 

5 

50.  00 

3460.00 

818.50 

0.473 

3.  TIME -SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

ULT 

0.02 

0.04 

0  ■  08 

0. 15 

0.31  0.61 

NO 

( YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS.)  (YRS.) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 

0.607 

0.018 

0.025 

0.035 

0.050 

0.069 

0.097 

4 

0.407 

0.013 

0.017 

0.024 

0.033 

0.046 

0.066 

5 

0.473 

0.015 

0.019 

0.027 

0.037 

0.053 

0.077 

TOTALS: 

1.487 

0.046 

0.061 

0.086 

0.  120 

0.  168 

0.240 

(SETTLEMENT 

IN  FEET 

AT  SPECIFIED 

TIMES) 

STRATA 

1.23 

2.45 

4.91 

9.82 

19.64 

39.28 

78.55 

NO 

(YRS.  ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 

0.  138 

0.  195 

0.276 

0.  386 

0.508 

0.587 

0.607 

4 

0.092 

0.  130 

0.  185 

0.  258 

0.341 

0.396 

0.407 

5 

0.  106 

0.  150 

0.214 

0.300 

0.395 

0.457 

0.473 

TOTALS* 

0.336 

0.475 

0.675 

0.944 

1.244 

1.440 

1.487 

Table  2.  Ultimate  Settlement  at  Center  of  Structure* 


Layer 

Number 

Mid-Depth 
of  Layer 
ft 

Layer 

Thickness 

ft 

In-Situ 

Stress 

psf 

Induced 

Stress 

Ultimate 

Settlement 

ft 

1 

25 

10 

1970 

1430 

0.600 

2 

35 

10 

2550 

1140 

0.377 

3 

50 

20 

3480 

810 

0.463 

Total 

Settlement  =  1 

.440  ft 

*  As  computed  in  Winterkorn  and  Fang. 


Example  Problem  5 


69.  Overview 

Example  5  was  taken  from  EM  11 10-2- 1904  (HQUSACE  1953).  The  descrip¬ 
tion  of  the  problem  is  given  below.  A  sketch  of  the  foundation  system  and  soil 
profile  is  included  as  Plate  1.  Input  data  are  shown  in  Figure  22.  Output 
is  generated  for  points  A,  B,  and  C  shown  in  Plate  1. 

70.  Computed  time-settlement  rates  and  ultimate  settlements  compare 
favorably  with  those  shown  in  this  EM.  Predicted  time-settlement  rates  are 
compared  in  Plate  12.  Ultimate  settlements  are  compared  in  Table  3. 

71.  The  following  is  quoted  from  EM  1110-2-1904: 

The  Problem.  Determine  the  total  settlement  resulting 
from  a  buried  clay  stratum  and  the  time-settlement  rate 
for  a  structure  supported  by  nine  4-  by  4-ft  footings 
located  at  an  elevation  5  ft  below  the  natural  ground 
surface.  The  foundation  plan  for  the  structure  and  the 
soil  conditions  beneath  the  structure  are  shown  in 
Plate  No.  1.  The  gross  unit  load  on  each  footing  is 
2  tons  per  square  foot.  The  construction  rate  of  load 
will  be  applied  uniformly  during  the  next  30  days.  The 
amount  of  rebound,  resulting  from  the  excavation  before 
the  construction  load  is  applied,  is  assumed  to  be 
negligible.  Consolidation  test  data  for  a  representative 
sample  of  the  clay  stratum  are  shown  on  Plates  Nos.  2  and 
8  on  which  are  plotted  pressure-void  ratio  and  time- 
settlement  data,  respectively,  for  the  sample  tested. 

Examination  of  the  consolidation  data  shows  the  clay  to 
be  normally  consolidated. 
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ENGINEERING  MANUAL  PART  CXIX  CHAPTER  4 


PLATE  NO  10 


1 


-  CSEJ T  SOLUTION 


005  TITLE 

010  EXAMPLE  PROBLEM  —  EM  1110-2-1904 

020  3DRE  1  0  .  1644  2531  (3-DIMENSIONAL  RECTANGULAR  LOAD  NL  TIMI  TIM2  PRES) 

030  0  0  4  4  (XI  Y1  XL  YL) 

040  3DRE  2  O  . 1644  2531 
050  13  0  22  4 
060  SORE  3  0  . 1644  2531 
070  36  0  40  4 
080  3DRE  4  0  . 1644  2531 
090  0  13  4  22 
100  3DRE  5  0  .1644  2531 
110  13  13  22  22 
120  SORE  6  0  . 1644  2531 
130  36  IS  40  22 
140  SORE  7  0  . 1644  2531 
150  0  36  4  40 
160  SORE  3  0  .1644  2531 
170  13  36  22  40 
ISO  3DRE  9  0  . 1644  2531 
190  36  36  40  40 
200  3DRE  10  . 1644 
210  0044 
220  3DRE  11  .1644 
230  IS  0  22  4 
240  SORE  12  .1644 
250  36  0  40  4 
260  3DRE  13  . 1644 
270  O  13  4  22 
230  3DRE  14  . 1644 
290  13  13  22  22 
300  3DRE  15  . 1644 
310  36  13  40  22 
320  3DRE  16  . 1644 
330  O  36  4  40 
340  SORE  17  .1644 
350  13  36  22  40 
360  3DRE  IS  . 1644 
370  36  36  40  40 
400  SOIL  1  -5  N  125  (SOIL  DATA  N  EL  NDC  GM) 

410  SOIL  2  -15  D  45  0  34.5  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

420  VOID  5  (VOID  RATIO  -  PRESSURE  CURVE  NP) 

430  1.15  500  1.12  1000  1.06  2000  .95  4000  .32  8000  (or.  AB.  ...  i  =  1  NP) 
440  SOIL  3  -35  N  63  (SOIL  DATA  N  EL  NDC  GM)  1  1  1 

450  BOUS  40  (BOUSSINESQ  SOLUTION  DMAX) 

460  TIMS  .0223  .0399  .205  .365  .563  .329  1.17  1.64  2.46  3. 27<SPEC1F1E» 
470  OUTPUT  0  20  20  20  20  (OUTPUT  XI  Y1  X2  Y2  SP)  TIMES  T1  T2. 

480  OUTPUT  40  0  50  0  20  (OUTPUT  XI  Y1  X2  Y2  SP) 

490  END  (TERMINATION) 


.2466  844 
.2466  844 
.2466  344 
.2466  844 
.2466  344 
.2466  344 
.2466  344 
.2466  844 
.2466  344 


Figure  22.  Input  data  for  Example  Problem  5 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/83  TIME:  15:56:20 


I.  INPUT  DATA 


1.  TITLE  -  EXAMPLE  PROBLEM  —  EM  1110-2-1904 


2.  BOUSSINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  40.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 


LOAD  NUMBER  1  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS. 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  ,  0. 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00,  4.00 


LOAD  NUMBER  2  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YRS 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS. 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  13.00,  0. 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00,  4.00 


LOAD  NUMBER 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 
ENDING  TIME  OF  APPLICATION  =  0.1644  YRS 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00.  0. 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  .  Y  =  40.00,  4.00 


LOAD  NUMBER  4  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  ,  18.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00,  22.00 


LOAD  NUMBER  5  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  18.00,  18.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00,  22.00 


LOAD  NUMBER  6  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00,  18.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  40.00,  22.00 


LOAD  NUMBER  7  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS 
PRESSURE  =  2531.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  ,  36.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00,  40.00 


1 


LOAD  NUMBER  8  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YRS.  | 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS.  t 

PRESSURE  =  2531.000  PSF  > 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT  * 

CORNER  OF  LOAD  REGION  -  X  ,  Y  =  18.00,  36.00  i 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00,  40.00 


LOAD  NUMBER  9  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 

ENDING  TIME  OF  APPLICATION  =  0.1644  YRS. 
PRESSURE  =  2531.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00,  36.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  40.00,  40.00 


LOAD  NUMBER  10  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  ,  0. 


i, 


► 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00,  4.00 


I 

I 

LOAD  NUMBER  11  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.2466  YRS.  i 

PRESSURE  =  844.000  PSF  * 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 

CORNER  OF  LOAD  REGION  -  X  ,  Y  =  18.00,  0.  i 

l 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT  j 

CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00,  4.00 


LOAD  NUMBER  12  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00,  0. 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  40.00,  4.00 


■«c*am  .j 


LOAD  NUMBER  1 3 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 
ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIN)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  .  18.00 

COORDINATE  PAIR  (PLAN  VEIN)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00.  22.00 


LOAD  NUMBER  14 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 
ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  18.00.  18.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00.  22.00 

LOAD  NUMBER  15  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 
PRESSURE  =  844.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00.  18.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  40.00.  22.00 


LOAD  NUMBER  16 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 
ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  0.  ,  36.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  4.00.  40.00 


LOAD  NUMBER  17 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 
ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  18.00.  36.00 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  22.00.  40.00 
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LOAD  NUMBER  18 


BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.1644  YRS. 
ENDING  TIME  OF  APPLICATION  =  0.2466  YRS. 

PRESSURE  =  844.000  PSF 


COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  36.00,  36.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  40.00,  40.00 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


7.  EXCAVATION  DATA 
NONE 


8.  SOIL  DATA 


STRATA  EL.  OF  TOP  DRAINAGE  EFF  UNIT  RECOMPR.  COEF.OF  POISSON'S 
NO.  OF  STRATUM  CONDITION  WEIGHT  INDEX  CONSOL.  RATIO 

(FEET  NGVD )  (PCF>  (SQFT/YR) 

1  -5.00  N  125.00 

2  -15.00  D  45.00  0.  34.50000  0.32000 


9.  STRESS-STRAIN  DATA 


STRATUM  NO.  1 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO. 


VOID  RATIO 

1.1500 
1.1200 
1 . 0600 
0.9500 
0.8200 


PRESSURE 
(PSF) 
500.0000 
1000.0000 
2000. 0000 
4000.0000 
8000. 0000 
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10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


TIME  RATE  OF  CONSOLIDATION  CALCULATIONS  WILL  BE  MADE 
AT  TIMES  <YRS>: 

0.02 
0.09 
0.21 
0.36 
0.57 
0.  S3 
1.  17 
1.64 
2.46 
3.27 


11.  OUTPUT  CONTROL  DATA 


PLANE  1  XA= 

0. 

FT 

YA= 

20 . 0000 

FT 

XB= 

20 . 0000 

FT 

YB= 

20 . 0000 

FT 

DELX  = 

20 . 0000 

FT 

PLANE  2 

XA= 

40. 0000 

FT 

YA= 

0. 

FT 

XB= 

50.0000 

FT 

YB= 

0. 

FT 

DELX  = 

20. 0000 

FT 

WILL  OUTPUT  GO  TO  TERMINAL,  FILE,  OR  BOTH? 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  5/33  TIME:  15:37:55 


II.  OUTPUT  SUMMARY. 


J 


.  -t: 


1.  TITLE-  EXAMPLE  PROBLEM 


EM  1110-2-1904 


POSITION:  X=  0.  Y=  20.0 

****#•#**■*•#■*■&***•#■************* 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 


MID-DEPTH 
OF  STRATA 
(FEET) 


IN-SITU 
OVERBURDEN 
(LB/SQ  FT) 


DELTA  ULTIMATE 

SIGMA  SETTLEMENT 
(LB/SO  FT)  (FEET) 


1 

9 


5.00  625.00  799.25  0. 

20.00  1700.00  127.89  0.071 


3.  TIME-SETTLEMENT  SUMMARY. 


STRATA 

NO 

(SETTLEMENT 
ULT  0.02 

( YRS. ) 

IN  FEET 

0.  09 
(YRS. ) 

AT  SPECIFIED 
0.21  0.36 

(YRS.)  (YRS.) 

TIMES) 
0.57 
(YRS. ) 

0.  Si 
(YRS. 

1 

0. 

0. 

0. 

0.  0. 

0. 

0. 

2 

0.071 

0. 

0 . 004 

0.013  0.024 

0.031 

0 . 039 

TOTALS: 

0.071 

0. 

0 . 004 

0.013  0.024 

0. 031 

0.  039 

( SETTLEMENT 

IN  FEET 

AT  SPECIFIED 

TIMES) 

STRATA 

NO 

1.  17 
(YRS. ) 

1.64 
(YRS. > 

2.46 
(YRS. ) 

3.27 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

2 

0.046 

0.054 

0.062 

0.066 

TOTALS: 

0.  046 

0.054 

0.062 

0.  066 

POSITION:  X=  20.0  Y=  20.0 

#  ******************************  * 

2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 

STRATA  MID-DEPTH  IN-SITU  DELTA  ULTIMATE 
NO.  OF  STRATA  OVERBURDEN  SIGMA  SETTLEMENT 
(FEET)  (LB/SQ  FT)  (LB/SQ  FT)  (FEET) 


1  5.00  625.00  1232.28 

2  20.00  1700.00  156.79 
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0. 

0.036 


r 


3.  TIME-SETTLEMENT  SUMMARY. 


STRATA 

NO 


(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 
LILT  0.02  0.09  0.21  0.36  0.57 

(YRS. )  (YRS. )  (YRS.)  (YRS. )  (YRS. > 


0.33 
(YRS.  ) 


1  0. 

2  0.086 

TOTALS:  0.0S6 


0.  0. 

0.  0.006 

0.  0.006 


0.  0. 

0.016  0.029 

0.016  0.029 


0.  0. 

0.038  0.048 

0.038  0.048 


(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 
STRATA  1.17  1.64  2.46  3.27 

NO  (YRS.)  (YRS.)  (YRS.)  (YRS.) 


1 


O.  0.  0.  0. 

0.060  0.068  0.073  0.081 


TOTALS: 


0.060  0.068  0.073  0.031 


POSITION:  X=  40.0  Y=  0. 

***************************** 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 

STRATA  MID-DEPTH  IN-SITU  DELTA  ULTIMATE 
NO.  OF  STRATA  OVERBURDEN  SIGMA  SETTLEMENT 
(FEET)  (LB/SQ  FT)  (LB/SQ  FT)  (FEET) 


1  5.00  625.00  524.05  0. 

2  20.00  1700.00  105.13  0.057 


3.  TIME-SETTLEMENT  SUMMARY. 


( SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

STRATA 

ULT 

0.02 

0.09 

0.21 

0.36 

0.57 

0.33 

NO 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS.  )' 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.057 

0. 

0.002 

0.010 

0.019 

0.027 

0.032 

TOTALS:  0.057  0 


0.002  0.010  0.019  0.027  0.032 


STRATA 

NO 


(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 
1.17  1.64  2.46  3.27 

(YRS. )  (YRS. )  <  YRS . )  (YRS. > 


1 

2 

TOTALS: 


0.  0. 

0.039  0.046 

0.039  0.046 


0.  0. 

0.051  0.055 

0.051  0.055 


Table  3.  Comparison 

of 

Ultimate  Settlements 

X 

X 

Ultimate  Settlement,  ft 

Point 

(ft) 

(ft) 

EM  110-2-1904  CSETT 

Hand  Solution 

A 

20 

20 

0.09  0.086 

B 

0 

20 

0.07  0.071 

C 

40 

0 

0.06  0.057 
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APPENDIX  A:  EXAMPLE  RUNS  WITH  HAND  VERIFICATIONS 


Example  Problem  A1 

1.  This  example  provides  a  solution  for  uniform  loads  extending  infi¬ 
nitely  in  all  directions.  Soil  stratification  and  loading  are  shown  in  Fig¬ 
ure  Al,  and  input  soil  properties  are  listed  in  Table  Al.  A  combination  of 
ramp  loads  and  an  instantaneous  load  are  applied  as  indicated  graphically  in 
Figure  A2. 

2.  The  problem  was  run  using  a  data  file  which  is  given  in  Figure  A3. 
Time-rate  calculations  are  performed  at  0.10,  0.25,  0.50,  0.80,  0.90,  1.00, 
1.20,  1.50,  1.75,  2.00,  4.00,  and  8.00  years.  Full  and  abbreviated  output 
listings  are  shown.  A  hand  solution  is  supplied  (Tables  A2-A5) ,  and  the 
results  are  compared  with  the  solution  by  CSETT  in  Table  A6. 


LOAD  ONE  p  3000  PSF 

LOAD  TWO  p  1000  PSF 

LOAD  THREE  p  2000  PSF 


EL  85.5 


EL  84.5 

7  =  115.0  PCF 

EL  81.5  ^ 

7' =52.5  PCF 

EL  81.0 

V*  57.5  PCF 

EL  78.0  f 

I 

7'=  52.5  PCF 

Figure  Al.  Soil  stratification  for  Example  Problem  Al 


3 

□X 

m 

“3- 

Al 


LOAD,  PSF 


TIME,  YRS 

Figure  A2.  Loading  sequence  for  Example  Problem  A1 

Table  A1 

Soil  Properties  for  Example  Problem  A1 


Consolidation  Data 


Stratum 

No. 

El  of  Top 
of  Stratum 
ft  NGVD 

Eff  Unit 
Weight 
pcf 

C 

2V 

ft  /yr 

e 

o 

P 

o 

psf 

CR 

cc 

1 

85.5 

115.0 

— 

— 

— 

— 

— 

2 

84.5 

52.5 

8.0 

1.6 

193.75 

0.04 

0.35 

3 

81.5 

57.5 

— 

— 

— 

— 

— 

4 

81.0 

52.5 

12.0 

1.3 

380.00 

0.045 

0.48 

5 


78.0 


57.5 


Q 

< 

o  5.000 


1.0 

TIME,  YRS 


Figure  A2.  Loading  sequence  for  Example  Problem  A1 


Table  A1 

Soil  Properties  for  Example  Problem  A 1 


El  of  Top 
of  Stratum 
ft  NGVD 


Eff  Unit 
Weight 
pcf 


115.0 


C 

2V 

ft  /yr 


Consolidation  Data 


'o  psf 


8.0  1.6  193.75  0.04  0.35 


12.0  1.3  380.00  0.045  0.48 


)  j  • 

I 


u 


I 


< 


► 


J 

i« 

J 

I 

! 

¥ 


SOLUTION  BY  CSETT 


LIST  UNI 1 


010  TITLE 

020  UNIFORMLY  DISTRIBUTED  LOADS 

030  2DPR  1  2  1.20  1.65  (2DPRESSURE  LOAD  DATA  NL  N  T1M1  T1M2) 

040  -9999.  2000  9999.  2000  (X±  Y  ...  ,  i=l,N) 

050  2DPR  22  .8  .8  (2DPRESSURE  LOAD  DATA  NL  N  T1M1  T1M2) 

060  -9999.  1000  9999.  1000  (Xi  Y  ...  ,  i=l,N) 

070  2DPR  3  2  0.0  .4  (2DPRESSURE  LOAD  DATA  NL  N  T1M1  T1M2) 

080  -9999.  3000  9999.  3000  (X.  Y  ...  ,  1=1, N) 

090  BOUS  7.5  (BOUSSINESQ  ANALYSIS  DMAX) 

100  SOIL  1  85.5  N  115.  (SOIL  DATA  N  EL  NDC  GM) 

120  SOIL  2  84.5  D  52.5  .04  8.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

130  INDEX  .35  193.75  1.6  (INDEX  PROPERTIES  CC  PO  EO) 

140  SOIL  3  81.5  N  57.5  (SOIL  DATA  N  EL  NDC  GM) 

160  SOIL  4  81.0  D  52.5  .045  12.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

170  INDEX  .48  380.0  1.3  (INDEX  PROPERTIES  CC  PO  EO) 

180  SOIL  5  78.0  N  57.5  (SOIL  DATA  IN  EL  NDC  GM) 

200  TIMS  .1  .25  .5  .8  .9  1.0  1.2  1.5  1.75  2.0  4.0  8.0  (SPECIFIED  TIMES  T1  T2) 

210  OUTP  0.0  5.0  10.0  (OUTPUT  XL  XR  SP) 

220  END  (TERMINATION) 

Figure  A3.  Input  data  file  for  Example  Problem  A1 


PROGRAM  C'SETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/S3  TIME:  10:  9:34 


I.  INPUT  DATA 


1.  TITLE  -  UNIFORMLY  DISTRIBUTED  LOADS 


2.  BOUSSINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  7.50  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 


LOAD 

NUMBER 

1  : 

NUMBER  OF  POINTS*  2 

BEGINNING  TIME  OF  APPLICATION  = 
ENDING  TIME  OF  APPLICATION  =  1. 

1.2000  YRS. 
6500  YRS, 

POINT  NO.  X 

(FT.  > 

PRESSURE 
( PSF ) 

1  -9999.00 

2  9999.00 

2000. 00 

2000 . 00 

LOAD 

NUMBER 

2  • 

NUMBER  OF  POINTS*  2 

BEGINNING  TIME  OF  APPLICATION  = 
ENDING  TIME  OF  APPLICATION  =  0. 

0.8000  YRS. 
8000  YRS. 

POINT  NO.  X 

(FT.  ) 

PRESSURE 
<  PSF ) 

1  -9999.00 

2  9999.00 

1000. 00 

1000. 00 

LOAD 

NUMBER 

3  : 

NUMBER  OF  POINTS*  2 

BEGINNING  TIME  OF  APPLICATION  = 
ENDING  TIME  OF  APPLICATION  =  0. 

0.  YRS. 

4000  YRS. 

Q 


POINT  NO 


X 

(FT.  ) 


PRESSURE 

(PSF) 


1 


-9999.00  3000.00 

9999.00  3000.00 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


7.  EXCAVATION  DATA 
NONE 


8.  SOIL  DATA 


POISSON'S 
RATIO 

'.32000 
0. 32000 

9.  STRESS-STRAIN  DATA 

STRATUM  NO.  1 
INCOMPRESSIBLE  STRATUM 


STRATA 

EL.  OF  TOP 

DRAINAGE 

EFF  UNIT 

RECOMPR 

.  COEF.OF 

NO. 

OF  STRATUM 

CONDITION 

WEIGHT 

INDEX 

CONSOL. 

(FEET  NGVD ) 

( PCF ) 

(SQFT/YR 

1 

85.50 

N 

115.00 

2 

84.50 

D 

52.50  0. 

04000 

8. 00000 

3 

81.50 

N 

57.50 

4 

81.00 

D 

52.50  0. 

04500 

12. 00000 

i 


STRATUM  NO 


o 


COMPRESSION  INDEX3  0.35000 
RECOMPRESSION  INDEX3  0.04000 
I NS ITU  VOID  RAT 10=  1.60000 
INSITU  0VERBURDEN=  193.75  PSF 


STRATUM  NO.  3 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  4 


COMPRESSION  INDEX3  0.4S000 
RECOMPRESSION  INDEX=  0.04500 
INSITU  VOID  RAT 10=  1.30000 
INSITU  OVERBURDEN3  380.00  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


TIME  RATE  OF  CONSOLIDATION  CALCULATIONS  WILL  BE  MADE 
AT  TIMES  (YRS): 

0.  10 
0.  25 
0.50 
0.80 
0.90 
1.00 
1.20 
1.50 
1.75 
2.00 
4.00 
8.00 


11. 


OUTPUT  CONTROL  DATA 


XXL=  0.  FT. 

XUL=  5.0000  FT. 

DELX=  10.0000  FT. 


DO  YOU  WANT  TO  EDIT  DATA?  ENTER  'Y'  OR  'N 


PROORAM  OSFTT  -  VFRT  T  i":Al  STRFSft  T  Nni.M":T  T  ON  AWO  SFTTI  FMFMT  PROORAM 
TiATF  :  11/17  /■?.  ?  T  T  MF  :  9 :  S7 :  39 


T  T  .  m  ITPI  IT  si  IMMARV . 


1.  TTTI  F-  UNIFORM!  Y  OTSTRTBI ITFP  I  ft  APS 


?.  SUMMARY  OF  TTMF  SFTTI  FMFNT  MATA. 


F:,l  AMF  OF  TNTFRFST:  XRT OHT=  O. 

XI  fft=  s.o 

rn-i  x  =  i n.o 


T  T  MF  X  = 

0. 

<  YR) 

1  l|.T »  1.  .STS 


0 .  1  0 

O.  IT:-  :/-. 

o ,  ss 

(") .  S09 

O  „  so 

O.  A  4  4 

0 .  so 

(") .  ASS 

0.90 

0,  SSS 

t  .  OO 

0 .  39S 

1  .  ?o 

0 .  9  1  4 

l  ..  so 

1  .  1  40 

1 . 7S 

1 . 349 

7. .  OO 

1.  .  373 

4.  OO 

1 . 37  S 

S .  00 

1 . 37 S 

nO  YOU  WANT  A  F'l  OT  OF  STRFSS  PI  II  PS  3 
FIMTFR  Y  OR  N. 

=N 


no  von  WANT  to  maff  anothfr  runs 
FNTFR  Y  OR  N 
=Y 


> 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE-*  7/12/83  TIME:  10:11:34 


II.  OUTPUT  SUMMARY. 


1.  TITLE-  UNIFORMLY  DISTRIBUTED  LOADS 


POSITION:  X=  0. 

* ■#***#■#*##*#•#*****# 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN  SITU 
OVERBURDEN 
(LB/SQ  FT  > 

DELTA 
SIGMA 
(LB/SQ  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

0.50 

57 . 50 

6000. 00 

0. 

2 

2.50 

193.75 

6000 . 00 

0 . 60S 

3 

4.25 

£;,c; 

6000.00 

0. 

4 

6. 00 

380 -  00 

6000 . 00 

0.  767 

3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

ULT 

0.  10 
(YRS. ) 

0 . 25 
(YRS. ) 

0.50 
(YRS. ) 

0  *  wo 

(YRS. > 

0.  90 
(YRS. > 

1 . 00 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.608 

0.034 

0.  127 

0.276 

0.302 

0.370 

0.  391 

3 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.767 

0.052 

0.  182 

0.368 

0.  383 

0.483 

0.  504 

TOTALS: 

1.375 

0.086 

0.  309 

0.644 

0.685 

0.  853 

0.895 

(SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

STRATA 

1.20 

1 . 50 

1.75 

2.00 

4.00 

8.00 

NO 

' YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

1 

I). 

0. 

0. 

0. 

0. 

0. 

2 

0.403 

0.501 

0.591 

0.  606 

0.608 

0.608 

3 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.511 

0.647 

0.758 

0.767 

0.767 

0.767 

TOTALS:  0.914  1.148  1.349  1.373  1.375  1.375 


A 


HAND  VERIFICATION:  SETTLEMENTS  UNDER  UNIFORMLY 
DISTRIBUTED  LOADS 


Internally  Generated  Stratification  For  Example  Problem  A1 


Input 

Stratum 

No. 


Internally 
Generated 
Stratum  No. 


Thickness  For 
Compression 
Calculations 
ft 


Eff  Thickness 
For  Consolidation 

Calculations  Drainage 

ft  Condition 


1 

2 

3 

4 

5 


1 

2 

3 

4 

5 


3.0 


3.0 


3.0 


3.0 


N 

D 

N 

D 

N 


Note:  Since  all  strata  thicknesses  are  less  than  or  equal  to  3  ft,  the 
internally  generated  stratification  is  the  same  as  the  input 
stratification. 


Table  A3 


Summary  of  Stresses  and  Ultimate  Settlements  for 
Example  Problem  A1 


Stratum 

No. 

P 

0 

psf 

e 

o 

apt  , 

Loads 

psf 

PFinal 

psf 

AH  i  . 

ult 

ft 

1 

— 

— 

— 

— 

— 

2 

193.75 

1.6 

6000.0 

6193.75 

0.608 

3 

— 

— 

— 

— 

— 

4 

380.00 

1.30 

6000.0 

6380.00 

0.767 

5 

—  — 

_ 

_ 

_ 

Calculations  for  Table  A3 

a.  Calculation  of  in  situ  overburden  (P  )  and  void  ratio 
(eo): 

Stratum  2 

P  -  1.0  x  115  +  1.5  x  52.5  *=  193.75  psf 
o 

Since  Pq  =  P0INPUT>  eQ  =  eo INPUT  1,6 
Stratum  4 

P  =  1.0  x  115  +  3.0  x  52.5  +  0.5  x  57.5  +  1.5  *  52.5  »  380. 
o 

Since  PQ  -=  P0INPUT  *  e0  “  eoiNPUT  1,3 


A12 


:;r3R 


b.  Calculation  of  ultimate  settlements 


Stratum  2 

AU  -  HCc  t  Pf inal  _  3.0(0.35)  ,  6193.75  _  ,no 

AHULT  1  +  e  log10  P  1  +  1.60  log10  193.75  ~  0,608  ft 

o  o 

Stratum  4 


_  3.0(0.48) 


6380.00 


ahult  -  nr#  logio  isoo  =  0,767  ft 


£.  Time-rate  of  consolidation  calculations: 
(1)  Time-dependent  loading 


C  t 

T  =-V 
c  h2 


where,  C  «=  coefficient  of  consolidation 
v 

t  =  time  of  interest 
tc  =  time  at  end  of  loading 

H  ■  length  of  longest  drainage  path 


for  T  >  T 


for  T  £  T 


u  ■ 1 '  N  -  ‘] 


EXP  (-M  T) 


u-(r)‘-  7  "  exp<-"2i)] 


X  100 


where  U  *  1  consolidation 


TT(2m  +  1) 


M  -  0,1,2,... 


(2)  Instantaneous  loading 


£ 


1  -  >  EXP  -M2T 

M2  C 


X  100 


M=0 


where  U  =  %  consolidation 


} 

f' 


t 


Table  A4 

Time-Rate  Calculations  for  Stratum  2 


T1M1 

T1M2 

T 

c 

Time 

yr 

Relative 

Time 

yr 

T 

1.20 

1.65 

1.60 

0.10 

0 

0 

0.25 

0 

0 

0.50 

0 

0 

0.80 

0 

0 

0.90 

0 

0 

1.0 

0 

0 

1.20 

0 

0 

1.50 

0.30 

1.07 

1.75 

0.55 

1.96 

2.0 

0.80 

2.84 

4.0 

2.80 

9.96 

0.80 

0.P0 

0.0 

0.10 

0 

0 

0.25 

0 

0 

0.50 

0 

0 

0.80 

0 

0 

0.90 

0.10 

0.36 

1.0 

0.20 

0.71 

1.2 

0.40 

1.42 

1.5 

0.70 

2.49 

1.75 

0.95 

3.38 

2.0 

1.20 

4.27 

4.0 

3.20 

11.38 

0.00 

0.40 

1.42 

0.10 

0.10 

0.36 

0.25 

0.25 

0.89 

0.50 

0.50 

1.78 

0.80 

0.80 

2.84 

0.90 

0.90 

3.20 

1.00 

1.00 

3.56 

1.20 

1.20 

4.27 

1.50 

1.50 

5.33 

1.75 

1.75 

6.22 

2.00 

2.00 

7.11 

4.00 

4.00 

14.22 

Note: 


Stratum  2 

C  ■  8.0  ft^/yr 
v 

H  -  1.5  ft 


U 

% 


0 

0 

0 

0 

0 

0 

0 

47.3 

91.6 
99.1 

100.0 

0 

0 

0 

0 

66.3 

86.0 

97.6 

99.8 
100.0 
100.0 
100.0 

il  .2 

41.7 

90.7 

99.3 

99.7 

99.9 
100.0 
100.0 
100.0 
100.0 
100.0 
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MICROCOPV  resolution  test  chart 


Table  A5 

Time-Rate  Calculations  for  Stratum  4 


Load 

T1M1 

T1M2 

T 

No. 

-1L- 

c 

1 

1.20 

1.65 

2.40 

2  0.80  0.80  0.0 


3  0.00  0.40  2.13 


Time 

yr 

Relative 

Time 

yr 

T 

0.10 

0 

0 

0.25 

0 

0 

0.50 

0 

0 

0.80 

0 

0 

0.90 

0 

0 

1.0 

0 

0 

1.20 

0 

0 

1.50 

0.30 

1.60 

1.75 

0.55 

2.93 

2.0 

0.80 

4.27 

4.0 

2.80 

14.93 

0.10 

0 

0 

0.25 

0 

0 

0.50 

0 

0 

0.80 

0 

0 

0.90 

0.10 

0.53 

1.00 

0.20 

1.067 

1.20 

0.40 

2.13 

1.50 

0.70 

3.73 

1.75 

0.95 

5.07 

2.0 

1.20 

6.40 

4.0 

3.20 

17.07 

0.10 

0.10 

0.53 

0.25 

0.25 

1.33 

0.50 

0.50 

2.67 

0.80 

0.80 

4.27 

0.90 

0.90 

4.80 

1.00 

1.00 

5.33 

1.20 

1.20 

6.40 

1.50 

1.50 

8.00 

1.75 

1.75 

9.33 

2.00 

2.00 

10.67 

4.00 

4.00 

21.3 

Note: 

Stratum  4 

C  -  12.0  ft2/yr 
v 

H  «  1.5  ft 


U 

1 


0 

0 

0 

0 

0 

0 

0 

53.0 

96.3 

99.9 

100.0 

0 

0 

0 

0 

78.3 
94.2 
99.6 

100.0 

100.0 

100.0 

100.0 

13.5 

47.5 

95.9 

99.9 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
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Example  Problem  A2 

3.  In  this  example,  settlements  are  computed  under  a  levee  and  berm 
fill.  Soil  stratification  and  load  cross  sections  are  shown  in  Figure  AA. 
Loads  are  assumed  to  be  instantaneous.  Void  ratio-log^  pressure  curves  used 
in  the  analysis  are  plotted  in  Figure  A5.  Soil  properties  are  summarized  in 
Table  A7. 

A.  The  input  data  file  used  to  generate  the  CSETT  solution  is  listed  in 
Figure  A6.  Output  is  requested  at  20-ft  intervals  extending  60  ft  from  the 
levee  centerline  in  both  directions.  A  geometric  progression  is  used  to  gen¬ 
erate  times  for  time-rate  calculations  using  the  default  parameters. 

5.  A  hand  solution  of  settlement  at  the  levee  centerline  is  provided  for 
comparison  (Tables  A8-A10) .  Results  of  the  hand  solution  and  the  CSETT  solu¬ 
tion  are  compared  in  Table  All. 


DISTANCE.  FT 

-60  -40  -20  0  20  40  60 


Figure  AA.  Soil  stratification  ard  loads  for  Example  Problem  A2 
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SOLUTION  BY  CSETT 


*LIST  LEVE 


010  TITLE 

020  SETTLEMENT  UNDER  LEVEE  SECTION 

030  2DS0  1  8  0.0  0.0  105.0  (2DSOIL  LOAD  DATA  NL  N  T1ML  TLM2  GAM) 

040  -9999.  100.0  -45.0  100.0  -7.5  112.5  7.5  112.5  20.0  108.33 

050  46.67  107.0  67.67  100.0  9999.  100.0  (Xi  Y  ....  1=1,  N) 

060  2DS0  2  6  0.0  0.0  115.0  (2DS0IL  LOAD  DATA  NL  N  TLML  TLM2  GAM) 

070  -9999.  100.0  -45.0  100.0  -7.5  112.5  7.5  112.5  45.0  100.0 

080  9999.  100.0  (X±  ...  i=l,  N) 

090  SOIL  1  100.0  N  120.0  (SOIL  DATA  N  EL  NDC  GAM) 

100  SOIL  2  98.5  S  52.5  .07  20.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

110  VOID  2  (VOID  RATIO  PRESSURE  CURVE  NP) 

120  1.6  200.0  .50  20000.0  (ORi  AB£  ....  i=l,  NP) 

130  SOIL  3  95.0  S  52.5  .05  10.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

140  VOID  2  (VOID  RATIO  PRESSURE  CURVE  NP) 

150  1.3  200.0  .40  20000.0  (OR.^  AB^  ....  1=1,  NP) 

160  SOIL  4  87.5  N  57.5  (SOIL  DATA  N  EL  NDC  GM) 

170  SOIL  5  82.5  D  52.5  .04  16.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

180  VOID  2  (VOID  RATIO  PRESSURE  CURVE  NP) 

190  1.2  200.0  .35  20000.0  (ORi  AB^  ....  1=1,  NP) 

200  SOIL  6  75.0  N  57.5  (SOIL  DATA  N  EL  NDC  GM) 

210  BOUS  25.0  (BOUSSINESQ  DMAX) 

220  TIMI  (GEOMETRIC  PROGRESSION  W/DEFAULT  VALUES) 

230  OUTP  -60  60  20  (OUTPUT  XL  XR  SP) 

240  END  (TERMINATION) 

Figure  A6.  Input  data  file  for  Example  Problem  A2 
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TS  T  NPI  IT  FROM  TFRMTNAt  OR  A  FT  IF 
FNTFR  T  OR  F 
=F 


FNTFR  OAT  A  F  TI..F  NAMF 
=1  FVF 


T  NPI  IT  COMPIFTF.  Pfl  YOU  WANT  T  NPI  IT  PAT  A 
FPHOPRTNTFP  TO  YOUR  TFRMTNAI  A  FTIF, 
ROTH-  OR  NFTTHFRv  ( FIvlTFR  T,  F-  R,  OR  N> 
=T 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/83  TIME:  9:41:55 

I.  INPUT  DATA 


1.  TITLE  -  SETTLEMENT  UNDER  LEVEE  SECTION 


2.  BOUSSINESG  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  25.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


* 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 


PROFILE  NUMBER  1  : NUMBER  OF  POINTS=  3 

BEGINNING  TIME  OF  APPLICATION  =  O.  YRS. 

ENDING  TIME  OF  APPLICATION  =  O.  YRS. 

EFFECTIVE  UNIT  WEIGHT  OF  SOIL  LOAD®  105.00  PCF 


POINT 

NO.  X 

Y 

(FT.  ) 

(FT.  ) 

• 

1 

-9999. 00 

100. 00 

'  i 

2 

-45.00 

100.00 

\* 

3 

-7.50 

112.50 

| 

4 

7.50 

112.50 

1  , 

i 

5 

20.00 

108.33 

V 

6 

46.67 

1 07 . 00 

7 

67.67 

1 00 . 00 

3 

9999. 00 

100. 00 

PROFILE  NUMBER  2  : NUMBER  OF  POINTS®  6 

BEGINNING  TIME  OF  APPLICATION  =  0.  YRS. 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 

EFFECTIVE  UNIT  WEIGHT  OF  SOIL  LOAD®  115.00  PCF 


POINT  NO.  X  Y 

(FT.)  (FT.) 

1  -9999.00  100.00 

2  -45.00  100.00 

3  -7.50  112.50 

4  7.50  112.50 

5  45.00  100.00 

6  9999.00  100.00 


S.  SOIL  DATA 


STRATA 

EL.  OF  TOP 

DRAINAGE 

EFF  UNIT  RECOMPR 

.  COEF.OF 

POISSON'S 

NO. 

OF  STRATUM 

CONDITION 

WEIGHT 

INDEX 

CONSOL. 

RATIO 

(FEET  NGVD) 

(PCF ) 

(SQFT/YR) 

1 

100.00 

N 

120.00 

2 

98.50 

s 

52.50 

0.07000 

20.00000 

0.32000 

3 

95.00 

S 

52.50 

0.05000 

10.00000 

0. 32000 

4 

87.50 

N 

57.50 

5 

82.50 

D 

52.50 

0.04000 

16.00000 

0.32000 

9.  STRESS-STRAIN  DATA 


STRATUM  NO.  1 


INCOMPRESSIBLE 

:  STRATUM 

STRATUM 

NO. 

2 

VOID  RATIO 

PRESSURE 

(PSF) 

1 . 6000 

200.0000 

0.5000 

20000.0000 

STRATUM 

NO. 

3 

VOID  RATIO 

PRESSURE 

(PSF) 

1 . 3000 

200.0000 

0.4000 

20000.0000 

STRATUM 

NO. 

4 

INCOMPRESSIBLE  STRATUM 


STRATUM 

NO.  5 

VOID  RATIO 

PRESSURE 

(PSF) 

1 . 2000 

200.0000 

0.3500 

20000. 0000 

10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


A  GEOMETRIC  PROGRESSION  WITH  AN  INITIAL  TIME  PERIOD  OF  0.019 
YEARS  AND  A  MULTIPLICATIVE  FACTOR  OF  2.0000  WILL  BE  USED  IN 
THE  TIME  RATE  OF  CONSOLIDATION  CALCULATIONS. 


11.  OUTPUT  CONTROL  DATA 


XXL*  0.  FT. 

XUL=  5.0000  FT. 

DELX*  10.0000  FT. 


DO  YOU  WANT  TO  EDIT  DATA?  ENTER  'V'  OR  'N'. 
=  N 


PROGRAM  CSFTT  -  VFRTTOAI 
OATF:  11/17/3? 


.-  -j-RF Thinl  1GTTOIJ  AND  •-.FTTI  FMFNT  PROGRAM 

ttmf:  9:  54:  ?/-. 


T  T  .  01 JTPI  IT  SUMMARY 


1.  TITI.F-  SFTTI  FMFNT  IJNDFR  I.FVFF  SECT T ON 


•2.  SUMMARY  OF  TTMF  SFTTI  FMFNT  DATA. 


PI  AMF  OF  T  MTFRFST :  <RTGHT=  O. 

VI  FFT=  5,0 
DPI  x=  10.0 


! 
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T  T  MF  y  = 
( YR  ) 


I II.  T  . 

o.  02 
0 . 04 

O .  O:-: 
o.  1  s 
o, 9  1 
O .  A 1 

1 . 23 

2.  4  s-', 
4.91 
9  .  .99 
19.  A 4 
39 . 90 
79,55 


1 . 995 
O.  265 
O.  >:7  4 
O .  599 
0.747 
1  ..  020 
1 . 329 
1  .  563 
1 . 759 
1.91  ? 
I  .  905 
1 . 9  9  Ft 

1 . 995 

1 . 995 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/83  TIME:  9:43:  3 


II.  OUTPUT  SUMMARY. 


1.  TITLE-  SETTLEMENT  UNDER  LEVEE  SECTION 


POSITION:  X=  0. 

*«*«***««#«'«'*#*«*«"» 


SUMMARY 

STRATA 

NO. 

OF  ULTIMATE 

MID-DEPTH 
OF  STRATA 
(FEET) 

SETTLEMENTS. 

IN— SITU 
OVERBURDEN 
(LB/SQ  FT) 

DELTA 
SIGMA 
(LB/SQ  FT) 

ULTIMATE 
SETTLEMENT 
<  FEET ) 

1 

0.75 

90.  00 

1437.40 

0. 

2 

3.25 

271.88 

1432.65 

0.610 

3 

8.75 

560.63 

1385.40 

0.879 

4 

15.00 

901.25 

1299. 10 

0. 

5 

21.25 

1241.88 

1204.80 

0.  506 

A26 


1 


3.  TIME-SETTLEMENT  SUMMARY. 


STRATA 

NO 

(SETTLEMENT 
ULT  0.02 

( YRS. ) 

IN  FEET  AT  SPECIFIED 
0.04  0.08  0.15 

(YRS.)  (YRS.)  (YRS.) 

TIMES) 
0.31 
(YRS. ) 

0.  61 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

O 

A- 

0.  610 

0.  122 

0.  173 

0.  244 

0.344 

0.  467 

0.  569 

o 

0.879 

0.058 

0.082 

0.  1 16 

0.  164 

0.  231 

0.  327 

4 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

5 

0.  506 

0 . 085 

0.  1 19 

0.  168 

0 . 239 

0 . 332 

0.433 

TOTALS: 

1 . 995 

0.265 

0 . 374 

0 . 528 

0.747 

1 . 030 

1 . 329 

STRATA 

NO 

(SETTLEMENT 
1.23  2.45 

(YRS.)  (YRS.) 

IN  FEET  AT  Si 
4.91  9.8 

(YRS.)  (YRS. 

PEC IF I ED 
2  19. 64 
)  (YRS.) 

TIMES) 

39.  28 
( YRS . ) 

78.55 
(YRS. > 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

Am 

0.  608 

0.  610 

0.610 

0 . 6 1 0 

0 .  6 1 0 

0 .610 

0 . 6 1 0 

o 

0.  462 

O  •  636 

0.  796 

0 . 869 

0. 879 

0 . 879 

0.879 

4 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

5 

0 . 493 

0 . 506 

0.506 

0 . 506 

0.506 

0 . 506 

0.506 

TOTALS: 

1 . 563 

1 . 75 2 

1.912 

1 . 985 

1 . 995 

1 . 995 

1 . 995 

5 
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HAND  VERIFICATION:  SETTLEMENTS  UNDER  £  LEVEE  SECTION 


Table  A8 

Internally  Generated  Stratification  for  Example  Problem  A2 


Thickness  for 

Eff  Thickness 

Input 

Internally 

Compression 

for  Consolidation 

Stratum 

Generated 

Calculations 

Calculations 

Drainage 

No. 

Stratum  No. 

ft 

ft 

Condition 

1 

2 

— 

— 

N 

3 

1.75 

3.5 

S 

2 

4 

1.75 

3.5 

S 

5 

2.50 

7.5 

s 

3 

6 

2.50 

7.5 

s 

7 

2.50 

7.5 

s 

8 

_ 

N 

4 

9 

— 

— 

N 

10 

2.50 

7.5 

D 

5 

11 

2.50 

7.5 

D 

12 

2.50 

7.5 

D 

Note:  The  stratification  generated  by  CSETT  will  be  used  in  this  solution. 


H 


Table  A9 

Summary  of  Stresses  and  Settlements  for  Example  Problem  A2 


Stratum 

P 

o 

AP  . 

Loads 

^Final 

ahult 

No. 

psf 

o 

psf 

psf 

ft 

2 

— 

— 

— 

— 

— 

3 

225.94 

1.571 

1436 

1662 

0.324 

4 

317.81 

1.489 

1430 

1748 

0.286 

5 

429.38 

1.151 

1414 

1843 

0.331 

6 

560.63 

1.099 

1388 

1949 

0.290 

7 

691.88 

1.057 

1354 

2046 

0.258 

8 

— 

— 

— 

— 

— 

9 

— 

— 

— 

— 

— 

10 

1110.63 

0.884 

1242 

2353 

0.184 

11 

1241.88 

0.863 

1204 

2446 

0. 168 

12 

1373.13 

0.844 

1168 

2541 

0.154 

Note : 

P  =  in  situ 
o 

overburden 

(see  previous 

calculations) 

e  =  In  situ  void  ratio  (see  previous  calculations) 
o 

AP,  .  =  induced  stresses  (from  CSETT) 

Loads 

P  =  p  +  AP 

Final  o  Loads 

AH^Li  =  Ultimate  settlements  (see  following  computations) 


Calculations  for  Table  A9 


a.  Computation  of  in  situ  overburden  (P  )  and  void  ratio  (e^) 


Stratum  3 

P  =  1.5  x  120.0  +  x  52.5  =  225.94  psf 

O  Z  •  u 


e 

o 


log 


P 

o 

10  P 


C 

c 


e  -  e 
1  2 


log10  p2 


log10  P1 


_ 1.60  -  0.50 

log1()  20,000  -  log1Q  200 


=  0.550 


e 

o 


1.60  -  0.55  log 


225.94 

10  200.0 


1.571 


Stratum  4 


P  =  225.94  +  1.75  (57.5)  =  317.81  psf 
o 

eo  -  1.6  -  0.55  l»i1()  ^  ■  1.489 


Stratum  5 


2  SO 

P  =  1.5  x  120  +  3.5  x  52.5  +  x  52.5  =  *29. 38  psf 

o  z .  'J 


C  = 


1.30  -  0.40 


c  l°g j o  20000  -  loglQ  200 


=  0.450 


eo  1.30  -  0.45  log.oHH4-  1.151 


Stratum  6 


F  =  429.38  +  2.5(52.5)  =  560.63  psf 
o 


«„  ■  ‘-3  -  °-45  ^lolnr*  >•»” 


Stratum  7 


P  -  560.63  +  2.5(52.5)  =  691.88  psf 
o 
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«o  ■ l-3  -  °-450  ‘'‘nif  ■  '-057 


Stratum  10 


P  -  1.5  x  120  +  11.0  x  52.5  +  5.0  *  57.5  +  1.25  x  52.5 
o 

r  =  1.20  -  0.35 

c  log10  20000  -  log1()  200  U*  3 


„  1  ->n  1  1110.63  _  oa/ 

o  =  1,20  ”  log10  200.0  *  °*884 


=  1110.63  psf 


Stratum  1 1 


P  =  110.63  +  2.5(52.5)  =  1241.88  psf 
o 


%  ■  ‘-2°  -  °-425  loeio  ifrotlr  •  °-863 


Stratum  12 


Pq  -  1241.88  +  2.5(52.5)  =  1373.13  psf 


eo  •  >-20  -  °-425  ^io  Tsttr  ■  °-844 


b.  Computation  of  ultimate  settlements: 


Stratum  3 


au  c  .  final  1.75(0.550)  .  1662  . 

AHULT  1  +  e  log10  P  1  +  1.571  log10  225.74  =  0,32A  ft 


1  +  1.571  e10  225.74 


Stratum  4 


ah  .  1*75(0.55)  .  1748  _ 

AHULT  1  +  1.489  1  810  317.8  °*286  ft 


Stratum  5 


AH  -  —  50(0.45)  j  B  o  331  ft 

ULT  1  +  1.151  810  429.4 


Stratum  6 


Au  2.5(0.45)  . 

AHULT  “  1  +  1.099  108 


10  560.6 


-  0.290  ft 


Stratum  7 


AH  — 

2. 

5(0.45) 

ahult 

1 

+  1.057 

Stratum 

10 

2. 

5(0.425) 

ahult 

1 

+  0.884 

Stratum 

11 

ahult“ 

2. 

5(0.425) 

1 

+  0.863 

Stratum 

12 

AH  m 

2. 

5(0.425) 

A  uLT 

1 

+  0.844 

log 


2046 


10  691.9 


log 


2353 

10  1110.6 


log 


2446 
10  1242 


log 


2541 
10  1373 


=  0 
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.258  ft 


0.184  ft 


.168  ft 


.154  ft 


Table  A10 


Time-Rate  of  Consolidation  Calculations  for  Example  Problem  A2 


Stratum* 

AH 

ULT** 

H 

C 

t 

U 

Settlement 

No. 

ft 

ft 

f t  /yr 

_  Xr 

T 

% 

ft 

3&4 

0.610 

3.5 

20.0 

0.019 

0.031 

19.9 

0.121 

(2) 

0.077 

0.126 

40.0 

0.244 

0.307 

0.501 

76.5 

0.466 

1.227 

2.003 

99.4 

0.606 

4.910 

8.016 

100.0 

0.610 

19.638 

32.06 

100.0 

0.610 

5,6,&  7 

0.879 

7.5 

10.0 

0.019 

0.003 

6.6 

0.058 

(3) 

0.077 

0.014 

13.2 

0.116 

0.307 

0.055 

26.4 

0.232 

1.227 

0.218 

52.6 

0.462 

4.910 

0.873 

90.6 

0.796 

19.638 

3.491 

100.0 

0.879 

10,11, 

0.506 

3.75 

16.0 

0.019 

0.022 

16.6 

0.084 

&  12 

0.077 

0.088 

33.4 

0. 169 

(5) 

0.307 

0.349 

65.8 

0.333 

1.227 

1.396 

97.4 

0.493 

4.910 

5.586 

100.0 

0.506 

19.638 

22.34 

100.0 

0.506 

*  Numbers  in  parentheses  i: licate  input  stratification. 

**  Total  settlement  for  indicated  stratum. 

Note:  For  instantaneous  loading  and  a  linear  initial  distribution  of  excess 
pore  pressure: 


U%  = 


x  100 


where 


Table  All 

Summary  of  Time-Settlement  Data  for  Example  Problem  A2* 


ahult 

ft 

Settlement,  ft. 

At  Indicated 

Time 

Stratum. 

No. 

0.019 

yr 

0.077 

yr 

0.307 

yr 

1.227 

_ZI _ 

4.910 

yr 

19.638 

_ Y1 _ 

2 

(0.610) 

0.610 

(0.122) 

0.121 

(0.244) 

0.244 

(0.467) 

0.466 

(0.608) 

0.606 

(0.610) 

0.610 

(0.6101 

0.610 

3 

(0.879) 

0.879 

(0.058) 

0.058 

(0.116) 

0.116 

(0.231) 

0.232 

(0.462) 

0.462 

(0.796) 

0.796 

(0.879) 

0.879 

5 

(0.506) 

0.506 

(0.085) 

0.084 

(0.168) 

0.169 

(0.332) 

0.333 

(0.493) 

0.493 

(0.506) 

0.506 

(0.506) 

0.506 

Total 

(1.995) 

1.995 

(0.265) 

0.263 

(0.528) 

0.529 

(1.030) 

1.031 

(1.563) 

1.561 

(1.912) 

1.912 

(1.995) 

1.995 

*  Results  of  the  CSETT  solution  are  shown  in  parentheses  for  comparison. 


6.  This  example  is  provided  to  illustrate  the  analysis  of  irregular  two- 
dimensional  loads  combined  with  a  construction  excavation.  Soil  stratifica¬ 
tion  and  load  cross  sections  are  shown  in  Figure  A7.  The  load  due  to  a  U-frame 
lock  system  is  shown  in  Figures  A8  and  A9.  Input  soil  properties  are  listed 

in  Table  A12.  The  void  ratio-log^  pressure  curve  for  stratum  1  is  plotted  in 
Figure  A10.  Ramp  loads  are  used  to  simulate  the  construction  and  backfill 
sequence  as  shown  in  Figure  All. 

7.  The  input  data  file  used  to  generate  the  CSETT  solution  is  listed  in 
Figure  A12.  Output  is  generated  at  30-ft  intervals  extending  120  ft  either 
side  of  the  centerline  of  the  U-frame  section.  Time-rate  calculations  are 
made  at  0.10,  0.25,  0.50,  1.00,  2.00,  2.25,  2.50,  5.00,  and  10.00  years. 

8.  A  hand  solution  of  settlements  under  the  centerline  of  the  U-frame  is 
provided  for  comparison  (Tables  A13-A17).  Results  of  the  hand  solution  and 
the  CSETT  solution  are  compared  in  Table  A18. 


Figure  A7.  Soil  stratification  and  loads  for 
Example  Problem  A3  (not  to  scale) 
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Table  A12 

Soil  Properties  for  Example  Problem  A3 


Consolidation  Data 


Stratum 

No. 

El  of  Top 
of  Stratum 
ft  NGVD 

Eff  Unit 
Weight 
pcf 

C 

ft2/yr 

e 

o 

P 

o 

psf 

c 

r 

c 

c 

1 

41.0 

52.5 

17.0 

: k 

* 

0.110 

* 

2 

35.0 

57.5 

— 

— 

— 

— 

— 

3 

34.0 

52.5 

16.0 

1.40 

4917.50 

0.072 

0.52 

4 

30.0 

57.5 

— 

— 

— 

— 

— 

5 

28.0 

52.5 

11.0 

1.60 

5268.75 

0.040 

0.41 

6 

23.0 

52.5 

11.0 

1.51 

5557.50 

0.035 

0.35 

7 

17.0 

52.5 

11.0 

1.30 

5872.50 

0.028 

0.30 

8 

11.0 

57.5 

_ 

_ 

__ 

*  Input  curve  (see  Figure  A10) . 
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VOID  RATIO, 


COEFFICIENT  OF  PERMEABILITY,  k2Q,  10'  CM/SEC 
0.1  0.2  0.3  0.40.5  1  2  3  4  5  10 


PRESSURE,  p,  T/SQ  FT 

Figure  A10.  Laboratory  compression  curve  for  stratum  1 


TIME,  YRS 


Figure  All.  Loading  sequence  for  Example  Problem  A3 


SOLUTION  BY  CSETT 


10  TITLE 

20  CAL ION  LOCK  SETTLEMENTS 

30  2DPR  1  10  0  2  (2DPRESSURE  LOAD  DATA  NL  N  TIM1  TIM2) 

40  -9999  0  -63  0  -63  3604.5  -42  3604.5  -42  1802.25  (Xi,  Yi _ ,  i  = 

50  42  1802.25  42  3604.5  63  3604.5  63  0  9999  0 

60  2DS0  2  14  2  2.2  130  (2DS0IL  LOAD  DATA  NL  N  TIM1  TIM2  GAM) 

70  -9999  41  -211  41  -100  78  -45  78  -47  63 
80  -63  61  -63  41  63  41  63  61 

90  47  63  45  78  190  78  116  41  9999  41  Yi 

100  EXCA  4  120  (EXCAVATION  DATA  N  GAM) 

110  -9999  41  116  41  190  78  9999  78  (Xi  Yi _  i  =  1,N) 

120  SOIL  1  41  D  52.5  .11  17  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

130  VOID  8  (VOID  RATIO  -  PRESSURE  CURVE  DATA  NP) 

140  1.44  30  1.4  500  1.39  1000  1.38  2000 

150  1.3  4000  1.15  7000  .9  16000  .82  32000  (OR.  AB . ,  i  =  1,NP) 

160  SOIL  2  35  N  57.5  (SOIL  DATA  N  EL  NDC  GM  CR  CV^T)1 

170  SOIL  3  34  D  52.5  .072  16  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

180  INDEX  .52  4917.5  1.4  (INDEX  PROPERTIES  CC  PO  EO) 

190  SOIL  4  30  N  57.5  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

200  SOIL  5  28  D  52.5  .04  11  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

210  INDEX  .  4 1  5268.75  1.6  (INDEX  PROPERTIES  CC  PO  EO) 

220  SOIL  6  23  C  52.5  .035  11  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

230  INDEX  .35  5557.5  1.51  (INDEX  PROPERTIES  CC  PO  EO) 

240  SOIL  7  17  C  52.5  .028  11  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

250  INDEX  .3  5872.5  1.3  (INDEX  PROPERTIES  CC  PO  EO) 

260  SOIL  8  11  N  57.5  (SOIL  DATA  N  EL  NDC  GM  CR  CV  ST) 

270  BOIJS  35  (BOUSSINESQ  DMAX) 

280  TIMS  .1  .25  .5  1.  2.  2.25  2.5  5.  10.  (SPECIFIED  TIMES  T1  T2  T3 _ ) 

290  OIJTP  0  5  10  (OUTPUT  XL  XR  SP) 

300  END  (TERMINATION) 


i  =  1 » N) 


Figure  A12.  Input  data  file  for  Example  Problem  A3 


1  • 


Ts  TNPIIT  FROM  TFRMTWAI  HR  A  FTI  F 
FNTFR  T  HR  F 
=F 


FNTFR  DATA  FT  IF  NANF 

=r.Ai  ton 


TIMF'IIT  I'flMPI  FTF.  Hfl  YhI  I  U1AIJT  ]  NPl  IT  TiATA 
FTHTlPR  T  WTFTi  TO  Ynl  |R  TFRI-NNAI  .  «  FIIF, 

RfiTH ,  flR  NFTTHFR’  (FNTFR  T,  F,  p.  I  iR  N ) 
=  T 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  6/ S3  TIME:  9:  1:60 


I.  INPUT  DATA 


1.  TITLE  -  CAL  I ON  LOCK  SETTLEMENTS 


2.  BOUSSINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  35.00  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 


LOAD  NUMBER  1  :  NUMBER  OF  POINTS-  10 

BEGINNING  TIME  OF  APPLICATION  =  0.  YRS. 

ENDING  TIME  OF  APPLICATION  =  2.0000  YRS. 


POINT 

NO.  X 

PRESSURE 

(FT.  ) 

<  PSF ) 

1 

-9999.00 

0. 

2 

-63. 00 

0. 

3 

-63.00 

3604.50 

4 

-42.00 

3604.50 

5 

-42.00 

1802.25 

6 

42.00 

1802.25 

7 

42.00 

3604.50 

S 

63.00 

3604.50 

9 

63.00 

0. 

10 

9999.00 

0. 
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**£  : 


« 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 


PROFILE  NUMBER  2  : NUMBER  OF  POINTS=  14 

BEGINNING  TIME  OF  APPLICATION  =  2.0000  YRS. 
ENDING  TIME  OF  APPLICATION  =  2.2000  YRS. 
EFFECTIVE  UNIT  WEIGHT  OF  SOIL  LOAD*  130.00  PCF 


POINT  NO.  X  Y 


(FT.  ) 

(FT.  ) 

1 

-9999. 00 

4  1 . 00 

9 

-211.00 

4 1 . 00 

9 

- 1 00 . 00 

78.00 

4 

-45.00 

73.00 

5 

-47.00 

63 . 00 

6 

-63. 00 

o 

o 

.■  A 

7 

-63. 00 

4  1 . 00 

o 

63. 00 

4 1 . 00 

9 

63.00 

61.00 

10 

47.  00 

63.  00 

1  1 

45.00 

78 . 00 

12 

190. 00 

78.  00 

13 

116. 00 

4 1 . 00 

14 

9999. 00 

4 1 . 00 

5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


7.  EXCAVATION  DATA 


NUMBER  OF  POINTS  DEFINING  EXCAVATED  PROFILE*  4 
EFFECTIVE  UNIT  WEIGHT  OF  EXCAVATED  MATERIAL*  120.00  PCF 


POINT  NO.  X  Y 

(FT. )  (FT.) 

1  -9999.00  41.00 

2  116.00  41.00 

3  190.00  73.00 

4  9999.00  73.00 
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3.  SOIL  DATA 


STRATA 

EL.  OF  TOP 

DRAINAGE 

EFF  UNIT 

RECOMPR 

.  COEF.OF 

POISSON 

NO. 

OF  STRATUM 

CONDITION 

WEIGHT 

INDEX 

CONSOL. 

RATIO 

(FEET  NGVD) 

( PCF ) 

(SQFT/YR) 

1 

41.00 

D 

52.50  0. 

1 1 000 

17.00000 

0.32000 

2 

35.  00 

N 

57 . 50 

3 

34.00 

D 

52.50  0. 

07200 

1 6 . 00000 

0. 32000 

4 

30 . 00 

N 

57 . 50 

5 

23.00 

D 

52 . 50  0 . 

04000 

1 1 . 00000 

0. 32000 

6 

23.00 

c 

52.50  0. 

03500 

1 1 . 00000 

0 . 32000 

7 

17.00 

c 

52.50  0. 

02800 

1 1 . 00000 

0.  32000 

9.  STRESS-STRAIN  DATA 


STRATUM  NO.  1 


VOID  RATIO 

PRESSURE 

(PSF) 

1.4400 

30.0000 

1 . 4000 

500.0000 

1 . 3900 

1000. 0000 

1 . 3300 

2000. 0000 

1 . 3000 

4000.0000 

1. 1500 

7000.0000 

0. 9000 

16000. 0000 

0.3200 

32000 . 0000 

STRATUM  NO. 

2 

INCOMPRESSIBLE 

STRATUM 

STRATUM  NO. 

3 

COMPRESSION  INDEX=  0.52000 
RECOMPRESSION  INDEX=  0.07200 

I NS ITU  VOID  RAT 10=  1.40000 
I NS ITU  0VERBURDEN=  4917.50 

STRATUM  NO. 

4 

INCOMPRESSIBLE  STRATUM 


A44 


STRATUM  NO.  5 


COMPRESSION  INDEX=  0.41000 
RECOMPRESSION  INDEX=  0.04000 
I NS ITU  VOID  RAT 10=  1.60000 
I NS ITU  uVERBURDEN=  5263.75  PSF 


STRATUM  NO.  6 


COMPRESSION  INDEX=  0.35000 
RECOMPRESSION  INDEX=  0.03500 
INSITU  VOID  RAT 10=  1.51000 
I NS ITU  OVERBURDEN=  5557.50  PSF 


STRATUM  NO.  7 


COMPRESSION  INDEX=  0.30000 
RECOMPRESSION  INDEX=  0.02300 
INSITU  VOID  RAT 10=  1.30000 
INSITU  OVERBURDEN=  5372.50  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


TIME  RATE  OF  CONSOLIDATION  CALCULATIONS  WILL  BE  MADE 
AT  TIMES  < YRS ) : 

0.  10 
0.25 
0.50 
1.00 
2.00 
2.25 
2.50 
5.  00 
10.00 


11.  OUTPUT  CONTROL  DATA 


XXL=  0.  FT. 

XUL=  5.0000  FT. 

DELX=  10.0000  FT. 


WILL  OUTPUT  GO  TO  TERMINAL,  FILE,  OR  BOTH? 
ENTER  T,  F,  OR  B 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/  6/83  TIME:  9:10:  8 


II.  OUTPUT  SUMMARY. 


1.  TITLE-  CAL I ON  LOCK  SETTLEMENTS 


POSITION:  X=  0. 


SUMMARY 

OF  ULTIMATE 

SETTLEMENTS. 

STRATA 

MID-DEPTH 

IN— SITU 

DELTA 

ULTIMATE 

NO. 

OF  STRATA 

OVERBURDEN 

SIGMA 

SETTLEMENT 

(FEET) 

(LB/SQ  FT) 

(LB/SC!  FT) 

(FEET) 

1 

3.  00 

4597 . 50 

1803. 05 

0.314  » 

~f 

6.50 

4783.75 

1806. 80 

0 . 

9 . 00 

4917.50 

1814.50 

0 . 080 

4 

12.00 

5080 . 00 

1829.20 

0. 

5 

15.50 

5268. 75 

1857. 75 

0.03? 

6 

21.00 

5557.50 

1924.25 

0 . 036 

7 

27.00 

5872.50 

2023.75 

0.027 

3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

ULT 

0.  10 
<  YRS.  ) 

0.25 
(YRS. ) 

0.  50 
(YRS. ) 

1 . 00 

(YRS. ) 

2.  00 
(YRS. ) 

2  m  2! 
(YRS. 

1 

0.314 

0.  005 

0 . 020 

0 . 053 

0.  129 

0. 286 

0 . 305 

2 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 

0.  080 

0 . 002 

0.  007 

0.017 

0.  036 

0.  076 

0.  080 

4 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

5 

0 . 039 

0. 

0. 

0 . 002 

0 . 006 

0.016 

0.019 

6 

0.036 

0. 

0. 

0.002 

0.005 

0.014 

0.  017 

7 

0.  027 

0. 

0. 

0 . 002 

0 . 004 

0 . 0 1 0 

0 . 0 1 2 

TOTALS: 

0.496 

0.007 

0.027 

0.076 

0.  180 

0.  402 

0.  433 
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(SETTLEMENT  IN  FEET  AT 


STRATA 

NO 

2.50 
(YRS. ) 

5 . 00 
(YRS. ) 

1 0 .  00 
(YRS. > 

1 

0.311 

0.314 

0.314 

2 

0. 

0. 

0. 

3 

0 . 080 

0 . 080 

0.  080 

4 

0. 

0. 

0. 

5 

0.  021 

0 . 03 1 

0.038 

6 

0.019 

0.  029 

0 . 035 

7 

0.  014 

0. 022 

0.  027 

TOTALS: 

0.445 

0.  47/-. 

0.  494 

DO  YOU  WANT  A  PLOT  OF  STRESS  BULB 
ENTER  Y  OR  N. 

=N 


DO  YOU  WANT  TO  MAKE  ANOTHER  RIJN? 
ENTER  Y  OR  N 


HAND  SOLUTION  OF  SETTLEMENTS  UNDER  £ 
OF  U-FRAME 


Table  A13 

Internally  Generated  Stratification  for  Example  Problem  A3 


Input 

Stratum 

No. 

Internally 
Generated 
Stratum  No. 

Thickness  for 
Compression 
Calculations 
ft 

Eff  Thickness 
for  Consolidation 
Calculations 
ft 

Drainage 

Condition 

2 

3.0 

6.0 

D 

1 

3 

3.0 

6.0 

D 

2 

4 

1.0 

— 

N 

5 

2.0 

4.0 

D 

3 

6 

2.0 

4.0 

D 

4 

7 

2.0 

— 

N 

8 

2.5 

17.0 

D 

5 

9 

2.5 

17.0 

D 

10 

3.0 

17.0 

D 

6 

11 

3.0 

17.0 

D 

12 

3.0 

17.0 

D 

7 

13 

3.0 

17.0 

D 

Note:  The  stratification  generated  by  CSETT  will  be  used  in  this  solution. 


* 


1 

i 
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Stratum  No. 

P 

o 

psf 

e 

o 

apexc 

psf 

p 

Post 

psf 

e 

post 

APLoad 

psf 

APFinal 

psf 

ahult 

ft 

2 

4519 

1.267 

4440 

79 

1.460 

1802 

1881 

0.185 

3 

4676 

1.258 

4440 

236 

1.401 

1804 

2040 

0.129 

4 

— 

— 

— 

— 

— 

— 

— 

— 

5 

4865 

1.400 

4440 

425 

1.476 

1811 

2236 

0.042 

6 

4970 

1.398 

4440 

530 

1.468 

1818 

2348 

0.038 

7 

— 

— 

— 

— 

— 

— 

— 

— 

8 

5203 

1.600 

4439 

764 

1.633 

1846 

2610 

0.020 

9 

S334 

1.598 

4439 

895 

1.629 

1903 

2798 

0.019 

10 

5478 

1.510 

4438 

1040 

1.535 

1903 

2943 

0.019 

11 

563* 

1.508 

4437 

1199 

1.532 

1946 

3145 

0.017 

12 

5794 

1.300 

4435 

1359 

1.318 

1996 

3355 

0.014 

13 

5951 

1.298 

4433 

1518 

1.315 

2052 

3570 

0.013 

Note:  P^  =  pre-excavation  overburden  (see  computations) 

eQ  =  pre-excavation  void-ratio  (see  computations) 

AP  =  stress  relief  due  to  excavation  (from  CSETT) 
lXO 

P„  =  post-excavation  stress  =  P  -  P_,,„ 

Post  o  EXC 

epost  =  Popt-excavat:*-on  ratio  (see  computations) 

AP,  ,  =  induced  stress  due  to  loads  (from  CSETT) 

Load 

AP„,  ,  *  stress  after  excavation  and  load  application  =  P„  _  +  P,  , 
Final  Post  Load 

AHyLT  *  total  settlement  under  applied  loads  (see  computations) 


Calculations  for  Table  A14 


Calculation  of  pre-excavation  overburden  (P  )  and  void  ratio  (e  ): 

o  o 


Stratum  2 

Pq  =  (78  -  41)  120  +  ^  (52-5)  =  4519  psf 


Interpolate 


between  appropriate  points  on  input  curve  for 


e 

o 


e 

Logio  : 

1.30 

3.6021 

1.30-x 

3.6550 

1.15 

3.8451 

1.30  -  (1.30-X)  _  3.6021  -  3 
1.30  -  1.15  3.6021  -  3 

X  =  0.0327 


e  =  1.30-X  =  1.267 
o 


Stratum  3 


Pq  =  4519  +  3.0  (52.5)  =  4676  psf 


e 

Logic 

1.30 

3.6021 

1.30-X 

3.6699 

1.15 

3.8451 

1.30  -  (1.30-X)  =  3.6021  -  3 
1.30  -  1.15  3.6021  -  3 

X  =  0.0418 


e  =  1.30-  X  =  1.258 
o 

Stratum  5 


Pq  =  4676  +  1.5  (52.5)  +  1.0  (57.5)  +  1.0  (52.5)  =  4865  psf 
Since  Pq  <  P0INp(JT  ,  compute  eQ  using  CR 

eo  =  e0INPUT  “  CR  Log10  P  =  1,40  “  0,072  log10  4917.5 

OINPUT 

Stratum  6 


1.400 


Po  =  4865  +  2.0  (52.5)  =  4970  psf 

Po  >  P0 INPUT  ’  30  Compute  e0  usi"8  Cc 


6o  "  ®0INPUT  Cc  log10  P 


OINPUT 


1.40  -  0.52  log 


4970 
10  4917.5 


1.398 


Stratum  8 

Pq  -  4970  +  1.0  (52.5)  +  2.0  (57.5)  +  1.25  (52.5)  =  5203  psf 


P  <  P 
o  OINPUT 


,  so  eQ  found  on  recorapression  slope 


6550 

8451 


6699 

8451 
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1.600 


*o  ■  u6  -  °-040  Io810  MfLs 


Stratum  9 


P  =  5203  +  2.5  (52.5)  =  5334  psf 
o 


Po  "  P0INPUT,  30  USe  CC 


1.6  -  0.41  log 


10  5268.75 


1.598 


Stratum  10 

P  =  5334  +  2.75  (52.5)  =  5478  psf 
o 

Po  "  P0INPUT’  USe  CR 
eQ  =  1.51  -  0.035  log  =  1.51 


Stratum  11 


P  =  5478  +  3.0  (52.5)  =  5636 
o 


e  =  1.51  -  0.35  log 
o 


10  5557 


1.508 


Stratum  12 


P  =  5636  +  3.0  (52.5)  =  5794 
o 

S794 

eQ  =  1.30  -  0.028  log1Q  m  1  • 300 


Stratum  13 


P  *  5794  +  3.0  (52.5)  =  5951  psf 
o 

S951 

eo  -  1.30  -  0.30  log^fgn  '  1  *298 


)  and  ultimate  settlement 


Calculation  of  post-excavation  void  ratio  (e post 


Stratum  2 


.  =  e  +  Ae,,  .  ,  =  e  -  C  log.  0  =  1 . 267  -  0.  1 1  log 

Post  o  Rebound  o  r  10  rQ  1 


79 

4519 


=  1.460 


H  _  ,  PFinal  3.0  „  ,,  1881  _  n 

AHULT  =  TTZ -  Cr  log10  ~T  77  =  1  +  1.460  0,11  Iog10  79 

post  post 


Stratum  3 


236 


'Post  -  K258  -  °-n  Io810  ffffi  “  ‘-401 


‘Vi-nr®  »•“  °-129 


Stratum  5 

425 

ep0st  =  1,400  “  °‘°72  1Og10  4865 


1.476 


2.0 


2236 


AHULT  =  1  +  1.476  0,072  1Og10  425 


0.042 


Stratum  6 

530 

ePost  =  1,398  "  °’072  logl0  4970 


1.468 


%LT  '  1  .2K468  °'°72  1O*10  T55  '  °'°38 


Stratum  8 

764 

ePost  =  1,600  ”  0,040  logl0  5203 


1.633 


AHULT 


2.5 

1  +  1.633 


0.040  log1Q 


2610 

764 


0.020 


Stratum  9 


Post 


895 


1.598  -  0.040  log1Q  =  1.627 


4IW  ■  T-ri7629  °-040  Io®io  tH  -  °-019 


.  185 
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Stratum  10 


ePost 

=  1 

.51 

-  0. 

AH 

3. 

0 

ULT 

1 

+  1 

.  535 

Stratum 

11 

6Post 

=  1 

.508 

-  0 

ahult 

=  — 

3. 

0 

1 

+  1 

.532 

Stratum 

12 

ePost 

=  1 

.300 

-  0 

AH 

3. 

0 

ULT 

T 

+  1 

.318 

Stratum 

13 

ePost 

=  1 

.298 

-  0, 

AH 

3J 

0 

ULT 

T 

+  1 

.315 

1040 
!10  5478 


=  1.535 


2943 
’10  1040 


0.019 


1199 
!10  5636 


=  1.532 


0.035  log  3145 
10  1199 


1357 


3355 


1518 
r10  5951 


=  1.315 


3570 


— i.315  °-°28  lo8!0  iflf  =  °-013 
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Table  A15 

Time-Rate  Calculations  for  Strata  2  and  3 


Load  No . 


T1M1 

yr 

0.0 


2.00 


T1H2 

yr  _ 

T 

c 

TIME 

yr 

Relative 

Time 

yr 

T 

U 

% 

2 

.00 

3.778 

0. 10 

0.10 

0.189 

1.6 

0.25 

0.25 

0.472 

6.4 

0.50 

0.50 

0.994 

17.0 

1.00 

1.00 

1.889 

41.3 

2.00 

2.00 

3.778 

91.3 

2.25 

2.25 

4.250 

97.3 

2.50 

2.50 

4.772 

99.2 

5.00 

5.00 

9.444 

100.0 

10.00 

10.00 

18.888 

100.0 

2 

.20 

0.378 

0.10 

0.0 

0.0 

0.0 

0.25 

0.0 

0.0 

0.0 

0.50 

0.0 

0.0 

0.0 

1.00 

0.0 

0.0 

0.0 

2.00 

0.0 

0.0 

0.0 

2.25 

0.25 

0.472 

58.1 

2.50 

0.50 

0.944 

87.0 

5.00 

3.00 

5.667 

100.0 

10.00 

8.00 

15.111 

100.0 

Note : 

Cv  =  17ft2/yr 


H  =  3.0  ft 


Table  A16 

Time-Rate  Calculations  for  Strata  5  and  6 


T1M1 

T1M2 

T 

Load  No. 

yr 

yr 

i 

C 

1 

0.00 

2.00 

8.00 

2 

2.00 

2.20 

0.80 

TIME 

yr 

Relative 

Time 

yr 

T 

u 

z 

0.10 

0.10 

0.400 

2.4 

0.25 

0.25 

1.00 

8.7 

0.50 

0.50 

2.00 

20.9 

1.00 

1.00 

4.00 

45.8 

2.00 

2.00 

8.00 

95.8 

2.25 

2.25 

9.00 

99.7 

2.50 

2.50 

10.00 

100.0 

5.00 

5.00 

20.00 

100.0 

10.00 

10.00 

40.00 

100.0 

0.10 

0.0 

0.0 

0.0 

0.25 

0.0 

0.0 

0.0 

0.50 

0.0 

0.0 

0.0 

1.00 

0.0 

0.0 

0.0 

2.00 

0.0 

0.0 

0.0 

2.25 

0.25 

1.00 

78.4 

2.50 

0.50 

2.00 

98.2 

5.00 

3.00 

12.00 

100.0 

10.00 

8.00 

32.00 

100.0 

Note : 

C  =  16.0  ft2/yr 

V  J 

H  =  2.0  ft. 


Table  A17 


Time-Rate  Calculations  for  Strata  8  Through  13 


T1M1 

T1M2 

T 

Load  No. 

yr 

yr 

C 

1 

0.00 

2.00 

0.304 

2 

2.00 

2.20 

0.03 

TIME 

Relative 

Time 

V 

yr 

yr 

T 

% 

0.10 

0.10 

0.015 

0.5 

0.25 

0.25 

0.038 

1.9 

0.50 

0.50 

0.076 

5.3 

1.00 

1.00 

0.152 

14.7 

2.00 

2.00 

0.304 

41.5 

2.25 

2.25 

0.343 

47.5 

2.50 

2.50 

0.381 

52.5 

5.00 

5.00 

0.761 

81.5 

10.00 

10.00 

1.522 

97.2 

0.10 

0.0 

0.0 

0.0 

0.25 

0.0 

0.0 

0.0 

0.50 

0.0 

0.0 

0.0 

1.00 

0.0 

0.0 

0.0 

2.00 

0.0 

0.0 

0.0 

2.25 

0.25 

0.038 

16.9 

2.50 

0.50 

0.076 

27.9 

5.00 

3.00 

0.457 

72.7 

10.00 

8.00 

1.218 

95.8 

Note: 

C  =  11.0  ft2/yr 
v 

H  =  8.5  ft. 
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c.  Calculations  for  Tables  A15,  A16,  and  A17.  Time-rate  of  consolidation 
computations  (time-dependent  loading) 


C  t 

T  =-V 
C  H2 


where 


for 


C  =  coefficient  of  consolidation 
v 

t£  =  length  of  construction  period  (i.e.,  TIM2-TIM1) 
H  *  length  of  drainage  path 


£  7  100 

c  U  =  j1"  (I-)  2  ^4  [EXP(M2Tc)  -lj  EXP(-M2T)J 


where 


T  >  T 


U  =  %  consolidation 
(2m  +  l)r 


x  100 


M 


,m=o •  ly  2 1  •  •  •  i 


For  time-dependent  loading,  it  is  assumed  that  the  settlement  attributed 
to  each  load  is  proportional  to  the  stress  induced  by  that  load,  or 


AH.  AP,  .. 
i  Loadi 


ahult 


AP 


Load 
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Stratum  designated  as  Input;  subscript  Indicates  Internal  stratification. 
Results  of  the  CSETT  solution  are  shown  In  parentheses  for  comparison. 


Example  Problem  A4 

9.  This  example  illustrates  the  use  of  a  three-dimensional  rectangular 
load  option.  The  soil  stratification  is  shown  in  Figure  A13.  A  strip  footing 
and  two  square  footings,  each  with  a  uniform  load,  are  located  as  shown  in 
Figure  A14.  Soil  properties  are  summarized  in  Table  A19. 

10.  The  input  file  is  listed  in  Figure  A15.  Output  is  generated  along 
the  centerline  of  the  strip  footing  and  along  a  diagonal  between  the  strip 
footing  and  the  square  footings  as  shown  in  Figure  A14.  Time-rate  calcula¬ 
tions  were  made  at  times  0.25,  0.50,  1.0,  2.0,  4.0,  and  8.0  years. 

11.  A  hand  solution  at  position  (X,Y)  =  (35.0,  32.5)  is  provided  for  com¬ 
parison  (Tables  A20-A22) .  Results  of  the  hand  solution  and  the  CSETT 
solution  are  compared  in  Table  A23. 
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DISTANCE,  X  (FT) 


Figure  AI4.  Plan  view  of  footing  loads  for  Example  Problem  A4 


Table  A19 

Soil  Properties  for  Example  Problem  A4 


El  of  Top  Eff  Unit 
of  Strata  Weight 
ft  NGVD  jpcf 


_ Consolidation  Data 

C  P 

2V  o 

ft  /vr  eo  nsf  LR 


THREE  FOOTING 


10  TITLE 

20  3-DIMENSIONAL  RECTANGULAR  LOADS  - 

30  3DRE  1  0.0  0  150.0 

40  10. 0  30 . 0  50 . 0  35 . 0 

50  3DRE  2  0.0  0  200.0 

60  10.0  7.5  17.5  15.0 

70  3DRE  3  0.0  0.0  200.0 

SO  45.0  7.5  52.5  15.0 

90  BOUS  7.5 

100  SOIL  1  35.5  N  115. 

110  SOIL  2  84.5  D  52.5  .04  3.0 

120  INDEX  .35  193.75  1.6 

130  SOIL  3  31.5  N  57.5 

140  SOIL  4  31.0  D  52.5  .045  12.0 

150  INDEX  .43  330.0  1.3 

160  SOIL  5  78.0  N  57.5 

170  TIMS  .25  .50  1.0  2.0  4.0  8.0 

180  OUTPUT  10.0  32.5  60.0  32.5  25.0 

190  OUTPUT  10.0  35.0  52.5  7.5  10 

200  END 


* 


Figure  A15.  Input  data  file  for  Example  Problem  A4. 


TS  TNPIIT  FRflM  TFRMTNAI  OR  A  FTIF 
FNTFR  T  OR  F 
=F 


FNTFR  HAT  A  FTI  F  NAMF 
=RFOT 


TMF'IIT  fOMPI  FTF.  Hfi  Yfil  I  WANT  T  MPI  IT  RATA 
F CHOPR T IMTFn  TO  Yfil  IR  TFRMTNAI  ,  A  FTI  F - 
ROTH,  HR  MFTTHFR  '  ( f  I JTFR  T ,  F.  Fi,  HR  IM ) 

-T 


PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/83  TIME:  9:55:29 


I.  INPUT  DATA 


1.  TITLE  -  3-DIMENSIONAL  RECTANGULAR  LOADS  -  THREE  FOOTINGS 


2.  B0U3SINESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES  . 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  7.50  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 


I 


A64 


LOAD  NUMBER  1  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YR 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 

PRESSURE  =  150.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  10.00,  30.00 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  50.00,  35.00 


LOAD  NUMBER  2  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YR 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 

PRESSURE  =  200.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  10.00,  7.50 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  17.50,  15.00 


LOAD  NUMBER  3  BEGINNING  TIME  OF  LOAD  APPLICATION  =  0.  YR 

ENDING  TIME  OF  APPLICATION  =  0.  YRS. 

PRESSURE  =  200.000  PSF 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  LOWER-LEFT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  45.00,  7.50 

COORDINATE  PAIR  (PLAN  VEIW)  DEFINING  UPPER-RIGHT 
CORNER  OF  LOAD  REGION  -  X  ,  Y  =  52.50,  15.00 


6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 
NONE 


7.  EXCAVATION  DATA 
NONE 
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8.  SOIL  DATA 


STRATA 

EL. 

OF  TOP 

DRAINAGE 

EFF  UNIT  RECOMPR. 

COEF . OF 

POISSON 

NO. 

OF 

STRATUM 

COND I T I ON 

WEIGHT  INDEX 

CONSOL. 

RATIO 

(FEET  NGVD ) 

( PCF  > 

(SQFT/YR) 

1 

85 . 50 

N 

115.00 

2 

84.50 

D 

52.50  0.04000 

8 . 00000 

0. 32000 

o 

81.50 

N 

57 . 50 

4 

81.00 

D 

52.50  0.04500  1 

2 . 00000 

0. 32000 

9.  STRE 

SS- 

■STRAIN  DATA 

STRATUM 

NO.  1 

I NCGMPRESS IDLE  STRATUM 


STRATUM  NO.  2 


COMPRESSION  INDEX=  0.35000 
RECOMPRESSION  IND£X=  0.04000 
I NS ITU  VOID  RAT 10=  1.60000 
I NS ITU  OVERBURDEN=  193.75  PSF 


STRATUM  NO.  3 
INCOMPRESSIBLE  STRATUM 

STRATUM  NO.  4 


COMPRESSION  INBEX=  0.48000 
RECOMPRESSION  INDEX=  0.04500 
I NS ITU  VOID  RAT 10=  1.30000 
I NS ITU  OVERBURDEN^  380.00  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


TIME  RATE  OF  CONSOLIDATION  CALCULATIONS  WILL  BE  MADE 
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AT  TIMES  ( YRS ) 


0.25 
0.50 
1 . 00 
2 . 00 
4 . 00 
8.00 


OUTPUT  CONTROL  DATA 


PLANE  1  XA= 

10. 0000 

FT 

YA= 

32. 5000 

FT 

XB= 

60 . 0000 

FT 

YB= 

32. 5000 

FT 

DELX= 

25 . 0000 

FT 

PLANE  2  XA= 

1 0 . 0000 

FT 

YA= 

35 . 0000 

FT 

XB= 

52 . 5000 

FT 

YB= 

7 . 5000 

FT 

DELX  = 

1 0 . 0000 

FT 

PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/33  TIME:  9:56: 33 


II.  OUTPUT  SUMMARY. 


1.  TITLE-  3-DIMENSIONAL  RECTANGULAR  LOADS  -  THREE  FOOTINGS 


POSITION:  X=  10.0  Y=  32.5 

*#*#*•*■&  ********  *****■***•»*■«■*** 


IJMMARY 

OF  ULTIMATE 

SETTLEMENTS. 

STRATA 

MID-DEPTH 

IN— SITU 

DELTA 

ULTIMATE 

NO. 

OF  STRATA 

OVERBURDEN 

SIGMA 

SETTLEMENT 

(FEET) 

(LB /SO  FT) 

(LB/SQ  FT) 

(FEET) 

1 

0 . 50 

57.  50 

74 .  SO 

0. 

2.50 

193.75 

6 1 . 40 

0. 048 

4.25 

♦“•O/.,  B  oo 

46.30 

0. 

4 

6. 00 

3S0. 00 

36 .  00 

0.025 

3.  TIME-SETTLEMENT  SUMMARY. 


STRATA 

NO 

(SETTLEMENT 
IJLT  0.25 

( YRS. ) 

IN  FEET 
0.  50 
(YRS. ) 

AT  SPECIFIED 
1.00  2.00 
(YRS.)  (YRS.) 

TIMES) 

4 . 00 

(YRS. ) 

8 .  00 
( YRS . ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.048 

0.044 

0. 048 

0.043 

0.048 

0.048 

0.  048 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.025 

0.024 

U  •  Oii'5 

0.  025 

0.  025 

0. 025 

0.  025 

TOTALS: 

0.073 

0.063 

0.073 

0.073 

0.073 

0.073 

0.073 

) 


Abbreviated  Output 


PROGRAM  GGFTT  -  VFRTTGAI 
PATF:  1 1/1  7/p? 


PTRFPP  TNiriljr  TT  ON  A  Kin  SFTTI  FMFKJT  PROGRAM 

TTMF:  1f>:  o:  7 


T  T  .  01  IT  PI  IT  P.IJMMARY . 


1  .  TTTI  F-  TMFNSTOIMAI.  RFGTAIJGI II  AR  I  ft  A  OP  -  THRFF  FOi'iT  T  NGS 


PI  IMMARY  OF  TTMF  SFTTI  FMFIvIT  DATA, 


PI  AI\IF  OF  IMTFRFST  1 


KA=  PS.  ft  YA  = 
yp-.  40.0  YR= 

ni-i  V-  io.o 


«•  ,  .  .  I 

? .  =1 


TTMF 

(YR> 


y- 

Y— 


s.  <t 


III  T. 

0.1? 

0 .  ?s 

0.  \  7 

o.  SO 

0 .  1 

1 . 00 

0 .  1 :? 

? .  Of) 

0.  1  ? 

4 .  Of) 

ft .  1  ? 

p.  Oft 

0.  1  ■ 
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POSITION:  X=  35.0  V=  32.5 

***************************** 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN-SITU 
OVERBURDEN 
(LB/SQ  FT) 

DELTA 
SIGMA 
(LB/SQ  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

0. 50 

57 . 50 

149.50 

0. 

Am 

2.  50 

193. 75 

122.70 

0.086 

4 . 25 

286. 88 

92 . 40 

0. 

4 

6 . 00 

380 . 00 

71 . 40 

0.047 

3.  TIME-SETTLEMENT  SUMMARY. 


STRATA  IJLT 
NO 


0 .  0 .  <: 
0.086  0.078  C 

o.  o.  <: 


4 

TOTALS: 


iCIFIED 

TIMES) 

2.  00 

4 . 00 

8.  00 

( YRS . ) 

(YRS. > 

(YRS. ) 

0. 

0 . 086 

0. 

0.  086 

0. 

0 . 086 

0. 

0. 

0. 

0.047 

0 .047 

0.047 

0.  133 

0.  133 

0.  133 

POSITION:  X=  60.0  Y=  32.5 

ft**************************** 


1.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 

STRATA  MID-DEPTH  IN-SITU  DELTA  ULTIMATE 
NO.  OF  STRATA  OVERBURDEN  SIGMA  SETTLEMENT 
(FEET)  (LB/SQ  FT)  (LB/SO  FT)  (FEET) 


0.50 

2.50 

4.25 

6.00 


57.  50 
193.75 
286. 88 
380 . 00 


0. 

0. 

0. 

0 . 00 1 


TIME-SETTLEMENT  SUMMARY 


STRATA 

NO 


1 


4 

TOTALS 

DO  YOU 
ENTER 
=  N 


( SETTLEMENT 

IN  FEET  AT  SPECIFIED 

TIMES) 

LILT 

0.25 

0.50 

1 . 00 

2.00 

4.00 

8.00 

—  —  — 

<  YRS.  ) 

(YRS. > 

(YRS.  ) 

(YRS. ) 

(YRS. ) 

(YRS. ) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 . 00  i 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0 . 00 1 

0 . 00 1 

0. 001 

0 . 00 1 

0.  001 

0.  001 

0 . 00 1 

WANT  A 

PLOT  OF  S 

TRESS 

BULBS  ? 

Y  -'  OR  -  N  '  . 


HAND  VERIFICATION  POINT  SETTLEMENT  AT  (X,Y)  =  (35,  32.5) 


i  ... 


. 1  ’IW 


Table  A20 

Internally  Generated  Stratification  for  Example  Problem  A4* 


Input 

Stratum  No. 

Internally 
Generated 
Strata  No. 

Thickness  for 
Compression 
Calculations 
ft 

Eff  Thickness 
for  Consolidation 
Calculations 
ft 

Drainage  Condition 

1 

1 

— 

— 

N 

2 

2 

3.0 

3.0 

D 

3 

3 

— 

— 

N 

4 

4 

3.0 

3.0 

D 

5 

5 

— 

— 

N 

*  Note:  Since  all  strata  thicknesses  are  less  than  or  equal  to  3  ft,  the 
internally  generated  stratification  is  the  same  as  the  input 
stratification. 
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Table  A21 


Summary  of  Stresses  and  Ultimate  Settlements 


Stratum  No. 

P 

0 

psf 

e 

o 

APT  j  * 
Loads* 

psf 

AP„.  , 
Final 

psf 

ahult 

ft 

1 

— 

— 

— 

— 

— 

2 

193.75 

1.60 

122.73 

316.48 

0.086 

3 

— 

— 

— 

— 

— 

4 

380.0 

1.30 

71.42 

451.42 

0.047 

5 

__ 

— 

— 

— 

*  From  CSETT 

Note:  Calculation  of  in  situ  void  ratio  (e  )  and  overburden  (P  ) 

o  o 

Stratum  2 

P  =  1.0  x  115  +  1.5  x  52.5  =  193.75  psf 
o 

Since  Pq  =  P0INpur  e0  *  e0INPUT  1,60 
Stratum  4 


P 

o 


1.0  x  115.0  +  3.0  x  52.5  +  0.50  x  57.5  +  1.5  x  52.5  =  380.0  psf 


Since  Pq  p0inPUT’  ®o  =  1,30 
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Table  A22 

Time-Rate  of  Consolidation  and  Ultimate  Settlement  Calculations 


Stratum  No. 

hult 

ft 

H 

ft 

t 

yr 

T 

V 

U 

% 

Settlement 

ft 

2 

0.086 

1.50 

18.0 

0.25 

0.889 

91.0 

0.078 

0.50 

1.78 

99.0 

0.085 

1.00 

3.56 

100.0 

0.086 

2.00 

7.11 

100.0 

0.086 

4.00 

14.22 

100.0 

0.086 

8.00 

28.44 

100.0 

0.086 

4 

0.047 

1.50 

12.0 

0.25 

1.33 

97.0 

0.045 

0.50 

2.67 

100.0 

0.047 

1.00 

5.33 

100.0 

0.047 

2.00 

10.67 

100.0 

0.047 

4.00 

21.33 

100.0 

0.047 

8.00 

42.67 

100.0 

0.047 

Note:  Calculation  of  ultimate  settlements. 
Stratum  2 


ahult 


-HCc...  log  !f  =  1^0.35) 

1  +  e  g10  P  1+1.60 
o  o 


log 


316.48 
10  193.75 


=  0.086  ft 


Stratum  4 


3.0(0.48)  451.42 

^JLT  1  +  1.30  iog10  380.0 


0.047  ft 


Table  A23 

Summary  of  Time-Rate  Calculations* 


Stratum  No. 

ahult 

ft 

Settlement 

in  Feet 

at  Specified  Time 

0.25 

yr 

0.50 

yr 

1.00 

yr 

2.00 

yr 

4.00 

yr 

8.00 

yr 

2 

(0.086) 

(0.078) 

(0.085) 

(0.086) 

(0.086) 

(0.086) 

(0.086) 

0.086 

0.078 

0.085 

0.086 

0.086 

0.086 

0.086 

4 

(0.047) 

(0.045) 

(0.047) 

(0.047) 

(0.047) 

(0.047) 

(0.047) 

0.047 

0.045 

0.047 

0.047 

0.047 

0.047 

0.047 

Total 

(0.133) 

(0. 123) 

(0.132) 

(0.133) 

(0.133) 

(0.133) 

(0.133) 

0. 133 

0. 123 

0.132 

0. 133 

0.133 

0.133 

0.133 

*  Results  from  the  CSETT  solution  are  shown  in  parentheses  for  comparison. 
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12.  This  example  provides  an  illustration  of  three-dimensional  irregular 
loads.  The  loading  consists  of  a  fill  and  a  storage  tank  as  shown  in  Fig¬ 
ure  A16.  The  soil  stratification  is  shown  in  Figure  A17,  and  the  soil  prop¬ 
erties  are  summarized  in  Table  A24. 

13.  The  input  data  file  is  listed  in  Figure  A18.  Output  is  generated  at 
20-ft  intervals  along  the  plane  of  interest  shown  in  Figure  A16.  Time-rate 
calculations  are  made  at  times  of  0.10,  0.20,  0.40,  0.80,  1.50,  3.0,  6.0,  and 
12  years. 

14.  A  hand  solution  at  position  (x,y)  =  (30,30)  is  provided  (Tables  A25- 
A27).  Results  of  the  hand  solution  and  the  CSETT  solution  are  compared  in 
Table  A28. 


DISTANCE  X  FT 

Figure  A16.  Plan  and  elevation  views  of  loads  for 
Example  Problem  A5 
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EL  84.5  ||  Ij 

^  T«  11 5.0  PCF 

* 

T  *52.5  PCF 

EL  81.5 

EL  81.0  pj] 

7' =57.5  PCF 

EL  78.0 

T-  525  PCF 

Figure  A17.  Soil  stratification  for  Example  Problem  A5 


Table  A24 

Soil  Properties  for  Example  Problem  A5 


Stratum 

No. 


El  of  Top 
of  Stratum 
ft  NGVD 


Eff  Unit 
Weight 
f 


2V 

ff  hr 


Consolidation  Data 
P 


193.75 


380.00 


SOLUTION  BY  CSETT 
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i 

i  4 


LIST  IRR 


10  TITLE 

20  3-DIM.  IRREGULAR  LOADS  -  FILL  AND  STORAGE  TANK 

30  3DIR  1  0  5  0.0  0.0  150.0  (3-D  IRREGULAR  LOAD  DATA  NL  NS  N  T1M1  T1M2  PRES) 

40  10.0  7.5  10.0  52.5  70.0  52.5  70.0  7.5  10.0  7.5  (X.  Y  ,  i=l,N) 

50  3DIR  2  1  5  0.0  0.0  300.0  (3-D  IRREGULAR  LOAD  DATA  NL  NS  N  T1M1  T1M2  PRES) 

60  12.5  10.0  12.5  50.0  67.5  50.0  67.5  10.0  12.5  10.0  (X.  Y  ,  i=l,N) 

70  3DIR  3  2  5  0.0  0.0  450.0  (3-D  IRREGULAR  LOAD  DATA  NL  NS  N  T1M1  T1M2  PRES) 

80  15.0  12.5  15.0  47.5  65.0  47.5  65.0  12.5  15.0  12.5  (X.  Y  ,  i=l,N) 

90  3DIR  4  0  17  0.0  0.0  1250.0  (3-D  IRREGULAR  LOAD  DATA  NL  NS  N  T1M1  T1M2  PRES) 

100  25.0  30.0  26.0  35.0  29.0  40.0  34.0  44.0 

110  40.0  45.0  46.0  44.0  51.0  40.0  54.0  35.0  X.  Y.  ,  i=l,N 

120  55.0  30.0  54.0  25.0  51.0  20.0  46.0  16.0  1  1 

130  40.0  15.0  34.0  16.0  29.0  20.0  26.0  25.0 

140  25.0  30.0 

150  BOUS  7.5  (BOUSSINESQ  SOLUTION  DMAX) 

160  SOIL  1  85.5  N  115.  (SOIL  DATA  N  EL  NDC  GM) 

170  SOIL  2  84.5  D  52.5  .04  8.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

180  INDEX  .35  193.75  1.6  (INDEX  PROPERTIES  CC  PO  EO) 

190  SOIL  3  81.5  N  57.5  (SOIL  DATA  NEL  NDC  GM) 

200  SOIL  4  81.0  D  52.5  .045  12.0  (SOIL  DATA  N  EL  NDC  GM  CR  CV) 

210  INDEX  .48  380.0  1.3  (INDEX  PROPERTIES  CC  PO  EO) 

220  SOIL  5  78.0  N  57.5  (SOIL  DATA  N  EL  NDC  GM) 

230  TIMS  .1  .2  .4  .8  1.5  3.0  6.0  12.0  (SPECIFIED  TIME  T1  T2  T3  ,..) 

240  OUTPUT  10  30  70  30  20  (OUTPUT  XI  Y1  X2  Y2  SP) 

250  END  (TERMINATION) 

Figure  A18.  Input  data  file  for  Example  Problem  A5 


Ts  TNPIIT  FROM  TFRMTNAI  mR  A  FTI  F 
FWTFR  T  HR  t- 
=F 


FNTFR  DATA  FT  IF  NAMF 
=  TRR 


tnpiit  r  riMPi  ftf.  nri  ym  i  uaijt  tkjpiit  fata 
FCHOPR T NTFP  Tfi  YOIIR  TFRMTNAI  .  A  FTIP, 
ROTH,  OR  NFTTHFR  T-  t'  FNTFR  T,  F,  R,  OR  |\j ) 


PROGRAM  C SETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/83  TIME:  10:  0: 17 


I.  INPUT  DATA 


1.  TITLE  -  3-DIM.  IRREGULAR  LOADS  -  FILL  AND  STORAGE  TANK 


2.  B0US3INESQ  SOLUTION  WILL  BE  USED  TO  COMPUTE  INDUCED  STRESSES 
THE  MAXIMUM  DEPTH  TO  WHICH  THE  ANALYSIS  WILL  BE  EXTENDED 
IS  7.50  FEET. 


3.  2-DIMENSIONAL  PRESSURE  LOAD  DATA 
NONE 


4.  2-DIMENSIONAL  SOIL  LOAD  DATA 
NONE 


5.  3-DIMENSIONAL  RECTANGULAR  LOAD  DATA 
NONE 
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6.  3-DIMENSIONAL  IRREGULAR  LOAD  DATA 


LOAD  NUMBER  1 


LOAD  NUMBER 


LOAD  NUMBER  3 


BEGINNING 

TIME  OF 

LOAD  af-plicat: 

ENDING  TIME  OF  APPLICATION  =  0. 

PRESSURE  = 

150. 

000  PSF 

POINT  NO. 

X 

Y 

(FT.  ) 

(FT.  ) 

1 

10.  00 

7.50 

2 

10.00 

52.50 

o 

■J 

70.  00 

52.50 

4 

70.00 

7.50 

5 

1 0 . 00 

7 . 50 

BEGINNING 

TIME  OF 

LOAD  APPLICAT 

ENDING  TIME  OF  APPLICATION  =  0. 

PRESSURE  = 

300 . 

000  PSF 

POINT  NO. 

X 

Y 

(FT.  ) 

(FT.  ) 

1 

12.50 

10.00 

2 

12.50 

50.00 

;T| 

67.50 

50.00 

4 

67.50 

10.00 

5 

12.50 

10.00 

BEGINNING  TIME  OF  LOAD  APPLICATION  =  0. 
ENDING  TIME  OF  APPLICATION  =  0.  YRS, 
PRESSURE  =  450.000  PSF 


POINT  NO. 


X 

(FT.  ) 

15.00 
15.  00 
65.  00 
65 . 00 
15.00 


Y 

(FT.  ) 

12.50 

47.50 

47.50 

12.50 

12.50 


LOAD  NUMBER  4  BEGINNING  TIME  OF  LOAD  APPLICATION 

ENDING  TIME  OF  APPLICATION  =  0. 
PRESSURE  =  1250.000  PSF 


POINT  NO.  X  Y 


(FT.  ) 

(FT.  ) 

1 

25.  00 

30.00 

2 

26.00 

35.00 

3 

29.  00 

40.00 

4 

34 . 00 

44 . 00 

5 

40.00 

45.00 

6 

46.00 

44.00 

7 

51.00 

40.00 

8 

54 . 00 

35.  00 

9 

55.00 

30 . 00 

10 

54.00 

25.00 

11 

5 1 . 00 

20.  00 

12 

46.  00 

16.00 

13 

40.  00 

1 5 . 00 

14 

34.00 

16.00 

15 

29.00 

20.  00 

16 

26.00 

25.00 

17 

25.  00 

30 . 00 

7.  EXCAVATION  DATA 

NONE 

8.  SOIL 

DATA 

STRATA 

EL.  OF  TOP 

DRAINAGE 

EFF  UNIT  RECOMPR 

.  COEF.OF 

NO. 

OF  STRATUM 

CONDITION 

WEIGHT  INDEX 

CONSOL. 

(FEET  NGVD) 

( PCF ) 

(SQFT/YR) 

1 

85.50 

N 

115.00 

2 

84.50 

D 

52.50  0.04000 

8.00000 

3 

81.50 

N 

57.50 

4 

81.00 

D 

52.50  0.04500 

12.00000 

9.  STRESS-STRAIN  DATA 

STRATUM  NO.  1 


INCOMPRESSIBLE  STRATUM 


0.  YRS. 
YRS. 


POISSON'S 

RATIO 

0.32000 

0.32000 
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STRATUM  NO 


COMPRESSION  INDEX=  0.35000 
RECOMPRESSION  INDEX=  0.04000 
I  NS  ITIJ  VOID  RAT  10=  1.60000 
I NS ITU  OVERBURDEN=  193.75  PSF 


STRATUM  NO.  3 


INCOMPRESSIBLE  STRATUM 


STRATUM  NO.  4 


COMPRESSION  I NDE X  =  0.43000 
RECOMPRESSION  INDEX=  0.04500 
I NS ITU  VOID  RAT 10=  1.30000 
I NS ITU  OVERBURDEN=  330.00  PSF 


10.  TIME  SEQUENCE  FOR  CONSOLIDATION  CALCULATIONS 


TIME  RATE  OF  CONSOLIDATION  CALCULATIONS  WILL  BE  MADE 
AT  TIMES  (YRS): 

0.  10 
0.  20 
0.40 
0 . 80 
1.50 
3.00 
6.00 
12.00 


11.  OUTPUT  CONTROL  DATA 


PLANE  1 

XA= 

1 0 . 0000 

FT. 

YA= 

30. 0000 

FT. 

XB= 

70.0000 

FT. 

YB= 

30. 0000 

FT. 

DELX  = 

20.0000 

FT. 

DO  YOU  WANT  TO  EDIT  DATA?  ENTER  'V'  OR  'N'. 
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PROGRAM  CSETT  -  VERTICAL  STRESS  INDUCTION  AND  SETTLEMENT  PROGRAM 
DATE:  7/12/83  TIME:  10:  2:25 


1 1 .  OUTPUT  SUMMARY . 


1.  TITLE-  3-DIM.  IRREGULAR  LOADS  -  FILL  AND  STORAGE  TANK 


POSITION:  X=  10.0  Y=  30.0 

****#**«-**»•••»*»***#******** 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


•TRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN-SITU 
OVERBURDEN 
(LB/S 0  FT) 

DELTA 
SIGMA 
(LB /SO  FT) 

ULTIMATE 
SETTLEMENT 
( FEET ) 

1 

0.50 

57.50 

75.30 

0. 

2 

2.50 

193.75 

92.  10 

0.068 

o 

4.25 

286 . 88 

115.80 

0. 

4 

6.00 

380. 00 

137.50 

0. 084 

3.  TIME-SETTLEMENT  SUMMARY. 


STRATA 

NO 

(SETTLEMENT 
IJLT  0.10 

( YRS -  1 

IN  FEET  AT  SPECIFIED 
0.20  0.40  0.80 

<YRS. )  (YRS. )  (YRS. ) 

TIMES) 

1 .50 
(YRS. > 

3 . 00 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

*-> 

0 . 068 

0.045 

0.059 

0.067 

0.  068 

0 . 068 

0 . 068 

3 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.084 

0 . 066 

0.079 

0 . 084 

0.084 

0. 084 

0.  084 

TOTALS: 

0.  152 

0.111 

0.  138 

0.  151 

0.  152 

0.  152 

0.  152 
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STRATA 

NO 

<  SETTLEMENT 
6.00  12.00 
(YRS. )  ( YRS. ) 

1 

0. 

0. 

O 

■C. 

0 . 068 

0 . 068 

3 

0. 

0. 

4 

0.084 

0 . 084 

TOTALS: 

0.  152 

0.  152 

IN  FEET  AT  SPECIFIED  TIMES) 


POSITION:  X=  30.0  Y=  30.0 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


:TRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN-SI TU 
OVERBURDEN 
(LB/ SO  FT) 

DELTA 
SIGMA 
(LB/ SO  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

0.50 

57.50 

1 699 . 60 

0. 

2.50 

193.75 

1661.00 

0.  396 

4.25 

'-I o  /  oo 

1 569 . 90 

0. 

4 

6 . 00 

380 . 00 

1 458 . 60 

0.429 

TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

ULT 

0.  10 
(YRS.  ) 

0.  20 
(YRS. ) 

0.  40 
(YRS. ) 

0 .  80 
(YRS. ) 

1 . 50 
(YRS. ) 

3.00 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.396 

0.  263 

0.341 

0.387 

0 . 396 

0 .  396 

0 . 396 

3 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.429 

0.335 

0.404 

0.  427 

0.429 

0.429 

0.429 

TOTALS: 

0.825 

0.598 

0.  745 

0 .814 

0.825 

0.  825 

0.825 

STRATA 

NO 


(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 

6 . 00  1 2 . 00 
(YRS.)  (YRS. ) 


1  0.  0. 


2 

0 . 396 

0.  396 

0. 

0. 

4 

0.429 

0.429 

TOTALS: 

0.825 

0.825 

POSITION:  X=  50.0  Y=  30.0 

***•*•&«•*****■#•**■***#*#•********* 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


TRATA 

NO. 

MID- DEPTH 
OF  STRATA 
(FEET) 

IN-SITU 
OVERBURDEN 
(LB /SO  FT) 

DELTA 
SIGMA 
(LB/SO  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

0.50 

57.50 

1699.60 

0. 

2.50 

193.75 

1661.00 

0.396 

3 

4.25 

286.88 

1569.90 

0. 

4 

6.00 

380.00 

1458.60 

0.  429 

3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

'ULT 

0.  10 
(YRS.  > 

0.20 
(YRS. ) 

0.40 
(YRS. ) 

0.80 
(YRS. ) 

1.50 
(YRS. ) 

3.00 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

2 

0.396 

0.  263 

0.341 

0.387 

0.396 

0.  396 

0.  396 

3 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.429 

0.335 

0.  404 

0.427 

0.429 

0.  429 

0.429 

TOTALS: 

0.825 

0.598 

0.745 

0.814 

0.825 

0.825 

0.825 

► 

► 


STRATA 

NO 


(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 
6.00  12.00 
(YRS.)  ( YRS. ) 


1 

0. 

0. 

2 

0.  396 

0.396 

0. 

0. 

4 

0.429 

0.429 

TOTALS: 

0.825 

0.825 

POSITION:  X=  70.0  Y=  30.0 

■&******■****•*#**■*•#■****■**■  ***■*■»* 


2.  SUMMARY  OF  ULTIMATE  SETTLEMENTS. 


STRATA 

NO. 

MID-DEPTH 
OF  STRATA 
(FEET) 

IN-SITU 
OVERBURDEN 
(LB /SO.  FT) 

DELTA 
SIGMA 
(LB/SG!  FT) 

ULTIMATE 

SETTLEMENT 

(FEET) 

1 

0.50 

57. 50 

75.30 

0. 

2 

2.50 

193. 75 

92.  10 

0.068 

3 

4.25 

286.88 

115.80 

0. 

4 

6.  00 

380. 00 

137.50 

0.084 

3.  TIME-SETTLEMENT  SUMMARY. 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 


STRATA 

NO 

ULT 

0.  10 
(YRS.  ) 

0.  20 
(YRS. ) 

0.40 
(YRS. ) 

0.80 
(YRS. ) 

1.50 
(YRS. ) 

3.00 
(YRS. ) 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4- 

0.  068 

0.045 

0 . 059 

0.067 

0.068 

0.068 

0.068 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

4 

0.  084 

0 . 066 

0.079 

0.084 

0.  034 

0.084 

0.084 

TOTALS: 

0.  152 

0.111 

0.  138 

0.  151 

0.  152 

0.  152 

0.  152 

A88 


r 


STRATA 

(SETTLEMENT  IN  FEET  AT  SPECIFIED  TIMES) 

6.00  12.00  1 

NO 

(YRS. ) 

(YRS.)  f 

1 

0. 

o.  i 

0.068 

0.068  ' 

0. 

0. 

4 

0.084 

0.084 

TOTALS: 

0.  152 

0.  152 

* ; 

N 


M 


f 


HAND  VERIFICATION  AT  SETTLEMENT  UNDER  (X,Y)  =  (30,30) 
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Table  A25 

Internally  Generated  Stratification  for  Example  Problem  A5 


Input 

Stratum  No. 

Internally 
Generated 
Stratum  No. 

Thickness  for 
Compression 
Calculations 
ft 

Eff  Thickness 
for  Consolidation 
Calculations 
ft 

Drainage  Condition 

1 

1 

— 

— 

N 

2 

2 

3.0 

3.0 

D 

3 

3 

— 

— 

N 

4 

4 

3.0 

3.0 

D 

5 

5 

_ 

_ 

N 

Note:  Since  all  strata  thicknesses  are  less  than  or  equal  to  three  (3)  ft, 
the  internally  generated  stratification  is  the  same  as  the  input 
stratification. 
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Table  A26 

Summary  of  Stresses  and  Ultimate  Settlements 


for  Example  Problem  A5 


Stratum  No. 

P 

o 

psf 

e 

o 

APt  ,  * 
Loads* 

psf 

PFinal 

psf 

ahult  j 

ft  ; 

| 

l  1 

i 

2 

o 

193.75 

1.60 

1661.00 

1845.75 

0.396 

1. 

4 

380.0 

1.30 

1458.60 

1838.60 

* 

0.429 

5 

— 

— 

— 

— 

y 

Note:  Calculation  of  in  situ  void  ratio  (e  )  and  overburden  (P  ). 

o  o 

Stratum  2 

P  =  1.0  x  115  +  1.5  x  52.5  =  193.75  psf 
o  r 

Since  Pq  =  P0INpuT.  eD  =  eoiNPUT  1,60 

Stratum  4 

P  =  1.0  x  115.0  +  3.0  x  52.5  +  0.5  x  57.5  +  1.5  x  52.5  =  380.0  psf 
o  v 

Since  ?q  -  P0INpuT  ,  eQ  =  eoiNPUT  1  * 30 


*  From  CSETT. 
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Table  A27 

Time-Rate  of  Consolidation  Calculations  for  Example  Problem  A5 


Stratum 

ahult 

H 

C 

t 

U 

Settlement 

No. 

ft 

ft 

ftZ/yr 

V 

% 

ft 

2 

0.396 

1.50 

8.0 

0.10 

0.356 

66.3 

0.262 

0.20 

0.711 

86.0 

0.340 

0.40 

1.422 

97.6 

0.386 

0.80 

2.844 

99.9 

0.396 

1.50 

5.333 

100.0 

0.396 

3.00 

10.67 

100.0 

0.396 

6.00 

21.33 

10O.0 

0.396 

12.00 

42.67 

100.0 

0.396 

4 

0.429 

1.50 

12.0 

0. 10 

0.533 

78.3 

0.336 

0.20 

1.067 

94.2 

0.404 

0.40 

2.133 

99.6 

0.427 

0.80 

4.267 

100.0 

0.429 

1.50 

8.00 

100.0 

0.429 

3.00 

16.00 

100.0 

0.429 

6.00 

32.00 

100.0 

0.429 

12.00 

64.00 

100.0 

0.429 

Note:  Calculation  of  ultimate  settlements. 


Stratum  2 


HC  Pf 

AHULT  =  I  +  e  log10  P~ 
o  o 


3.0(0.35)  .  1854.75 

1  +  1.60  810  193.75 


=  0.396  ft 


Stratum  4 


Table  A28 

Summary  of  Time-Settlement  Data  for  Example  Problem  A5* 


ahult 

Settlement,  FT,  at  Specified  Time 

0. 10 

0.20  0.40  0.80  1.50 

3.00 

Stratum  No. 

ft 

yr 

yr 

yr 

yr 

yr 

yr 

O 

(0.396) 

(0.263) 

(0.341) 

(0.387) 

(0.396) 

(0.396) 

(0.396) 

L. 

0.396 

0.262 

0.340 

0.386 

0.396 

0.396 

0.396 

4 

(0.429) 

(0.335) 

(0.404) 

(0.427) 

(0.429) 

(0.429) 

(0.429) 

0.429 

0.336 

0.404 

0.427 

0.429 

0.429 

0.429 

Total 

(0.825) 

(0.598) 

(0.745) 

(0.814) 

(0.825) 

(0.825) 

(0.825) 

0.825 

0.598 

0.744 

0.813 

0.825 

0.825 

0.825 

*  Results  of  the  CSETT  solution  are  shown  in  parentheses  for  comparison. 


APPENDIX  B:  GENERAL  EQUATIONS  FOR  STRESS  INDUCTION* 


INTRODUCTION 


Although  much  literature  exists  on  settlement  caused  by  stresses  induced 
when  a  load  is  placed  on  the  surface  of  a  soil  mass,  few  texts  contain 
expanded  formulas  for  computing  the  induced  vertical  stresses.  The  equation 
found  in  the  classical  texts  on  soils  mechanics  and  foundation  design  are  nor¬ 
mally  those  for  point  loads,  distributed  line  loads,  "nice"  geometric  shapes, 
and,  in  general,  very  limited  loading  conditions.  These  equations  were  cre¬ 
ated  by  integration  of  the  Boussinesq  or  Westergaard  point  load  formulas  over 
the  shapes  to  be  considered.  Applying  the  same  techniques  to  more  general 
elemental  regions,  the  formulas  can  be  expanded  to  do  the  most  general  two-  or 
three-dimensional  problems. 

In  recent  years,  many  digital  programs  have  been  created  using  either 
finite  element  methods  or  finite  difference  testers  to  approximate  induced 
stresses.  No  one  is,  nor  should  be,  concerned  about  the  approximations,  since 
even  "exact"  solutions  are  merely  academics.  However,  the  user  is  burdened  by 
having  to  input  much  of  the  bookkeeping  necessary  to  adequately  describe  the 
loading  condition.  Also,  a  program  burdened  with  the  required  secretarial 
logic  can  do  little  more  than  compute  vertical  stresses.  Thus,  the  ultimate 
goal  of  settlement  analysis  must  be  done  by  subsequent  activities,  if  all  of 
the  aspects  of  consolidation  such  as  ultimate,  time-rate,  graphics,  etc.,  are 
to  be  included. 

The  purpose  of  this  paper  is  to  present  the  general  equations  of  vertical 
stress  induction  for  the  two-  and  three-dimensional  Boussinesq  and  Westergaard 
theories,  discuss  techniques  for  their  application,  and  present  FORTRAN  sub¬ 
routines  which  can  be  used  in  the  computation  of  vertical  stresses.  These 
subroutines  are  excerpted  from  an  all-purpose  settlement  analysis  program 
designed  and  coded  by  the  authors.  The  program  handles  stratified  soil  and 
nonhomogeneous  loads,  computes  ultimate  and  time-rate  of  consolidations,  pro¬ 
vides  graphics  capability  for  displaying  lines  of  equal  stress  intensity,  and 
allows  the  user  to  select  either  Boussinesq  or  Westergaard  analysis. 

*  J.  G.  Flock  and  M.  E.  Pittman.  1979.  "Vertical  Stress  Induction  and 
Settlement  Analysis  Program",  U.  S.  Army  Engineer  District,  New  Orleans, 
Corps  of  Engineers,  New  Orleans,  Louisiana. 

Bl 


Discussion 


The  Boussinesq  theory  yields  the  following  formula  for  a  point  load,  P, 
acting  on  a  homogeneous  soil: 


and  the  Westergaard  theory  gives 


where 

p2  =  <l-2p)/2(l-p) 

and 

P  is  the  Poisson's  ratio  for  Che  soil. 

If  we  substitute  Z'  for  pZ, 


For  the  two-dimensional  case,  these  equations  must  be  integrated  from 
y«_tx>toy«  +  °°  to  achieve  the  equations  for  a  distributed  line  load  of 

B2 


infinite  extent.  This  has  been  done  by  many  authors  (at  least  for  the 
Boussinesq  theory)  and  the  following  equations  are  obtained: 
Boussinesq,  distributed  line  load: 


,  2PZ 

1  ir(X2  +  Z2)2 


Westergaard,  distributed  line  load: 


Ao 


PZ' 


Z  IT  (X2  +  z’2) 


Also,  for  the  two-dimensional  case,  the  most  general  loading  condition 
may  appear  as: 


In  this  sketch,  the  load  is  not  homogeneous  nor  geometrically  "visa,"  and  .  r,e 
virgin  ground  is  not  horizontal  since  there  was  an  excavation  prior  to 
subsequent  loading.  The  design  engineer  would  either  have  to  approximate  the 
load  or  prepare  to  input  the  load  as  a  large  set  of  finite  elements. 

Consider  the  natural  ground  prior  to  subsequent  loading  activity,  which 
includes  the  excavation  as  a  loading  activity.  From  some  datum  below  the  natu¬ 
ral  ground  (conveniently,  the  bottom  of  the  proposed  excavation),  if  one  were 
to  remove  all  the  above  material  (assign  negative  stress  values  to  the  soil 
continuum)  and  add  all  but  the  excavated  material  back,  the  natural  ground 
appears  as  a  homogeneous  irregular  area  on  a  horizontal  surface  (datum'* .  Each 
subsequent  load,  using  judicial  superpositioning  techniques,  appears  the  same 
way  because  each  layer  shows  a  common  region  with  the  layer  immediately  below 
it.  Only  the  effective  density  of  the  load  is  altered.  Thus,  superposition 
techniques  can  be  applied  vertically  as  well  as  horizontally,  which  esta¬ 
blishes  the  game  plan  for  general  analysis. 
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If  there  are  i  layered  loads,  each  piecewise-linearly  defined,  having  ji 
homogeneous  areas  between  adjacent  coordinates  defining  the  region,  then 
Ao  is  the  stress  contribution  by  this  elemental  region  at  a  point 

i.ji 

(X  ,Z)  in  the  soil  continuum.  The  total  vertical  stress,  Act  ,  is  therefore 
c  Z 

given  by: 


A0Z  = 


1  Ji 

EE 

n=l  m=l 


i.Ji 


Further  discussion  will  develop  the  equations  for  Aa  for  the  most 

general  case.  For  facility,  the  symbol  will  be  used  for  the  stress  contri¬ 
bution  by  an  elemental  region. 

The  elemental  region  is  defined  to  be  a  homogeneous  region  between 
adjacent  coordinates  defining  a  load  layer.  Since  the  element  extends  down  to 
the  datum,  it  shows  areas  with  previously  defined  regions  of  different 
densities;  but  the  difference  in  areas  and  densities  can  be  taken  care  of  by 
application  technique.  The  elemental  region  extends  infinitely  through  the 
page  and  appears  in  the  sketch  below. 


P  is  a  linear  function  of  X,  given  by 


y  X  -  X  y.  -  y 

P  =  v  01  L-2  +  _I _ v 

y  X,  -  X  X.  -  X  x 

1  o  1  o 


where  y  is  the  unit  weight  of  the  soil. 

Writing  the  equations  for  a  line  load  with  the  substitution  for  P  and 
using  the  elemental  figure,  we  get  the  following  formulas  for  the  line  loads. 
Boussinesq  line  load: 


2  3 

dN 

Z  TT 


Westergaard  line  load: 


dNz  = 


t1-  »  spr  x][<v>2  * z’2] 


Integrating  both  equations  with  respect  to  X  between  Xq  and  X^,  yields  the 
equations  for  the  elemental  trapezoidal  region.  Thus,  we  finally  achieve  the 
following  two  equations  for  a  general  two-dimensional  homogeneous  region  of 
trapezoidal  shape. 

— ou-ss-n-s-a  f  (X,X  -  X  X  -  X  X.  +  X  2  +  z2) 

N  =  _ ll _  Y  +Y  -  Y  -  1-° _ £-2 _ _ £ _ 

Z  r(XrXo)  1  o  1  (  )2  +  z2 

1  c 

Y  (XlXo  "  XcXo  ~  XcXl  +  Xc2  +  z2) 

O  -  -  -  + 

(X  -X.T  +  Z^ 

O  1 


y  (X.-X  )  -  y.(X  -X  )  <  -1  X..-X 

7o  1  c  ■'l  o  c  vtan  1  c  tan 


-1  X  -X  ) 
_o _ c  l 

Z  I 


Westergaard 


y  X.  ■  y.X  y.  -  y  X  .  -1  X.  -  X  .-IX  X 

N  vz  o  1  1  o  J 1  o  c  tan  1  c  -  tan  o  c 

Z  “  TT  2 (X. -X  )  X.  -  X^  Z'  Z’ 

1  o  1  o 

*  <l,(Vxc>2  *  z’2>  -  “«vxc>2  ♦  z'2) 

1  o 


The  following  two  FORTRAN  subroutines  are  provided  for  computing  values 
of  vertical  stress  (SIG)  for  the  Boussinesq  (subroutine  B0US2)  or  Westergaard 
(subroutine  WEST2)  theories  for  the  two-dimensional  case.  It  is  assumed  that 
a  mainline  will  be  written  by  the  users  to  supply  the  calling  arguments  whose 
definition  are  given  preceding  the  listings. 

It  should  be  noted  that  in  the  event  of  an  excavation,  the  amount  of  mate¬ 
rial  above  the  datum  (normally  selected  to  be  the  bottom  of  the  excavation)  is 
relieved  from  the  stress  values  by  a  value  equaling  the  height  of  material  in 
feet  times  its  effective  density.  This  is  a  calling  argument.  The  profile 
above  the  datum  is  then  considered  as  a  load  placed  back  upon  the  soil  as  a 
load  (minus  the  excavated  region).  Also,  for  excavations,  the  X-coordinates 
will  tend  to  infinity,  represented  as  9999  and  given  the  appropriate  sign. 

The  equations  have  been  computed  for  these  limiting  conditions  and  are  built 
into  the  subroutines. 


Calling  Argument  Definitions 


X  -  abscissa  value  (feet)  defining  point  in  soil  continuum  at  which 
vertical  stress  is  to  be  computed. 

Z  -  depth  in  feet  from  the  datum  at  which  stress  is  to  be  computed. 

XL(i,j)  -  abscissa  (feet)  value  for  the  jth  coordinate  defining  the  ith  load 
strata . 

YL(i,j)  -  ordinate  value,  defined  as  above. 

3 

GAM(i)  -  effective  density  (lb/ft  )  of  the  ith  load  strata. 

LL  -  total  number  of  layered  loads. 

KS(i)  -  number  of  coordinates  defining  ith  load. 

3 

SIGU  -  value  of  stress  (lb/ft  )  to  be  relieved  from  calculated  value 
due  to  excavation-type  loading. 

3 

SIG  -  returned  value  of  stress  (lb/ft  )  computed  at  (X,Z). 

UMU  -  Poisson's  ratio  for  the  soil;  Westergaard  routine  only. 

The  subroutines  assume  a  mainline  exists  which  has  read  the  input  into 
the  argument  arrays. 

When  using  these  subroutines,  the  load  layers  should  be  entered  into  the 
coordinate  arrays  defining  the  load  with  the  topmost  layer  first.  Otherwise, 
the  equations,  which  account  for  the  effects  of  layering  on  the  effective  den¬ 
sity,  will  not  yield  the  proper  results.  The  soil  continuum  may  be  strati¬ 
fied;  this  presents  no  difficulty  with  the  Boussinesq  equations  since  soil 


properties  are  not  a  factor,  but  with  the  Westergaard  analysis  the  mainline 

program  should  provide  bookkeeping  to  account  for  changes  in  the  Poisson's 

ratio.  The  mainline  and  subroutines  are  small  enough  to  fit  on  even  very 

small  processors,  if  only  the  values  of  vertical  stress  are  desired.  However 

large  programs  can  be  written  using  these  routines  to  do  a  multitude  of  compu 

tation  and  bookkeeping  to  do  most  desirable  end  products  for  settlement  analy 

sis  and  stress  induction  activity. 

SUBROUTINE  B0US2  (X,  Z,  XL,  YL,  GAM,  LL,  KS ,  SIGU,  SIG) 

DIMENSION  XL  (20,  50),  YL  (20,  50),  GAM  (20),  KS  (20) 

2-DIMENS 10NAL  BOUSSINESQ 
GAML  =0.0 
SI  =  0.0 

D0  200  I  =  1,  LL 

L00P  PERF0RMED  F0R  EACH  L0AD  LAYER 
SIG  =  0.0 
KM  =  KS (I) 

D0  100  J  =  2,  KM 

L00P  PERF0RMED  F0R  EACH  LAYER  ELEMENT 
CH  =  ABS  (XL  (I,  J)  -  XL  (I,  J-l)) 

IF  (CH.  LT.  0.001)  G0  T0  100 
IF  (XL (I ,  J-l).  LE.  -9998.)  G0  T0  50 
IF  (XL (I,  J).  GE.  9998)  G0  T0  80 
NEITHER  X0  0R  XI  G0  T0  INFINITY 

A  =  XL(I ,  J)  *  XL (I,  J-l)  -  X*  XL(I ,  J-I)  -  X*XL(I ,  J)  +  X**2  +  Z**2 

B  =  YL(I ,  J)  *  A/ ( (XL (I ,  J)  -  X)**2  +  Z**2) 

C  =  YL(I ,  J-l)  *A/ ( (XL (I ,  J-l)  -  X)  **2  +  Z**2) 

D  =  ATAN2  ( (XL (I ,  J)-X),  Z)  -  ATAN2((XL(I,  J-l)-X)),  Z) 

E  =  (YL(I ,  J-l)  *  ( XL ( I ,  J)  -  X)  -  YL  (I,  J)  *  (XL(I,  J-l)-X))/Z 

F  =  YL(I,  J-l)  +  YL(I ,  J) 

G  =  Z/(3. 14159  +  ( XL ( I ,  J)  -  XL(I,  J-l))) 

G  =  G/ (GAM(I)  -  GAML) 

SIG  =  G* (F-B-C+E*b)  +  SIG 
G 0  T0  100 

LEFT  X  G0ES  T0  MINUS  INFINITY 
50  A  =  (XL ( I ,  J)  -  X) 

B  =  YL(I ,  J) 

C  =  AT AN 2  (A,  Z)  +  1.5708 
A  =  -A 
G0  T0  90 

RIGHT  X  G0ES  T0  INFINITY 
80  A  =  (XL (I ,  J-l)  -  X) 

B  =  YL(I ,  J-l) 

C  =  1.5708  -  ATAN2  (A,  Z) 

90  G  =  Z/3. 14159  *  (B/Z*C  -  B*A/(A**2  +  Z**2)) 

G  =  G*  (GAM(I)  -  GAML) 

SIG  =  SIG  +  G 
100  C0NTINUE 

GAML  =  GAM(I) 

SI  +  SI  +  SIG 
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200  C0NTINUE 

SUBTRACT  RELIEVED  PRESSURE  FR0M  C0MPUTED 

SIG  =  SI  -  SIGU 

RETURN 

END 

SUBR0UTINE  WEST2  (X,  Z,  XL,  YL,  GAM,  LL,  KS,  SIGU,  UMU,  SIG) 

DIMENSION  XL  (20,  50),  YL  (20,  50),  GAM  (20),  KS  (20) 

C  2-DIMENSIONAL  WESTERGAARD 

GAML  =  0.0 

ETA  =  SQRT  ((1.  -  2 . *UMU) / (2 . * ( 1 .  -  UMU))) 

SI  =  0.0 

D0  100  I  =  1,  LL 

L00P  PERF0RMED  F0R  EACH  L0AD  LAYER 
SIG  =  0.0 
KM  =  KS(I) 

D0  100  J  =  2,  KM 

C  L00P  PERF0RMED  F0R  EACH  ELEMENTAL  REGI0N 

IF  (XL(I,  J-l).  LE.  -9998.)  G0  T0  80 
IF  (XL  (I,J)  .GE.  9998.)  G0  T0  50 
C  IF  HERE,  REGI0N  FINITE 

IF  (ABS  (XL(I ,  J)  -  XL(I ,  J-l)).  LT.  0.001)  G0  T0  100 

A  =  (YL(I ,  J-1)*XL(I,  J)  -  YL(I,  J)  *  XL (I,  J-1))/(XL(I,  J)  -  XL(I,  Jl)) 

B  =  (YL(I,  J)  -  YL(I ,  J-1))/(XL(I,  J)  -  XL(I ,  J-l)) 

C  =  (A  +  B*X) / (Z*ETA) 

D  =  ATAN2  ( (XL(I ,  J)  -X),  (Z*ETA))  -  ATAN2  ((XL(I,  J-l)  -X),  (Z*ETA) ) 

E  =  B/2 .*  (AL0G( (XL(I,  J)  -  X)**2  +  (Z*ETA)**2)  -  AL0G((XL(I,  J-l)  -  X)**2 
$  +  (Z  *ETA) **2) ) 

W  =  (GAM(I)  -  GAML)  *ETA/3. 14159 
SIG  =  W  *  Z*(C  *D  +  E)  +  SIG 
G0  T0  100 

C  LEFT  X  G0ES  T0  INFINITY 

50  SIG  =  (GAM (I)  -  GAML)*YL(I ,  J-l)/3. 14I59*(1 .570796  -  ATAN2((XL(I,  J-l)  -X) , 
$  (Z  *  ETA)))  +  SIG 
G0  T0  100 

C  RIGHT  X  G0ES  T0  INFINITY 

80  SIG  =  (GAM(I)  -  GAML)  *  YL(I,J)/3. 14159  *(1.570796  +  ATAN"' ( (XL(I ,  J)  -X), 

$  (Z  *ETA) ) )  +  SIG 
100  C0NTINUE 

GAML  =  GAM(I) 

SI  -  SI  +  SIG 
200  C0NTINUE 

SIG  =  SI  -  SIGU 

RETURN 

END 

For  the  three-dimensional  case,  the  starting  point  is  the  equation  for  a 
point  load  vertical  to  a  horizontal  surface,  and  again  superposition  will  be 
used.  The  load  need  not  be  homogeneous,  although  the  analysis  will  assume 
that  it  is,  simply  since  a  non-homogeneous  load  can  be  separated  into  homogen¬ 
eous  regions  and  superpositioning  will  allow  summing  of  the  regions.  Each 
homogeneous  region  is  assumed  to  be  definable  by  (X,Y)  coordinate  pairs  in 


plan  view;  i.e.,  the  region  boundaries  are  piecewise  linear.  The  most  general 
case  is  therefore  a  homogeneous  polygon  as  a  horizontal  surface  as  shown  by 
the  sketch  below: 


To  analyze  this  case,  consider  the  plan  view  to  uncover  the  elemental  region. 


Using  superpositioning,  the  region  can  be  analyzed  by  considering  each  region 

defined  by  (X.,Y.),  (X.  .,Y.  ,),  (X  ,Y  )  incrementally.  Thus,  the  incremental 
11  j+lj+l  cc 

region  to  be  considered  is  further  reduced  to  triangular  regions.  In  fact. 
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since  any  triangular  region  can  be  considered  as  a  combination  of  two  right 
triangles,  the  right  triangle  becomes  the  most  elemental  region  for  formulae 
derivations. 


Hence,  the  three-dimensional  equation  will  be  developed  using  a  right- 
triangular,  homogeneous  load  on  a  horizontal  surface  for  a  point  immediately 
beneath  one  corner  of  the  load  and  at  a  depth  Z. 


The  total  vertical  stress  is  the  algebraic  sum  of  contributions  by  all  such 
triangular  regions,  summed  over  all  regions  of  homogeneous  density. 

For  the  three-dimensional  Boussinesq  case,  the  following  equation  holds. 
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For  the  three-dimensional  Westergaard  case,  the  following  equation  holds. 
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In  the  Westergaard  case,  if  we  replace  yZ  by  Z',  as  was  previously  done 
in  the  two-dimensional  problem,  the  equation  simplifies. 

„  P  ,a  ,  ,b  x/a  ,,X^+y^  .  ,N3/2  , 

N_  =  - 7  /  dx  /  ( - +1)  dy 
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